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Preface 


The theory of numbers is the last great uncivilized continent of 
mathematics. It is split up into innumerable countries, fertile 
enough in themselves, but all the more or less indifferent to one 
another’s welfare and without a vestige of a central, intelligent 
government. If any young Alexander is weeping for a new world to 
conquer, it lies before him. 


Eric Temple Bell 


The elementary theory of numbers should be one of the very best 
subjects for early mathematical instruction. It demands very little 
previous knowledge; its subject matter is tangible and familiar; the 
processes of reasoning which it employs are simple, general, and 
few; and it is unique among the mathematical sciences in its appeal 
to natural human curiosity. A month’s intelligent instruction in the 
theory of numbers ought to be twice as instructive, twice as useful, 
and at least ten times as entertaining as the same amount of 
“calculus for engineers.” 


Godfrey Harold Hardy 


Number theorists are like lotus-eaters—having once tasted of this 
food they can never give it up. 


Leopold Kronecker 


Some time ago | was looking at several textbooks for the undergrad- 
uate number theory course. I was struck by how few illustrations were 
included in many of those textbooks. A number—specifically a positive 
integer—can represent many things: the cardinality ofa set; the length 
of a line segment; or the area of a plane region. Such representations 
naturally lead to a variety of visual arguments for topics in elementary 
number theory. Since the number theory course usually begins with 
properties of the positive integers, the texts should have more pictures. 
That observation became the motivation for this book. 

Work on this book began when I was invited to give a talk at the 
MAA’s MathFest in Albuquerque in August 2005, in a session entitled 
“Gems of Number Theory” organized by Arthur Benjamin and Ezra 


ix 
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Brown. The title of that talk was “Some visual gems from elementary 
number theory.” Later a version of the talk appeared as an article in the 
February 2008 issue of Math Horizons. This article was subsequently 
included in the 2009 MAA book Biscuits of Number Theory, edited by 
Benjamin and Brown. 

Nuggets of Number Theory: A Visual Approach is not a textbook. Al- 
though it is designed to be used by the instructor of an undergraduate 
number theory course as a supplement to a standard textbook, it will 
be of interest to anyone who loves number theory. Certain chapters in 
the book may be of interest to those who teach discrete mathematics, 
abstract algebra, and teacher preparation courses. The book can also 
be used as a resource for group projects or extra-credit assignments. 

A nugget is a lump of precious metal or, more generally, something 
of significance or of great value. The nuggets in this book are topics 
in number theory for which I believe a visual approach is appropriate 
and beneficial, with a chapter devoted to each one. Chapter |1] is de- 
voted to figurate numbers—numbers that can be represented by ob- 
jects such as pebbles arranged in geometrical patterns—that were stud- 
ied by the early Greeks. Chapter (2) deals with the important concept of 
congruence, and includes visual demonstrations of Fermat's little theo- 
rem and Wilson’s theorem. Visual approaches to the solutions to linear 
Diophantine equations and Pell equations are the subject of Chapter Bl 
Pythagorean triples—solutions to the Diophantine equation a? + b? = 
c*—are represented naturally by integer-sided right triangles and ex- 
plored in Chapter lal Certain irrational numbers and their rational ap- 
proximations, including ones based on continued fractions, appear in 
Chapter |5, Many identities for the Fibonacci and Lucas numbers are 
presented visually in Chapter i and similarly for properties of perfect 
numbers in Chapter iA Each chapter includes a set of exercises, with 
solutions to all the exercises following the final chapter. The book con- 
cludes with references and a complete index. 

Acknowledgments. I would like to express my appreciation and grat- 
itude to Susan Staples and the members of the editorial board of Class- 
room Resource Materials for their careful reading of an earlier draft of 
this book, and for their many helpful suggestions. I would also like to 
thank Stephen Kennedy, Carol Baxter, Beverly Ruedi, and Bonnie Ponce 
of the MAA publication staff and Christine Thivierge, Sergei Gelfand, 
Courtney Rose, and Jennifer Wright Sharp of the AMS production staff 
for their encouragement, expertise, and hard work in preparing this 
book for publication. 


Roger B. Nelsen 
Lewis & Clark College 
Portland, Oregon 


CHAPTER 1 


Figurate Numbers 


Why are numbers beautiful? It’s like asking why is 
Ludwig van Beethoven’s Ninth Symphony beautiful. 

If you don’t see why, someone can’t tell you. I know 
numbers are beautiful. If they aren’t beautiful, nothing is. 


Paul Erdés 


Mighty are numbers, joined with art resistless. 


Euripides 


The figurate numbers are positive integers that can be represented 
geometrically by arrangements of points, or physically by arrangements 
of objects like pebbles. As such they form a connection between geom- 
etry and number theory—indeed, number words such as square and 
cube reflect this relationship. Figurate numbers appear in some of the 
works by early Greek geometers, and today they regularly appear in un- 
dergraduate number theory texts. In fact, figurate numbers appear in 
every chapter in this book. Our approach to the figurate numbers re- 
flects their origin in geometry. 

We begin with plane figurate numbers known as polygonal num- 
bers, represented by arrangements of objects in the plane in polygonal 
patterns. We conclude this chapter with several three-dimensional fig- 
urate numbers. 


1.1. Polygonal numbers 


The triangular numbers are sums of consecutive positive integers, 
e.g., 1, 1+2, 1+2+3, etc., and can be represented by arrangements of 
objects in a triangle, as shown in Figure We use the notation 
T, = 1+2+-::+k for the kth triangular number for k = 1 (and de- 
fine Ty) = 0). The shape of the triangle is arbitrary, usually it is drawn 
as an equilateral, isosceles, or right triangle. Similarly, the squares are, 
of course, represented by objects in a square pattern, as shown in Fig- 


ure 
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(a) Ty: 1, 3, 6, 10, 15,... (b) k?: 1, 4,9, 16, 25,... 


FIGURE 1.1 


Observe in Figurel.1(b)\ that k? = 1+3+5 +--+ (2k —1), that 
is, k? is the sum of the first k odd positive integers, analogous to T, = 
1+2+-::-+k for the triangular numbers. 


Raphael and the fourth triangular number 


In his classic painting The School of Athens, the Italian renaissance 
painter Raphael (1483-1520) included Pythagoras in the lower left cor- 
ner, writing in a book, with a slate at his feet. At the bottom of the slate 
is the tetractys, the representation of T, as ten objects (Roman numerals) 
arranged in a triangular pattern, with the sum X. 


The School of Athens, Pythagoras, and the tetractys 


The drawing of the slate is from the book Descrizione delle immagini 
dipinte da Raffaelle d’Urbino by Giovanni Bellori, published in 1751. 


Carl Friedrich Gauss and the 100th triangular num- 
ber 
Nearly every biography of the great mathematician Carl Friedrich Gauss 


(1777-1855) relates the following story. When Gauss was about ten years 
old, his arithmetic teacher asked the students in class to compute the sum 
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1+2+3+---+100, anticipating this would keep them busy for some time. 
He barely finished stating the problem when young Carl came forward and 
placed his slate on the teacher’s desk, void of calculation, with the correct 
answer: 5050. When asked to explain, Gauss admitted he recognized the 
pattern 1+100 = 101,2+99 = 101,3+98 = 101,andsoonto50+51 = 
101. Since there are fifty such pairs, the sum must be 50-101 = 5050. The 
pattern for the sum (adding the largest number to the smallest, the second 
largest to the second smallest, and so on) is illustrated below, along with 
a portrait of Gauss on a pre-euro 10 Deutsche Mark note. 
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Gauss and his computation 


Two more types of polygonal numbers are the oblong numbers O;, 
(objects in a rectangular array with k rows and k + 1 columns) and the 
pentagonal numbers P, (objects in a pentagonal pattern), as shown in 
Figure 

Clearly O, = k(k +1),and0O, =2+4+6+-: + 2k, that is, O, is 
the sum of the first k even positive integers. 

In Figure we illustrate two simple but important relationships 
involving the triangular, square, and oblong numbers. 

k(k+1) 


Figure [1.3(a)lyields a formula for the triangular numbers T, = —>—, 


while Figure foreshadows a method we shall use to find formulas 
for other polygonal numbers. 
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(b) Py: 1,5, 12, 22,35... 


FIGURE 1.2 


(a) O, = 2T, (b) k* = Ty + Th_1 


FIGURE 1.3 


In Table f1.1}we exhibit the values of the first few polygonal numbers 
based on regular n-gons (polygons with n sides), i.e., the values of the 
kth n-gonal numbers for3 <n < 8and1<k < 10. We refer to k as 
the rank of the n-gonal number. 

Table yields some important information about a general for- 
mula for the kth n-gonal number. The formula should be linear in n, 


TABLE 1.1 

k= |1/2}3)4],5]6] 7 8 9 | 10 
n = 3 (Triangular) |}1/3] 6 |10) 15) 21] 28 | 36 | 45 | 55 
n = 4 (Square) 1|4) 9 |16 > 25)36| 49 | 64 | 81 | 100 
n =5(Pentagonal) |1|)5]12)22)/35/51] 70 | 92 | 117) 145 
n = 6 (Hexagonal) |1/6)]15) 28) 45) 66/ 91 | 120/153 | 190 
n = 7 (Heptagonal) | 1 | 7 | 18 | 34) 55 | 81 | 112 | 148 | 189 | 235 
n= 8 (Octagonal) | 18] 21) 40) 65 | 96/133 | 176 | 225 | 280 
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since the numbers in each column after the first increase linearly, and 
quadratic in k, since the differences of adjacent numbers in the rows 
increase linearly. To find the formula, we “triangulate” the kth n-gonal 
number, partitioning it into a collection of triangular numbers. We il- 
lustrate with the 5th hexagonal number 45 in Figure 


FIGURE 1.4 


Hence the 5th hexagonal number 45 equals 7; + 37,. In general, the 
kth n-gonal number equals 7;, + (n— 3)T;,_, , since it can be partitioned 
into one copy of T; and n — 3 copies of T;_,, one for each of the n — 3 
diagonals emanating from each vertex of the n-gon. Other expressions 
for the kth n-gonal number include k+(n — 2) Ty_ andk?+(n—4)T;_1. 
Each of these formulas is equivalent to =[(n —2)k?-(n—4)k] fork = 
landn > 3. 


Example 1.1. A test for n-gonal numbers. Is 903 a triangular number? 
Is 783 a heptagonal number? A test can be constructed from the for- 
mula at the end of the preceding paragraph. Set the candidate x equal 


to =[(n — 2) k* —(n— 4)k] and complete the square on k. To do so we 
first multiply x by 8(n — 2): 


8(n — 2) x = 4[(n — 2)°k2 — (n—4) (n—2)k]; 
and then add (n — 4)’: 
8(n—2)x+(n—4)* =4 [cn — 2)°k2 —(n— 4) (n—2) I| Gedy 
=[2(n-2)k-(n—4)]’. 
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So if x is an n-gonal number, then 8 (n — 2) x + (n—- 4)° is a perfect 
square. To find its rank k, we add n — 4 to the root 2 (n — 2) k —(n—4) 
and divide by 2 (n — 2). 

Isx = 903 triangular? Heren = 3 so we consider 8x+1 = 8 (903)+ 
1 = 7225 = 857, so 903 is a triangular number. Which one? We add 
—1 to 85 and divide by 2: k = (85 — 1)/2 = 42, hence 903 = Tq. 

Is x = 783 heptagonal? Here n = 7 so we consider 40x + 9 = 
40 (783) + 9 = 31329 = 177’, so 783 is a heptagonal number. Which 
one? We add 3 to 177 and divide by 10: k = (177 + 3)/10 = 18, hence 
783 is the 18th heptagonal number. 


Example 1.2. Many integers appear more than once in Table ft} For 
example, 36 is both triangular and square, every hexagonal number is 
triangular, and extending the table shows that 210 is both triangular 
and pentagonal, and 225 is both square and octagonal. Are these coin- 
cidences or patterns? We shall answer these questions in later sections 
in this chapter and in subsequent chapters. 


Polygonal numbers in number theory 


Polygonal numbers have long played a role in the theory of numbers. 
Their properties were investigated by Nicomachus of Gerasa (circa 60- 
120 CE) and Diophantus of Alexandria (circa 200-284 CE). Pierre de Fer- 
mat (1601-1665) wrote that every positive integer is the sum of three 
or fewer triangular numbers, four or fewer squares, and, in general, n or 
fewer n-gonal numbers. His proof (if it existed) has never been found. Carl 
Friedrich Gauss proved the triangular case and on July 10, 1796 wrote 
in his diary “EYPHKA! Num = A+A+ A.” The case for squares was 
proven in 1770 by Joseph Louis Lagrange (1736-1813) and is known as 
Lagrange’s four-squares theorem. In 1813 Augustin-Louis Cauchy (1789- 
1857) proved the general case of Fermat’s claim. France has honored the 
three Frenchmen on postage stamps—Fermat in 2001, Lagrange in 1958, 
and Cauchy in 1989. 
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1.2. Triangular number identities 


There are a great many identities relating the triangular numbers. 
All can be established with elementary algebra, but many can be illus- 
trated visually. Here are several examples. 


Example 1.3. Consecutive sums of consecutive integers. Consider the 
sequence of positive integers: 
123456789 10 11 12 13 14 15:. 


If one has an unlimited supply of + and = signs and commas, it is easy to 
turn the sequence into a sequence of true arithmetic statements: 
14+2=3, 44+5+6=7+8, 94+104+114+12=134+14+415,.... 
The above sums—3, 15, and 42—are multiples of the first three trian- 
gular numbers, i.e., 37,, 57>, and 7T3. In fact, the identity is 
n?+(n?4+1)4+--4+(n? +n) =(n?+n+1) +--+ (n? + 2n) 
= (2n + 1)T;. 

Figure gives an illustration for n = 4 by considering two arrange- 
ments of nine 7,’s in the form of collections of unit squares. The num- 


bers along the right side of the figure enumerate the number of unit 
squares in that row. 


FIGURE 1.5 


Example 1.4. Summing integer cubes. Sums of integer cubes are related 
to triangular numbers in the formula 
13423 4--4n3 =(1424+--4+n) =7? 


often encountered in precalculus mathematics. Figure f1.1(a)jillustrates 
k =T, —T,-_1, and Figure illustrates k* = T,, + T;,_,. Multiplying 
these together yields 


(1.1) k? =k k? = (Ty — Ty-1) Te + Te-1) = Te — Ty 
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Summing ((1.1)) and observing that the sum telescopes yields the for- 
mula for the sum of the first n positive integer cubes (recall that Tyg = 0): 


n n 
=) @-7p=7-=2 
k=1 k=1 


In Figure (Golomb, 1965]] we see an illustration of this result, where 


we represent k? by k copies of squares with area k? from 1 to n (shown 
here for n = 6). When k is even, two squares overlap, but the area of 
the overlap is the same as the area of the square (in white) not covered 
by the shaded squares. 


FIGURE 1.6 


Example 1.5. Summing squares of consecutive triangular numbers. 
While the difference of the squares of two consecutive triangular num- 
bers is a cube, the sum of the squares of two consecutive triangular 
numbers is triangular, T? + T7_, = Tyz, as illustrated in Figure /1.7| for 
k = 4. We represent T? with T, copies of T, and similarly for T?. 


Example 1.6. Sums of consecutive triangular numbers. Analogous to 
Example |1.3) we consider the sequence of triangular numbers: 


1 3 6 10 15 21 28 36 45 55 
66 78 91 105 120 136 153 171-.. 
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+O-O-D0 O-O-O-O 
>O0-0-D 0C—OoO-O-0 C—O-O-— 


FIGURE 1.7 


If, as in Example we have an unlimited supply of + and = signs and 
commas, we can also turn this sequence into a sequence of true arith- 
metic statements: 


1+3+6=10, 
15+ 214+ 28+ 36=45+55, 
66+ 78+ 914+ 105 + 120 = 1364+ 1534+ 171,.... 


The pattern is 
Th2-n-1 + Th2-n $0 + Tp2-1 = Tn2 + Tr24a + + Tr24n-2: 
For n = 4 the pattern is T,; + 7,2 + 7,3 +Ty4 + M15 = Tig + 717 + Ti, 
which we now illustrate [Nelsen and Unal, 2012]. Figure f-dshows that 
Ty6 +T17 + Tig = Ty5 +T14 +713 +1-47 43-47 45-4. 


But 1-424+3-424+5-42 =9.42 = 12” =7,, + Tz, as illustrated in 
Figure which completes the illustration. 


FIGURE 1.8 
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FIGURE 1.9 


In the following theorem we relate the factors of one triangular num- 
ber to a pair of larger triangular numbers with a triangular sum. The 
theorem is surprisingly useful, as the two examples following the proof 
illustrate. 


Theorem 1.1. Let n, p, and q be positive integers. Then 
(1.2) Ty = pq ifand only if Tr4p + Tn+q = Tn+p+q: 


Proof. In Figure we illustrate T,4y+q forn = 6,p = 3,and q = 7. 
Counting the dots using the inclusion-exclusion principle yields 
Thap+q = Tn+p + Tr+q — Tr + Pa, from which follows. The the- 
orem can also be established algebraically. 


FIGURE 1.10 


Figure shows that T, = 21 = 3: 7 is equivalent to Tg + T,3 = 
T,¢. Similar but perhaps more interesting sums are encountered in the 
following examples. 


Example 1.7. Triangular sums and Pythagorean triples. Let k be a pos- 
itive integer. Then ([1.1)) yields the following triangular sums of two 


1.2. TRIANGULAR NUMBER IDENTITIES 11 


triangular numbers. Let k be a positive integer; then: 
(a) with T2, = k(2k + 1), we have 


T3y + Taaa = T5x413 
(b) with Ty, = k(8k + 2), we have 


Tsk + Ti2k+2 = Ti3k+2: 
(c) with Ten42 = (2k + 1)(9k + 3), we have 


Tgx+3 + Tiskss = Ti7K-+6: 
(d) with Tioe43 = (8k + 2)(9k + 3), we have 


Thox+s + To1x+6 = Trox-+8- 


Each triple of coefficients of k in the ranks of the triangular numbers 
in the sums—(3,4,5), (5,12,13), (8,15,17), and (20, 21,29)—is a 
Pythagorean triple, i.e., a triple (a, b,c) where a? + b? = c?. Observe 
some patterns in the factors of the given triangular number T,, and the 
triples: 


i) If the coefficient of k in one factor of T,, is 1, then the even leg 
and hypotenuse of the triple are consecutive integers. 
ii) If the coefficient of k in one factor of T,, is 2, then the odd leg 
and hypotenuse of the triple differ by 2. 
iii) Ifthe coefficients of k in the factors of T,, are consecutive inte- 
gers, then the legs of the triple are consecutive integers. 


Pythagorean triples are the subject of Chapter lah 


Example 1.8. Square triangular and oblong triangular numbers. Some 
numbers, such as 1 and 36, are both square and triangular; and others, 
such as 6 and 210, are both oblong and triangular. Are there others? 


Setting g = p in ({1.2) yields 
Ty = p? if and only if Th42p = 2Tntp = (n+ p)(n+p+1), 


that is, each square triangular number corresponds to a larger oblong 
triangular number. Setting q = p + 1 in (1.2) yields 


. 2 
Ty = p(p + 1) if and only if Ths2p+1 = Tnap + Tnepsi = (n+p +1), 
that is, each oblong triangular number corresponds to a larger square 
triangular number. Combining the two equivalences in two ways yields 
(1.3) Ty = p? ifand only if Tsn44p41 = (2n + 3p + 1)” 

and 


(1.4) 
T, = p(p + 1) ifand only if T3n44543 = (2n + 3p + 2)(2n + 3p + 3). 
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Equivalence generates infinitely many square triangular numbers 
T, = 12,T3 = 6%, Tig = 35°, Toag = 204”, etc, and generates 
infinitely many oblong triangular numbers T3 = 2 - 3, To9 = 14: 15, 
Ty19 = 84: 85, Tgog = 492 - 493, etc. We shall continue our study of 
multi-polygonal numbers (numbers polygonal in two or more ways) in 
the next two chapters. 


You may have noticed that the nth triangular number is a binomial 


coefficient, i.e., T, = >), as illustrated with a portion of Pascal’s tri- 
angle in Figure anak One explanation for this is that each is equal to 
n(n + 1)/2, but this answer sheds little light on why it is true. Here is a 
better explanation. 


| 2 4 - 
(3 = 
14641 
1 5 1010 5 1 
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1 7 21 3535217 1 re 3 n+l 


FIGURE 1.11 


Theorem 1.2. There exists a one-to-one correspondence between a set of 
T,, objects and the set of two-elements subsets of a set with n + 1 objects. 


Proof. See Figure /1.11(b) [[Larson, 1985], and recall that the binomial 


coefficient (*) is the number of ways to choose 2 elements from a set of 
k objects. The arrows denote the correspondence between an element 
of a set with T,, objects and a pair of elements from a set with n + 1 
objects. 


We conclude this section with a recurrence relation for the triangu- 
lar numbers that generalizes to a recurrence relation for all the n-gonal 
numbers [[Edgar, 2017|]. We begin with a general inclusion-exclusion 
principle for triangular numbers. 

Theorem 1.3. [f0 < a,b,c <nanda+b+c 2 2n, then 
(1.5) Ty = Ta + Ty + Te — Ta+p-n — Th+c-n — Te+a-n + Ta+b+c-2n: 
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FIGURE 1.12 


Proof. See Figure and compute 7, using the inclusion-exclusion 
principle. 


Corollary 1.4. Forallk = 2, T,., = 3T;, — 3T,_1 + Ty_-2. 


Proof. Set (n,a,b,c) = (k + 1,k,k,k) in (11.5), noting thata+b+c= 
3k > 2(k+1) = 2nsincek = 2. 


With initial conditions Tj = 0, T,; = 1, and T, = 3 the identity in 
Corollary produces all the triangular numbers. Now let pz denote 
the kth n-gonal number. In the preceding section we saw that every n- 
gonal number is a sum of triangular numbers, ie., pe = T,+(n—3)Tx-1, 
which yields the following corollary. 


Corollary 1.5. For allk = 2 andn = 3, pei, = 3pk — 3pR_-1 + De-2- 


With initial conditions p>} = 0, pf = 1, and p} = n this identity 
generates all the n-gonal numbers. 


Example 1.9. To compute the heptagonal numbers p; we use the iden- 
tity in the corollary with pj = 0, p/7 = 1,and p3 = 7: 


pl =3-7-3-140=18, 
pl =3-18—-3-741= 34, 
pe = 3-34-—3-18+7=55, etc. 


1.3. Oblong numbers and the infinitude of primes 


The earliest proof that there are infinitely many prime numbers is 
probably Euclid’s in the Elements (Book IX, Proposition 20). After more 
than 2000 years, it is difficult to find a better one. Here is a version of 
Euclid’s proof employing the dimensions n and n + 1 of the nth oblong 
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number O,, created by Ernst Eduard Kummer (1810-1893) in 1873 
[Ribenboim, 1996]. In the proof we show that the dimensions n and 
n + 1 of O,, have no common prime factors. Kummer’s proof is indirect 
like virtually all proofs of the infinitude of primes, since one definition 
of “infinitely many” is “not finitely many.” 


Theorem 1.6. There are infinitely many prime numbers. 


Proof. Assume there are only k primes, p; = 2, pz = 3, p3 = 5,..., De: 
Let N = pyp2p3 °°: Px. Since N + 1 is larger than p, it is not prime, and 
thus N + 1 has a prime divisor p; in common with N. Since p; divides 
both N and N +1, it divides (V + 1)—N = 1, whichis impossible. Hence 
there are infinitely many primes. 


Euclid primes 

Numbers of the form N; = p1P2P3 °** Dx are called primorials (from prime 
and factorial), and the number E, = N; +1 is called a Euclid number. The 
first five Euclid numbers 3, 7, 31, 211, 2311 are prime (and called Euclid 
primes), however E, = 30031 = 59: 509. It is not known if the number 
of Euclid primes is finite or infinite. 


1.4. Pentagonal and other figurate numbers 


The kth pentagonal number P, can be computed by T;, +2T;,_; and by 
k? + T,_1. There are two other less obvious but useful ways, illustrated 


in Figure 1.13} 


(a) 


FIGURE 1.13 


Distorting the pentagonal shape of P, into the trapezoid in Figure 
shows that P, = Tox, — Te-1, and Figure [1.13(b)} shows how 


three copies of the trapezoid yield a triangle, so that T3,_, = 3F, and 
hence P, = T3,_;/3 (in the next chapter we show that T3;_, is always 
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divisible by 3). We shall use this result in Chapter B to find infinitely 
many triangular pentagonal numbers. 

Examination of Table |1.1) shows that every hexagonal number H,, is 
a triangular number with an odd rank, and that their common value is 
the product of their ranks, i.e., Hy, = Top_1 = k(2k —1). See Figure 
for an illustration when k = 5, i.e, H; = Tg = 5-9, on the cover of the 
January 2014 issue of The College Mathematics Journal. 
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FIGURE 1.14. (This illustration originally appeared on the 
cover of the January 2014 issue of The College Mathemat- 
ics Journal. © The Mathematical Association of America, 
2014. All rights reserved.) 


Analogous to our expression of a pentagonal number as the differ- 
ence of two triangular numbers, each octagonal number is the differ- 
ence of two squares. In Figure [1.15] we illustrate Oct, = k (3k —2) = 
(2k —1)* —(k—1)° fork = 4. 


FIGURE 1.15 


16 1. FIGURATE NUMBERS 


Chinese checkers and star numbers 


Chinese checkers is a board game usually played with marbles on a 
wooden board with holes arranged in a star-like pattern, as illustrated be- 
low on the left. The standard board has S; = 121 holes, where S,, is the 
nth star number. For n 2 2, S, = 12Tp_; +1 = 6n(n —1) + 1, seen in 
the middle image below. The second star number S, = 13 appears in the 
Great Seal of the United States on the green side of the one-dollar bill, as 
shown below on the right 


Chinese checkers and the star numbers S; and S, 


In spite of its name, Chinese checkers did not originate in China and is 
unrelated to ordinary checkers. The game was invented in Germany in the 
1890s and first produced in the U.S. in 1928. The name “Chinese checkers” 
was chosen for marketing purposes. 


1.5. Polite numbers 


An integer is polite if it can be written as the sum of two or more con- 
secutive positive integers in at least one way. For example, every trian- 
gular number (greater than 1) is polite, and every pentagonal number 
(greater than 1) is polite (see Figure f1.5). Some integers have more than 
one polite representation, for example, 15 has three polite representa- 
tions: 7+8,4+5+6,and1+2+3+4+5. Numbers that can be written 
as the sum of two or more consecutive positive integers greater than 1 
are called trapezoidal, e.g., all pentagonal numbers greater than 1 are 
trapezoidal, and 15 is both triangular and trapezoidal. So polite num- 
bers are triangular or trapezoidal (or both). So we ask—which integers 
are polite? 


Theorem 1.7. A positive integer n is polite if and only if it is not a power 


of 2. 


Proof. Assume n is nota power of 2, i.e.,n = 2°(2b + 1), where a = 0 
and b > 1,andletm = min (Zee 2b+ 1} and M = max (20% 2b+ 1}. 
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FIGURE 1.16. (This illustration originally appeared in 
C. L. Frenzen, “Proof without words: Sums of consecu- 
tive positive integers”, Mathematics Magazine 70 (1997), 
294. © The Mathematical Association of America, 1997. 
All rights reserved.) 


Thus 2n = mM, which we represent by 27 balls in an m X M array, as 
shown in Figure 1.16 [Frenzen, 1997]. 


M+m-1 M-m+1 M+m-1 
Then M > 3,m = 2, > 2, > 1, and 
M-m+1 


2 2 


> 1so that 


fam rd in cma reer M+m-1 

aa 2 2 2 
Now assume that n is polite, thatis,n = (a+1)+(a+2)+--+(at+k) 
witha => Oandk = 2,sothat 2n = k(k+2a+1). Sincek andk+2a+1 


have opposite parity, one of them is odd and hence n cannot be a power 
of 2. 


When 2b + 1 is a composite odd number in the above proof, n = 
2°(2b + 1) has polite representations different from the one in Fig- 
ure For example, the figure illustrates 18 = 3+4+5+6. Butif we 
partition a set of 36 balls ina 3 x 12 array, we have 18 = 5+6+7. When 
2b + 1 is prime, 7 has only the polite representation in Figure e.g., 
when n = 14, the only polite representation is 14 =2+3+4+5. 


Example 1.10. Non-trapezoidal polite numbers. If n = 2°(2b + 1), 
where a = 0 and b = 1, is a non-trapezoidal polite number, then nis a 
triangular number with 2b+1 prime. In this case the balls in Figure [1.16] 
represent an oblong number partitioned into two triangular numbers, 
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as in Figure f1.3(a)} There are only two forms for such triangular num- 
bers, Tox = 2*~1(2* + 1) with 2* + 1 prime, and T,x_, = 2*~1(2* — 1) 
with 2* — 1 prime . When 2* + 1 is prime, k is a 
power of 2 and 2* + 1 is a Fermat prime. When 2* — 1 is prime, so is k 
(and 2* — 1 is known as a Mersenne prime) and T>x_, is an even perfect 
number, the subject of Chapter|7} As a consequence, the non-trapezoidal 


polite numbers are rather rare, the only ones less than 2,000,000,000 
are 


3,6, 10, 28, 136, 496, 8128, 32896, and 33550336. 


1.6. Three-dimensional figurate numbers 


Figure shows how cannonballs were stored at Fort Monroe in 
Hampton Roads, Virginia, a Union military installation in 1861. The 
piles of cannonballs had both triangular and rectangular bases. 


FIGURE 1.17 


Example 1.11. Tetrahedral numbers. How many cannonballs are there 
in a tetrahedral pile? In Figure we see cannonballs in a pile 
whose shape is that of a regular tetrahedron. The cannonballs in each 
layer are arranged in a triangle, and the number of balls in each layer is 
a triangular number. The number of cannonballs in such a pile with n 
layers is the tetrahedral number Tet,. Figure [1.18(b)| shows that Tet; = 
1+3+6+10+15 = 35. 

In Figure we see aT, X (n + 2) rectangle of unit 
squares partitioned into three copies of T; + Tz + -:- + Tn, and hence 
3(T, + T2 +--+ + T,) = T,(n + 2), and thus 


n 
n(n+1)(n+2 
Tet, =) 7, = StF) 
k=1 
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FIGURE 1.19 


Note that like T,,, Tet, is also a binomial coefficient: Tet,, = ("*?). 


Example 1.12. Pyramidal numbers. Analogous to the tetrahedral num- 
bers we have the pyramidal (or square pyramidal) numbers, represented 
by the number of cannonballs stacked in a pyramidal pile with a square 


base, as seen in Figure 


© 
<p LP 
(a) (b) 
FIGURE 1.20 


The number of cannonballs in each layer is a square, and so the pyra- 
midal number Pyr,, is the sum of the first n squares. Figure [1.20(b 


shows that Pyr, = 1+4+9+16+25 =55. 

In Figure we see arepresentation of Pyr, with unit cubes rather 
than cannonballs. The white cubes represent Tet; while the gray cubes 
represent Tet,, so that Pyr, = Tet, + Tets. 
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FIGURE 1.21 


Hence in general we have 


n 
Pyr,, = >. ke 
k=1 


= Tet,_; + Tet, 
(n—1)n(n+1) n(n+1)(n+2) 
—— + a 
6 6 
_ n(n + 1)(2n + 1) 
= ; : 
For another derivation of the formula for Pyr,, as a sum of squares, see 
Figure (Gardner, 1973; |Kalman, 1991]. Three sets of the first n 


squares (for n = 5) are shown in Figure 
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FIGURE 1.22 
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If we slice up the squares in the third set as indicated by the different 
shadings, all the unit squares can be rearranged into the rectangle in 
Figure |1.22(b)|with basen+1+n = 2n+1 and height 1+2+--+n =T,, 
so that 
n(n + 1)(2n + 1) 

, 
or 3 Pyr,, = (2n + 1) T,, from which the formula for Pyr,, follows upon 
division by 3. 


3(17 +27 +--+n7) =(2n+1)T, = 


Example 1.13. Sums of pentagonal numbers. For our final example of 
three-dimensional figurate numbers, we sum pentagonal numbers. In 
Figure [1.23(a)] we illustrate P, = k? + T,_1 witha vertical arrangement 
of k* unit cubes joined to a horizontal arrangement of T,_, unit cubes 
for k = 1, 2,3, and 4 (= n). The resulting configurations of cubes can 


now be nested to form a pile with base area T,, and height n as shown 
in Figure |1.23(b)j [Miller, 1993] so that 


n 

n?(n +1) 
> Pe = Tn = ——. 
k=1 
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FIGURE 1.23 


For a comprehensive overview of figurate numbers and their prop- 
erties, see [[Deza and Deza, 2012]. 


1.7. Exercises 


1.1 In Figure 1.3] we saw illustrations of two relationships between tri- 
angular, oblong, and square numbers. Two more relationships are 
k2 + Oy = To, and (k +1)? + Og = T2441. Can you illustrate these 
with figures similar to Figure {1.3} 

1.2 Show that the sum of the first n (n = 0) powers of 9 is a triangular 
number, e.g., 1 +9 = T,,1 +9 + 81 = 743, etc. 

1.3 One triangular number can be twice another, as in T3 = 2T>. Are 
there other instances of this phenomenon? 
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1.4 In Section we saw that each pentagonal number is the differ- 
ence of two triangular numbers. Prove that the same is true of odd 
squares, e.g., 5% = T, —T>, 77 = Tyg —Ts, etc. Can you illustrate the 
result with a picture? 

1.5 In Table ft. djit appears that the only prime n-gonal numbers are the 
2nd p-gonal numbers for p a prime. Prove that this is indeed the 
case. 

1.6 Centered triangular numbers. The centered triangular number cy, 
enumerates the number of dots in an array with one central dot 
surrounded by dots in triangular borders. The number cy = 46 is 
illustrated in Figure 


FIGURE 1.24 


In Figure we see that c, = 1+ 37,. Show that another 
formula for Cy is Ch = Ty_-1 + Tn + Tn forn = 2. 

1.7 Show that the sum of the first n odd squares is a tetrahedral num- 
ber. [Hint: See Figure where we represent integers by unit 
cubes. | 


FIGURE 1.25 
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1.8 The sequence {Tet,, }?-_, of tetrahedral numbers is {1, 4, 10, 20, 35, 
56, 84, 120, ...}, while the sequence {Pyr,,};~1 of pyramidal num- 
bers is {1, 5, 14, 30, 55, 91, 140, 204, ...}. Observe that the sequence 
{4 Pyr, }n=i is {4, 20, 56, 120, ...}. Formulate a theorem and prove 
it. Can you illustrate the theorem by stacking and re-stacking can- 
nonballs? 

1.9 The octahedral numbers are defined as the number of cannonballs 
in octahedral piles, as illustrated in Figure for the first five 
octahedral numbers 1, 6, 19, 44, and 85. 


+ & & & 


FIGURE 1.26 


Find a formula for the nth octahedral number. 


CHAPTER 2 


Congruence 


Mathematics is the queen of the sciences and 
number theory is the queen of mathematics. 


Johann Carl Friedrich Gauss 


No one who lacks an acquaintance with congruences 
can claim to know much about number theory. 


Underwood Dudley 


One of the most important relations in number theory—and one of 
the most beautiful—is congruence. We say that a is congruent to b mod- 
ulo m, and write a = b (mod m), if and only if m divides a— b where m, 
called the modulus, is a positive integer. Equivalently, a = b (mod m) 
if and only ifa = b + km for some integer k. The congruence symbol 
= was first employed by Gauss about 1800. For example, congruences 
modulo 2 express whether an integer is even (congruent to 0 modulo 2) 
or odd (congruent to 1 modulo 2). When a number k divides a number 
n one often writes k|n, i.e.,a = b (mod m) is equivalent to m|(a — b). 


2.1. Congruence results for triangular numbers 


In this section we consider some congruence theorems for triangu- 
lar numbers that will be useful in later chapters, for example, when we 
consider Pythagorean triples in Chapter/4]and perfect numbers in Chap- 
ter {7 Here we evaluate 7, modulo m for m equal to 2, 3, and 5. In each 
instance we partition triangular numbers to illustrate identities that es- 
tablish the congruences. 

In the n = 3 row of Table /1.1) the first two triangular numbers are 
odd, the next two are even, the next two odd, etc., which leads to the 
following theorem. 


Theorem 2.1. Forn = 1 we have 
T= 0 (mod 2), n=0Oor3 (mod 4), 
"(1 (mod 2), n=1or2 (mod 4). 
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Proof. In Figure |2.1| we illustrate that when k = 6, T,, = 10T, + 6T,_4, 
T4441 = 14T, + ra ae + 1, T4K42 = 14T, + 21 ped + 1, and T4143 = 
107; + 6T;,4,. These four equations can now be verified by simple alge- 
bra using T, = n(n + 1)/2. Hence we have Ty, = T4x43 = 0 (mod 2) 
and Typ44 = T4x42 = 1 (mod 2). 


T4441 = 147, + 214 +] T4k43 = 107; + 6744 
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T4x = 10T; + 671 


FIGURE 2.1 


Theorem 2.2. Forn => 1 we have 


_ 1 (mod 3), n= 1 (mod 3), 
~ (0 (mod 3), n#1 (mod 3). 


Th 


Proof 1. In Figure |2.2) we show that T3,_1 = 6T,_1 + 3T,, T3x = 6T, + 
3T,-1, and 73,4, = 9T, + 1 for k = 7, and the general case can be 
verified by algebra. Hence T3,_, = T3, = 0 (mod 3) and 73,4, = 1 
(mod 3). 
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T3441 = 6T;_] + 3T, 135 = 67; + 3Ty_] T3544 = 9T;, +1 


FIGURE 2.2 
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Observe that, as a bonus, we have 73,4; = 1 (mod 9). This congru- 
ence will be useful in Chapter 7|when we consider perfect numbers. 


Proof 2. Corollary {1.4] states that forall k = 2, Ty44 = 3T, — 3T,y-1 + 
Tz, SO that Ty4, = Ty-2 (mod 3). Hence any two triangular numbers 
whose ranks differ by a multiple of 3 are congruent modulo 3, and the 
theorem follows from T,; = 1 (mod 3), Tz = T; = 0 (mod 3). 


Theorem 2.3. Forn > 1 we have 


0 (mod 5), n=0Oor4 (mod 5), 
Tn = 41 (mod 5), n=10r3 (mod 5), 
3 (mod 5), n=2 (mod 5). 


Proof. In Figure b.3]we show that T;, = 157; +107;_1, T5441 = 20T;, + 
5Te-1 +1, Ts 42 _ 25T, oy Ts 43 = 20T, + 5T 41 + 1, and Ts44 = 
157, + 107T;4, for k = 5, and again the general case can be verified by 
algebra. Hence we have Ts, = Tsx44 = 0 (mod 5), Tsx44 = Tsx43 = 1 
(mod 5), and 75,42 = 3 (mod 5). 


T5% = 15T; + 107,_1 T5k41 = 20T, + STy_1 +1 
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FIGURE 2.3 


2.2. Congruence results for other figurate numbers 


In Example [1.1] we presented a test for triangular numbers: 8 times 
a triangular number plus 1 is a square. Now we present a visual illus- 
tration of that relationship and express it as a congruence result for odd 
squares. 


Lemma 2.4. Forn = 1 we have (2n + 1)° = 8T,,+1, so that (2n+1)? = 
1 (mod 8). 
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(2n+1)2 = 8T,, +1 
FIGURE 2.4 


Proof. See Figure for n = 10. Algebra verifies the result in general. 


In Example in the next chapter we evaluate squares modulo 3 
using triangles, similar to our evaluation of triangular numbers modulo 


3 in Theorem ?2.21 
Lemma {2.4 enables us to prove the following two theorems. 


Theorem 2.5. If the positive integer m is not divisible by 2 or 3, then 
m? = 1 (mod 24). 

Proof. The statement is clearly true form = 1 and 5. Every positive 
integer m = 6 has one of six forms: 6k, 6k + 1, 6k + 2, 6k + 3, 6k + 4, 
or 6k + 5 for some positive integer k. If m is not divisible by 2 or 3, 
then m = 6k + 1 orm = 6k + 5 since each of 6k, 6k + 2, 6k + 3, and 
6k + 4 is divisible by 2 or by 3. For the case m = 6k +1 settingn = 3k in 
Lemmal2.4lyields m? = [2(3k) + 47° = 873, +1. But from Theorem[2.2| 
we have 73, = 67; + 3T,_; and hence (6k + iy = 24(2T, +Ty-1) + 
1. For the case m = 6k + 5 letting n = 3k + 2 in Lemma yields 
m* = [2(3k +2) + i" = 873(x41)-1 + 1. But from Theorem we 
have T3x41)-1 = 67, +37, 41 and hence (6k + 5)° = 24(2T, + Trait 
1. Thus when m is not divisible by 2 or 3, m? = 1 (mod 24). 


In Example 1.8] we saw an example of multi-polygonal numbers, e.g., 
infinitely many numbers that are simultaneously square and triangular. 
Here is another proof of that result. 


Theorem 2.6. Infinitely many square triangular numbers exist. 

ee yields Tap, = Stn Stat) - 
4T,,(2n + he So if T, is a square, so is Tgr,. Since T; = 17, we have 
the following infinite sequence of square triangular numbers: T, = 17, 


Tg = 62, Tog = 204”, ware 


Proof. Setting k = 87, inT, = 
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However, the above proof only shows that there are infinitely many 
even square triangular numbers. There are also infinitely many odd 
square triangular numbers, such as Ty) = 35° and T1681 = 11897, but 
that requires a different proof. See Example £.8)in the next chapter. 

In Section we illustrated how pentagonal and triangular num- 
bers are related, e.g., P, = Ton-1 — Tn. That relationship leads to the 
next result. 


Theorem 2.7. Forn = 1 we have P, =n (mod 3). 


Proof. The result is clearly true for P, = 1. In Figure we see that 
P, = 3T,-1 + nforn = 2 (illustrated for n = 8). 
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Pr = Ty-1 - Th-1 = 3T,1+n 
FIGURE 2.5 


The fact that every polygonal number is a sum of triangular num- 
bers enables us to use the congruence results for triangular numbers in 
Section to evaluate similar congruences for n-gonal numbers with 
n => 4, as in the following example. 


Example 2.1. Heptagonal numbers modulo 3. From Section |1.1)the kth 
heptagonal number is Hep, = Ty + 47, _1. Thus Theorem {2.2} yields 


Hepsy_4 =T3x-1 + 47 3x-2 = 37 + 42Tx-1 + 4, 
Hep,, =T3x + 473x-1 =18T, + 27T,-1, 
Hepa = 13x41 + 473, = 337, + 12T,-1 +1, 


which proves 


ie a 0 (mod 3), n=0 (mod 3), 
"(1 (mod 3), n #0 (mod 3). 


Integer cubes are also related to triangular numbers, as we show in 
the following theorem. 


Theorem 2.8. Forn > 1 we have (n+ 1)° — n3 = 67, + 1, so that 
(n+ 1)* —n3 = 1 (mod 6). 
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Proof. In the first row of images in Figure 2.6] we see the difference of 
consecutive cubes as a “shell” of unit cubes. In the second row of im- 
ages we replace each unit cube by a small sphere, and project that array 
of spheres onto the plane as a hexagonal array at the right end of the 
second row consisting of a central black sphere surrounded by six tri- 
angular arrays of gray spheres. 


FIGURE 2.6 


In Exercise 3.13 in the next chapter you can study the existence of 
consecutive cubes whose difference is a square or, equivalently, if 67,,+1 
can equal a square. 


2.3. Fermat’s little theorem 


The statement in the following theorem first appeared in a letter 
Pierre de Fermat wrote to Bernhard Frénicle de Bessy (1605-1675) in 
1640. It is called his “little” theorem to distinguish it from his “last” or 
“great” theorem, the one about the lack of solutions to x” + y” = z” 
pictured on the postage stamp with Fermat at the end of Section 

There are a variety of proofs of this theorem. Ours is a combinatorial 


proof due to Solomon Golomb (1932-2016) [(Golomb, 1965; 
1971] 


Fermat’s Little Theorem 2.9. [fp is prime and n is a positive integer, 
thenn? = n (mod p). 


Proof. We shall show that p|(n? — n). Suppose we have a collection of 
beads inn different colors, and we wish to make necklaces (or bracelets) 
by stringing together exactly p beads. Our first step is to put p beads on 
a string. Since each bead can be chosen in n ways, there are n” possi- 
ble strings of p beads. Exactly n of these strings will be monochromatic, 
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and we discard these strings, leaving n” — n strings, each of which has 
beads of two or more colors. 

Next we join seamlessly the ends of each string to form a necklace. 
However, many of the resulting necklaces will be indistinguishable, as 
illustrated in Figure 2.7| when p=5andn=3. 


FIGURE 2.7 


So for each string S in the collection, there will be another p — 1 
strings that yield a necklace indistinguishable from S. Thus the num- 
ber of distinguishable necklaces is (n? — n)/p, which must be an inte- 
ger. 


Where in the proof did we use the hypothesis that p is prime? It was 
in the last step where we asserted that for each string there were p — 1 
additional strings that yield an indistinguishable necklace. If p is not 
prime, this does not hold. For example, let p be the composite number 
6, and consider the strings in Figure 


@ 00 0 
00-00 - 
0-0 -0-0- 


FIGURE 2.8 


Here a set of three beads is repeated twice, and so that there are only 
6/2 — 1 = 2 additional strings yielding the same necklace. 

Fermat’s little theorem is often stated in the form of the following 
corollary. 


Corollary 2.10. [fp is prime and nis a positive integer not divisible by p, 
thenn?~1 = 1 (mod p). 

Proof. From the little theorem we have p|(n” — n) or, equivalently, 
p|n(n?~1 — 1). Since the prime p does not divide n, it divides n?~1 — 1 
and the conclusion follows. 
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2.4. Wilson’s theorem 


In 1770 the English mathematician Edward Waring (1734-1798) 
published without proof a conjecture (the statement in the next theo- 
rem) of one of his students John Wilson (1741-1793). In the same year 
Lagrange gave the first proof. However, Gottfried Leibniz knew of the 
result before either Waring or Wilson was born. There are many proofs 
of the theorem; the following one is from 
1971]. 


Wilson’s Theorem 2.11. [fp is prime, then (p — 1)!= —1 (mod p). 


Proof. We shall show that p|[(p — 1)! +1]. The theorem is clearly true 
for p = 2 and p = 3 (2|(1! +1) and 3|(2! +1)), so assume that p = 5. 
Now choose p points on a circle, one directly above the center and the 
remaining p — 1 equally spaced around the circle. These points define 
the vertices of a regular p-gon (a polygon with p congruent sides and an- 
gles). We now construct stellated p-gons by connecting the vertices in 
every possible order (we allow edges to cross one another). The num- 
ber of distinct stellated p-gons is : (p —1)-(p — 2)-++-2-1 = (FE —1)!/2. 
Since the first vertex is the vertex directly above the center, the second 
can be chosen in p — 1 ways, the third in p — 2 ways, and so on to the 
last in one way; we divide by 2 since we can choose one or the other of 
the two edges that meet at the first vertex as the first edge (i.e., there 
are two ways in which to join the vertices to form the stellated p-gon). 
See Figure bdlfor the 4!/2 = 12 stellated pentagons. 


OWerwe 
& 2OHRa 


FIGURE 2.9 


Of the total number of stellated p-gons, (p — 1)/2 are regular stel- 
lated p-gons, in that the angles and side lengths are equal. These are 
constructed by connecting every kth vertex of the p-gon forl < k < 
(p — 1)/2. For p = 5 these are the two regular stellated pentagons in 
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FIGURE 2.10 


Figure and for p = 11 the five regular stellated 11-gons are illus- 
trated in Figure 

The [(p — 1) !/2] —[(p — 1)/2] non-regular stellated p-gons fall into 
subsets of size p, since for each such p-gon P there are p—1 others which 
can be obtained from P by rotating P through an angle of (360/p)’. 
When p = 5 there are two such subsets, the two rows of stellated pen- 


tagons in Figure 2.9(b)}. Thus [(p — 1)!— (p — 1)]/2p is an integer, so 


that p divides (p — 1)! +1. 
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2.6 


21 


2.5. Exercises 


Prove that in base 9, every integer consisting solely of 1’s is a tri- 
angular number, e.g., 11, =9+1=10=7,,111, = 814+9+1= 
91 = T;3, etc. (here the subscript on an integer indicates a base 
other than 10). [Hint: See Exercise 1.2. 

Examination of the columns of Table |1.1| leads to the observation 
that when k is odd, the kth n-gonal number is congruent to 0 mod- 
ulo k. Prove that this is indeed the case. 

Give a third proof of Theorem by proving the following 
three identities: 

(a) T3x-1 = 3 Tox — Te+1) + To; (b) T3~ = 3 Tox — Th); and (c) 
T3x41 = 3 (Tox41 — Tr+1) + Th. 

Prove that there exist infinitely many pairs (a, b) of integers such 
that for every positive integer t the number at + b is a triangular 
number if and only if t is a triangular number. [This is problem B6 
from the 1988 William Lowell Putnam Competition. ] 

Prove the following corollary to Fermat’s little theorem: [fp is an 
odd prime and nis a positive integer, thenn” =n (mod 2p). 

In Theorem 2.8|we proved that the difference of consecutive cubes 
is congruent to 1 modulo 6. Now prove that each cube n? is con- 
gruent to n modulo 6. 

Let n be a positive integer. Prove that in base 10, n and n° have the 
same units digit, e.g., isa = 16807 and 12° = 248832. [Hint: Two 
numbers have the same units digit if they are congruent modulo 
10.] 
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In Example 1.10|we mention Fermat primes. More generally, a Fer- 
mat number is an integer of the form f, = 22" + 1forn > 0. The 
first few Fermat numbers are fp = 3, f, = 5, fo = 17, fg = 257, 
fa = 65,537, and f; = 4,294,967, 297. Prove that for n = 2 (a) 
the units digit of f,, is 7, and (b) the tens digit of f,, is odd. 

Show that the nth Fermat number satisfies f, = 5 (mod 12) for 
n2>1. 

Observe that1 +3+6+10+15+ 21+ 28 = 84 = 3-28. Are 
there other instances where ¥,_, T; = 0 (mod Tp)? 

For what values of nis P, = 0 (mod 5)? 

Consider the determinant |a;;| of order 100 with a;; = i xj. Prove 
that if the absolute value of each of the 100! terms in the expan- 
sion of this determinant is divided by 101, then the remainder in 
each case is 1. [This is problem B1 from the 1957 William Lowell 
Putnam Competition. ] 

Show that p is the smallest prime that divides (p — 1)! +1. 


CHAPTER 3 


Diophantine Equations 


Equations are important to me, because politics is for 
the present, but an equation is something for eternity. 


Albert Einstein 


I do believe in simplicity. When the mathematician would 
solve a difficult problem he first frees the equation from 
all encumbrances, and reduces it to its simplest terms. 


Henry David Thoreau 


A Diophantine equation is an equation with integer coefficients in 
which the unknowns are also integers. They are named after Diophan- 
tus of Alexandria (circa 200-284 CE), who made an extensive study of 
equations whose coefficients and solutions were integers or rational 
numbers. In this chapter we discuss and illustrate linear Diophantine 
equations such as ax + by = c (equivalent to the congruence ax = c 
(mod b)) in Sections and 8.3} quadratic Diophantine equations of 
the form x? — dy? = 1in Sections 8.4 through and quadratic Dio- 
phantine equations of the form a? + b* = c? in the next chapter. 


Example 3.1. On the left in FigureB.1|we see a graph of the linear equa- 
tion 2x — 3y = 1 with some of its infinitely many integer solutions 
marked with ; a graph of the quadratic equation x? — 2y? = 1 with 
some of its infinitely many integer solutions in the center; and a graph 
of the quadratic equation x? + y? = 13 with all eight of its integer so- 
lutions on the right. 


FIGURE 3.1 
35 
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Diophantine equations and David Hilbert 


On August 8, 1900, David Hilbert (1862-1943) presented an address at 
the Second International Congress of Mathematicians in Paris. In this ad- 
dress, Hilbert presented 23 problems (in the printed version, but only 10 
in the oral presentation), the solutions to which he believed would pro- 
foundly impact the course of mathematical research in the twentieth cen- 
tury. In his 10th problem, Hilbert asked: Does there exist an algorithm to 
determine whether a given polynomial Diophantine equation has a solution 
in integers? This question was answered in the negative (no such algo- 
rithm exists) in the work of Martin Davis, Hilary Putnam, Julia Robinson, 
and Yuri Matiyasevich over a span of more than 20 years, ending in 1970. 


David Hilbert 


The above photograph was taken in 1912 for postcards of faculty mem- 
bers at the University of Gottingen. 


We begin in the next section with a way to illustrate some properties 
of squares using triangles. 


3.1. Triangles and squares 


In some instances it will be useful to use equilateral triangles to rep- 
resent integral squares. To do so, we first recall that every square is the 
sum of consecutive odd numbers, e.g., 


n?=14+34+5+-+-+(2n-1), 


as illustrated in Figure and in Figure for n = 8 (count the 
dots in each of the L-shaped regions). 


This enables us to represent n? using an equilateral triangle subdi- 
vided into n rows containing 1, 3, 5, ..., 2n — 1 identical smaller equi- 
lateral triangles (A and V), as seen in Figure S2{0) As an application, 
we color the small triangles as shown in Figure to illustrate the 
identity T,, + T,-, = n? (for n = 8) with T,, small gray triangles and 
Tn—1 Small white triangles. 
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FIGURE 3.2 


Example 3.2. Squares modulo 3. The n = 4 row in Table [1.1] yields the 
hypothesis that every square is either a multiple of 3 or 1 more than a 
multiple of 3. In Theorem R.2|we used a triangular array of dots to eval- 
uate triangular numbers modulo 3. With triangular diagrams similar to 
those in Figure B.2|we can evaluate squares modulo 3: 


0 (mod 3), n=O (mod 3), 
(mod 3), n#0 (mod 3). 


2 


See Figure and note that a set of small triangles in the shape of a 
trapezoid represents the difference of two squares. The figure is drawn 
for k = 3, representing n? forn equal to 8, 9, and 10. 


(3k)2 = 3[(2k)2 — k2] 


III 
LISSA 
KITA N\A 


(3k-1)2 = 143[(2k-1)2 — (k-1)2] (3k+1)2 = 143[(2k+1)2 — (k+1)2] 
FIGURE 3.3 
Example 3.3. Fermat’s little theorem for p = 3. Figure also illus- 


trates the little theorem for p = 3 in the form expressed in Corollary 
i.e., ifn is an integer not divisible by 3, then n? = 1 (mod 3). 


38 3. DIOPHANTINE EQUATIONS 


Example 3.4. An identity for seven squares. For positive integers a, b, 
and c, we have 


a2+b2+c24+(atbtc) =(atb) +(atc)+(b+c). 


In Figure B.4|we use the inclusion-exclusion principle to count the small 
A and V triangles, obtaining 


(a+b+c) =(atb)’ +(atc)+(b+c)* — a? — b? - c?, 


from which the desired result follows. Simple algebra shows that the 
identity actually holds for all real numbers. 


NAAN 
KIN/\TV. 
WAVAV, 


\ 


a b c 


FIGURE 3.4 
In the figure (a, b,c) = (3,4,5), which yields 37 + 47 + 57 + 127 = 
7? + 8% + 92. The change of variables 
(a,b,c) = (—*e+y42,%-—yr+2,%*«+y7—2) 
yields 


(xtytz)*+(xty—z)*+(x-ytz)*+(—xt+yt4z)* = 4(x*+y*%4+2z7). 


Other identities for squares can be obtained in a similar fashion. 


3.2. Linear Diophantine equations 


A linear Diophantine equation is an equation of the form ax + by = 
c, where a, b, and c are given integers and x and y are integer-valued 
unknowns. The task is to solve the equation, finding integer values of x 
and y so thatax + by = c. 
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A few observations: 


1. Letd = gcd(a,b), the greatest common divisor of a and b. 
Hence d is always a factor of ax + by. So if d is not a factor 
of c, then ax + by = c is not solvable in integers. 


2. Assume d is a factor of c. Then we can divide through by d 
ag us and < are integers and 
d’d d 

gcd(S, 7) = 1. Furthermore, any integer values of x and y that 


- a b Cc 
to obtain qx + ao = ai where 


satisfy a - -y = - also satisfy ax + by = c, and vice versa. 

3. Hence solving ax + by = c is reduced to the case of solving 
the equation when gcd(a, b) = 1. But as a consequence of the 
Euclidean algorithm for finding gcd(a, b), there exist integers 
r ands suchthatar+bs = 1,henceax+by = c hasasolution 
Xp =rcand yp = Sc. 


The Euclidean algorithm 


The Euclidean algorithm is an arithmetic procedure for finding the 
greatest common divisor of two positive integers a and b, so named since 
it appears in Book VII of Euclid’s Elements. Here is a geometric version 
iKiizek et al., 2001; Walser, 2001]. Draw an a-by-b (a > b) rectan- 
gle, as shown below. Cut off as many squares with side b as possible, 
until a rectangle remains. Again cut off squares from the new rectangle, 
and continue until no rectangle remains. The side length of the smallest 
square is then gcd(a, b). The figure below illustrates the procedure with 
(a, b) = (64,27) and gcd(64, 27) =1. 


64 


A method for finding all solutions of ax + by = c (when gcd(a, b) = 
1) is based on the following lemma. 


Lemma 3.1. /f a,b, x, and y are integers, thenax + by = a(x —b) + 


b(y + a). 
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b 


FIGURE 3.5. ax + by = a(x —b)+b(y +a) 


Proof. The result follows from simple algebra. In Figure B.5|we illustrate 
the result for positive integers by counting the number of small squares 
in two different ways as indicated by the shading. 


Similar algebra establishes ax + by = a(x —kb) + b(y + ka) for 
any integer k. Alternatively, you can also derive this from the formula in 
the lemma by replacing x by x + b, x + 2b, etc., and similarly replacing 
ybyy+a,y + 2a, etc. This leads to the following theorem. 


Theorem 3.2. Suppose a and b are non-zero integers with gcd(a, b) = 1, 
and that x9, Yo isa solution toax+by = c. Thenallsolutions toax+by = 
c are given by x = Xy9 — kb, y = yo + ka for some integer k. 


The proof that all solutions have the stated form can be found in 


most number theory textbooks, e.g., [Marshall et al., 2007]. 


Linear Diophantine equations are often encountered in recreational 


mathematics puzzles, such as the one in the following example adapted 
from [Gardner 2001]. 


Example 3.5. Three sailors and a monkey were shipwrecked on a 
desert island and spent a day gathering coconuts, as illustrated in a Ja- 
vanese drawing in Figure 3.6 

They piled the coconuts together, and agreed to divide them up after 
a nights rest. During the night one sailor wakes up, divides the coconuts 
into three equal piles with one left over, which he gives to the monkey. 
He hides his share, puts the remaining coconuts in a pile, and goes back 
to sleep. The second sailor does the same thing, as does the third. In 
the morning the sailors divide the remaining pile of coconuts, with one 
left over for the monkey. How many coconuts were there in the original 
pile? 

Let x denote the number of coconuts in the original pile, a, b, and 
c the number taken by each sailor during the night, and y the number 
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FIGURE 3.6 


each received in the morning. Then x = 3a+1, 2a = 3b+1,2b = 3c+1, 
and 2c = 3y+1. Eliminating a, b, andc yields the Diophantine equation 
8x — 81y = 65. Since (x, y) = (—10, —1) is a solution to 8x — 81y = 1, 
(Xo, Yo) = (—650,—65) is a solution to 8x — 81y = 65. Hence the 
general solution is (x,y) = (—650+ 81k, —65+ 8k). Fork = 9 we have 
the positive solutions (x, y) = (79, 7), (160, 15), (241, 23), etc. 


This problem originally appeared in the October 9, 1926 issue of 
The Saturday Evening Post, with five sailors rather than three 
2001]. You are invited to solve the original version in Exercise 3.1. 


3.3. Linear congruences and the Chinese remainder theorem 


Knowing how to solve a linear Diophantine equation enables one to 
solve a linear congruence, a congruence of the form ax = c (mod b) for 
x, since the congruence is equivalent to the linear Diophantine equation 
ax + by = c. So we now consider solving a set of simultaneous linear 
congruences, e.g., finding values of x that satisfy x = c; (mod m;) for 
1 <i<k. The key to solving such systems is the following theorem, 
whose proof can be found in a number theory textbook. 


The Chinese Remainder Theorem 3.3. The system of k linear congru- 
ences x = c; (mod m;) for 1 < i < k, where the moduli are relatively 
prime in pairs, has a unique solution modulo mym «+: Mx. 


Note that the theorem tells us when the system has solutions, but 
does not provide a procedure for finding them. The name of the the- 
orem honors the contributions of early Chinese mathematicians to the 
theory of congruences. There are a variety of methods for finding solu- 
tions to systems of linear congruences. The following two examples il- 
lustrate a tabular procedure applicable to systems with relatively small 
moduli. 
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Example 3.6. Consider a system of two congruences: x = 2 (mod 3) 
and x = 3 (mod 5). The Chinese remainder theorem guarantees the 
existence of a solution mod 15. To find it, we set up a 3-by-5 table with 
row labels mod 3 and column labels mod 5. Highlighting the row and 
column corresponding to the given congruences identifies the cell that 


will contain the solution; here it is cell (2,3). See Table B.11 


TABLE 3.1 


We now successively enter the numbers 0, 1, 2, ..., starting in the 
(0,0) cell and proceeding down the diagonal until we reach an edge of 
the table. When we reach the bottom edge we jump up to the top edge, 
and when we reach the right edge we revert to the left edge of the table, 
as if the table were the map ofa torus, where the top and bottom edges 
of the table are adjacent, as are the right and left edges. See Figure 


FIGURE 3.7 


When we reach cell (2,3) in the table, we stop entering numbers as 
we have found the solution: x = 8 (mod 15). 


Example 3.7. Sun Tzu’s puzzle. Perhaps the earliest use of the Chinese 
remainder theorem appears in Problem 26 of Chapter in Sun Tzu Suan 
Ching (Sun Tzu’s Mathematical Manual), circa 250 CE (see Figure B.8): 


“Suppose we have an unknown number of objects. 
When counted in threes, 2 are left over, when counted 
in fives, 3 are left over, and when counted in sevens, 
2 are left over. How many objects are there?” 
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FIGURE 3.8. Facsimile of Qing dynasty edition of Sun Tzu’s 
Mathematical Manual 


Solving Sun Tzu’s puzzle requires solving the system 


x = 2 (mod 3), 
x =3 (mod 5), 
x = 2 (mod 7). 


We can replace the first two equations with x = 8 (mod 15) from Ex- 
ample and solve the system 
x = 8 (mod 15), 
x =2 (mod 7), 
using our tabular method, as illustrated in Table 
When we reach cell (2,8) we stop entering numbers, as we have 


found the solution x = 23 (mod 105). So the number of objects in Sun 
Tzu’s puzzle is one of the numbers 23, 128, 233, 338,.... 


TABLE 3.2 


mod 7 
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3.4. The Pell equation x* — 2y* = 1 


Some of the most often studied quadratic Diophantine equations are 
the so-called Pell equations x* — dy* = c, where c is an integer and d is 
a positive non-square integer. Such equations may have infinitely many 
solutions. When d is a negative integer or a positive integral square, the 
Pell equation has at most finitely many solutions. 

These equations were mistakenly attributed to a 17th century En- 
glishman John Pell by Leonard Euler in 1730. This was unfortunate, for 
as the mathematical historian E. T. Bell [Bell, 1990] writes, “Pell math- 
ematically was a nonentity, and humanly an egregious fraud. It is long 
past time that his name be dropped from the textbooks.” 

Some authors refer to these equations as Fermat equations, since 
Pierre de Fermat initiated their modern history with a letter in 1657 to 
several mathematicians challenging them to solve x?—61y” = 1 ininte- 
gers. However, these equations were studied in India by Brahmagupta 
(circa 598-665 CE) and others long before Fermat. 

The standard way of solving Pell equations is via the theory of con- 
tinued fractions. Our approach to solving x? — 2y? = 1 is more combi- 
natorial and is based on counting the small squares in Figure 


x+y x+2y x+y 


(b) 


FIGURE 3.9 


Counting the number of small squares in Figure using the 
inclusion-exclusion principle yields (x + 2y)* =2(x+ y) —x?42y?, 
from which it follows that x2 — 2y? = 2(x+y)* —(x+2y)?. Simi- 
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larly counting the small squares in Figure yields (3x + 4y)? = 
Z2( 2x + 3y)° = (-F 2y)* + 2(x + y), so that 2(x + yy —(«+ 2y)° a 
(3x + 4y)° —2(2x + 3y)”. Hence x?—2y? = (3x + 4y)* —2(2x + 3y)’ 
and thus 


(3.1) x?-—2y? = 1 ifand only if (3x + 4y)* — 2 (2x + 3y)’ =1. 


We nowuse (/3.1)) to generate a sequence (xX, yp) Satisfying x2—2y2 = 1 
using the recursions X¥y41 = 3Xy + 4Vn, Vn41 = 2Xn + 3y¥y and start- 
ing with an initial solution (x9, yo) = (1,0). The results are shown in 
Table for n from 0 to 5. 


TABLE 3.3 


n|O 12 3 4 #5 
X,|1 3 17 99 577 3363 
yn |0 2 12 70 408 2378 


We now find explicit expressions for the solutions x, and y,. The 
recursions Xn41 = 3Xy) + 4¥n and Yn44 = 2X, + 3y, imply that 
Xny1 = 3X_ + 4(2Xp_1 + 3Yn-1) 
= 3%, + 8Xnai +3 (4%_— 3X p4) 
= 6%, = Xa 
Yn+1 = 2(3Xn-1 + 4Vn-1) + 39 
= 3 (Yn — 3Yn-1) + BYn-1 + 3Yn 
= 6Yn — Yn-1) 
1.€., X44 = OX,—Xp_1 aNd Vp 44 = OVn —Vn_1- Evaluating the generating 
functions G,(t) = ¥,-9 Xnt™ and Gy(t) = 0 Ynt™ for the sequences 
{x,} and {y,,} yields 
G,(t) = (1 - 3t)/(1—-6t +t?) and G,(t) = 2t/(1—6t + t?), 
from which we obtain 
1 n n 
mn = 51 +2V2) +(3-2v2) ] 
and 


1 n n 
Yn = 5 gIG + 2V2) — (3 -2v2) ] 


via partial fraction decompositions and geometric series. 


46 3. DIOPHANTINE EQUATIONS 


Example 3.8. Square triangular numbers. In Example 1.8 we saw that 
infinitely many square triangular numbers exist. We obtained the same 
result in Theorem?.6] but the proof of that theorem generated only even 
square triangular numbers. In this example we generate the infinite se- 
quence of square triangular numbers T,,, = k* using the Pell equation 
x? — 2y* = 1. Completing the square for both m and k shows that 
Tm = k? is equivalent to (2m + 1) a 2(2k)° = 1, hence each solu- 
tion (Xn, Yn) in Table B.3] yields a square triangular number. For exam- 
ple (x2,¥2) = (17,12) yields (m,k) = (8,6) so that Tg = 36 = 67, 
(x3, ¥3) = (99,70) yields (m,k) = (49,35) so that Tyg = 1225 = 
35°, (x4, V4) = (577,408) yields (m,k) = (288,204) so that Togg 
41616 = 2047, and so on. 


Solutions to the Pell equation x* — 2y” = c for another integer c are 
found similarly. If we have x*—2y? = —1 we have the same recursion as 
above but with an initial solution (x9, yp) = (1, 1), yielding the results 


in Table B.4) for n from 0 to 5. 


TABLE 3.4 


n{O 12 3 4 5 
Xn |1 7 41 239 1393 8119 
yn |1 5 29 169 985 5741 


ple |1.8, infinitely many positive integers are both oblong and triangu- 
lar: Tz = 6 = 2-3 and Tz) = 210 = 14- 15, for example. As in the 
preceding example, we can generate the sequence of oblong triangular 
numbers by considering T,, = k(k + 1) and completing the square to 
obtain (2m + 1)° —2(2k + 1) = —1. Thus the solutions (%y, yp) in 
Table 3.4; to the Pell equation x* — 2y? = —1 yield oblong triangular 
numbers. For example (x3, y3) = (239, 169) yields (m, k) = (119, 84), 
and T;,5 = 7140 = 84-85. 


isin: 3.9. Oblong triangular numbers. As we observed in Exam- 
1.8 


3.5. The Pell equation x* — 3y* = 1 


To solve this Pell equation, we first count the small triangles in Fig- 
ure using the inclusion-exclusion principle, recalling that each 
large triangle contains a square number of small triangles, as illustrated 
in Figure 

This yields (2x + 3y)* = 3(x + 2y)* + x2 — 3y2, and hence 


(3.2)  x?—3y? =1ifand only if (2x + 3y)* — 3(x + 2y)* = 1. 
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x+y y x+y 


FIGURE 3.10 


We now use ((3.2)) to generate a sequence (xy, Y,) Satisfying x2—3y? = 
using the recursions Xn44 = 2X» + 3Vy, Vn41 = Xn + 2¥n and starting 
with an initial solution (x9, yo) = (1, 0), yielding the results in Table 
for n from 0 to 6. 


TABLE 3.5 
012 3 4 5 6 
x,'1 2 7 26 97 362 1351 
01 4 15 56 209 780 


We now find explicit expressions for the solutions x, and y,. The 
recursions Xn44 = 2Xy + 3p and Yn44 = Xp + 2y, imply that x,41, = 
AXn — Xn—1 ANd Vn44 = 4Vn — Yn_1- Evaluating the generating functions 
G,(t) = 28 x,t" and Gy(t) = aa y,t” for the sequences {x,,} and 
{Yn} yields G,,(t) = (1 — 2t)/(1 — 4t + t?) and G,(t) =t/(1 — 4t + t*), 
from which we obtain 

1 n n 
Xm = 5[(2+V3) + (2- V3) |] 
and i 
n n 
= — [(2+ V3) -(2-v3 
eo — ait | ) | 


via partial fraction decompositions and geometric series. 


Example 3.10. Almost equilateral Heronian triangles. A Heronian tri- 
angle is one whose sides (a, b,c) and area K are integers (some sources 
replace “integers” with “rational numbers”). An almost equilateral 


48 3. DIOPHANTINE EQUATIONS 


Heronian triangle is a Heronian triangle whose sides are consecutive 
integers, such as (3, 4,5) with area 6 and (13, 14, 15) with area 84. Are 
there others? In answering, we recall that the area K of an arbitrary 
triangle (a, b,c) is given by Heron’s formula: 


K =.s(s—a)(s—b)(s—c), 


where s denotes the semiperimeter, s = (a+b+c)/2. 

Letting (b — 1,b,b + 1) denote the sides of the triangle yields 
s = 3b/2 and K* = (3b/2) (b/2) [(b — 2)/2] [(b + 2)/2] or 16K? = 
3b7(b? — 4). Thus b must be even, so we set b = 2x, hence K? = 
3x?(x?-—1). Then x*—1 must be 3 times a square, so that x?—1 = 3y?, 
the Pell equation solved in this section. So pairs (x, y,,) in Table 
yield triangles with sides (2x, — 1, 2x,,2x, + 1) and area K = 3xyyy, 
(thus the altitude h to the even side 2x,, is the integer h = 3y,,). Hence 
there are infinitely many almost equilateral Heronian triangles. For ex- 
ample, n = 3 yields the triangle (51, 52,53) with area K = 1170. The 
inradius (the radius r of the inscribed circle) of an almost equilateral 
Heronian triangle is also an integer; in Exercise 3.16 you can prove that 
Yr = yp. In Figure B.11] we see an illustration of the almost equilateral 
Heronian triangle (2x, — 1, 2x, 2x, + 1) and its altitude and inscribed 
circle for the case (x3, y2) = (7,4). 


FIGURE 3.11 


Solutions to the Pell equation x* — 3y” = c for another integer c are 
found similarly. Ifwe have x*—3y* = —2 we have the same recursion as 
above but with an initial solution (x9, yp) = (1, 1), yielding the results 
in Table B.6 for n from 0 to 6. 
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TABLE 3.6 
n | 01 2 3 4 5 6 
X,i'1 5 19 71 265 989 3691 
y,|1 3 11 41 153 571 2131 


Example 3.11. Pentagonal triangular numbers. In Chapter /1} we en- 
countered numbers that were simultaneously pentagonal and triangu- 
lar: P, = 1 = T, (which is rather uninteresting) and P,, = 210 = 
T29, for example. Are there more? How many? As we did with square 
triangular numbers, we look to solve P, = Ty, i.e, m(3m—1)/2 = 
k(k + 1)/2 for m and k. Completing the square on m and k yields 
(6m — 1)? — 3(2k +1)” = —2, the Pell equation x? — 3y? = —2 with 
x = 6m—1andy = 2k + 1. Hence the pairs (x, y,) in Table 
with x,, = —1 (mod 6) yield solutions, e.g., (x3, y3) = (71,41) yields 
(m,k) = (12,20) with Pjy = 210 = Tro, (%s, ys) = (989,571) yields 
(m, k) = (165, 285) with Pig; = 40755 = Togs, etc. 


3.6. The Pell equations x? — dy” = 1 


We now consider x* — dy” = 1, where d is positive and non-square. 
Note that (x,y) = (1,0) is a solution for all d, and the only solution 
when d is square. A proof that x? — dy* = 1 (for d positive and non- 
square) has infinitely many solutions in positive integers whenever it 
has at least one such solution is based on the following simple lemma. 


Lemma 3.4. For u,v > 0, (u+ v) = (u — v) + 4uv. 


Proof. Without loss of generality assume u = v and see Figure for 
a visual proof. Of course the lemma can also be proved by elementary 
algebra. 


Theorem 3.5. Ifx? — dy” = 1, where d is positive and non-square, has 
a solution in positive integers, then it has infinitely many. 


Proof, Set u = x? and v = dy? in LemmaB4|to yield (x? + dy?) = 
(27 ~ dy?) + 4dx*y? or, equivalently, (x? + dy?) — d(2xy) 
(x? = dy?) So if x* — dy? = 1, then (x? + dy?) —d(2xy)” = 1, 
and hence if (x, y) is a solution to x* — dy? = 1 in positive integers, so 
is (x? + dy”, 2xy). 


The recursion in the proof of Theorem {3.5} generates an infinite se- 
quence of solutions, but the sequence may not contain all solutions. For 
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FIGURE 3.12 


example, d = 2 and (x,y) = (3,2) generates (17,12), (577,408), and 
so on, omitting half the solutions to x? — 2y* = 1 in Table 

It is not obvious that x? — dy” = 1 always has a solution in positive 
integers (when d is positive and non-square). For example, the smallest 
solution in positive integers to x? — 29y? = 1 is (x,y) = (9801, 1820), 
and the smallest solution in positive integers to x*—61y? = Lis (x,y) = 
(1766319049, 226153980). Joseph-Louis Lagrange proved in 1766 
that x? — dy? = 1 (where d is positive and non-square) has infinitely 
many solutions in positive integers ee, 

The combinatorial approach in the previous two sections (counting 
small squares and triangles) to solving Pell equations x* — dy* = 1 
becomes difficult for d > 5. But examination of the results for d = 2 
and d = 3 reveals a pattern: 


x? — 2y? = (3x + 4y)? — 2(2x + 3y)’, 
x? — 3y* = (2x + 3y)* — 3(x + 2y). 


Replacing y by 2y in x? — 2y? = (3x + 4y)° — 2(2x + 3y)’ yields 
x? — By? = (3x + By)? — B(x + 3y)?. 


The terms to the right of the equal sign each have the form (ax +dby) — 
d(bx + ay)”, where a? — db? = 1, thatis, (a, b) is the smallest solution 
to x*—dy? = 1 in positive integers from Tables.3|andB.5 and (a, b) = 
(3,1) ford = 8. This is readily verified for d = 5, 6, and 7. For a 
5y? = 1 wehave (a, b) = (9,4), forx?—6y? = 1 we have (a,b) = (5, 2), 
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and for x* — 7y* = 1 we have (a, b) = (8,3) so that 
x? — 5y* = (9x + 20y)? — 5(4x + 9y)?, 
x? — 6y* = (5x + 12y)? — 6(2x + 5y)?, 
x? — 7y? = (8x + 21y)? — 7(3x + 8y)?, 
and so on. 

We now generate a sequence (Xn, yp) satisfying x2 — dy? = 1 using 
the recursions X744 = AXn + dbDyYn, Vn41 = DXn + Ayy and starting with 
an initial solution (Xo, yo) = (1, 0). Note that (a, b) = (x4, y,), so that 

x? — dy? = (xx + dy) — d(x + xy)”. 
The recursions Xn44 = 1X, + dyi¥n and Yn41 = W1X, + X1Vy_ With 
x? — dy? = 1 imply that 
Xnt1 = XX + dy 1Xn-1 + X1Yn-1) 
= X4Xy + dYPXn—-1 + X1 (Xn — X1Xn-1) 
= 2X4Xy_ + (dy? eae x7) %n-1 
= 2X4Xp — Xn-1 


and 


Yn+1 = Vi%1%Xn-1 + AV1Yn-1) + X1Yn 

= X1(Yn — X1Yn—1) + AYP Yn—1 + X1Yn 

= 2x1 + (dy? — xf) Yn-1 

= 2X1Yn — Yn-v 
Le. Xn41 = 2X1X, — Xp-1 ANd Yny1 = 2X1V_ — Yn-1. Evaluating the 
generating functions G,(t) = aan X,t” and Gy(t) = ae y,t” for the 
sequences {x,} and {y,} yields G,(t) = (1 — x,t)/(1 — 2x,t + t”) and 
G,(t) = y,t/(1 — 2x,t + t*), from which we obtain 


1 n n 
m= 5lGa+uvd) +a -yVd) | 
and 


1 n n 
Yn = sal + yvd) —(x1-yVa) | 


via partial fraction decompositions and geometric series. 


Example 3.12. Square pentagonal numbers. At first glance it appears 
that the trivial case P, = 1 = 1? may be the only square pentago- 
nal number. But P,, = k? is equivalent to m(3m—1) = 2k? which, 
upon completing the square, yields (6m — 1)° — 6(2k)* = 1. Using 
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TABLE 3.7 


n|O 12 3 4 5 
Xn|1 5 49 485 4801 47525 
yn |0 2 20 198 1960 19402 


x2 — oy? = (5x + 12y)* — 6(2x + 5y)’ to solve x? — 6y? = 1 with the 
initial solution (Xo, yo) = (1, 0) yields the results in Table }3.7| for n from 
Oto 5. 


Since x = 6m — 1 we require x, = —1 (mod 6), so that n must be 


odd. Thus (x3,y3) = (485, 198) yields Pj, = 9804 = 997, (xs, ys) = 


47525, 19402) yields P. = 94109401 = 97017, and so on. 
( y 7921 


3.7. Exercises 


3.1 Solve the “coconuts” problem in Example with five sailors 


rather than three. 


3.2 Do there exist 1,000,000 consecutive integers each of which con- 


tains a repeated prime factor? (This is Problem B4 from the 1955 
William Lowell Putnam Competition.) 


3.3 Prove that if (xp, ¥,) is a solution to x? — 2y? = 1, then x;y, isa 


multiple of 6. 


3.4 Prove that if (x,, ¥p,) isa solution to x? — 2y* = 1 or to x* — 2y? = 


—1, then (Xn¥n) is a square triangular number. 


3.5 Ramanujan’s house number problem. The following problem, 


which appeared in the December 1914 issue of the British 
magazine The Strand, is often mentioned in biographies of the In- 
dian mathematician Srinivasa Ramanujan (1887-1920). The bi- 
ographies tell us that Ramanujan solved the problem while stirring 
vegetables in a skillet as a colleague read the problem to him. 

In a certain street there are more than 50 but less than 

500 houses in a row, numbered from 1, 2, 3, etc., con- 

secutively. There is a house in the street, the sum of all 

the house numbers on the left side of which is equal to 

the sum of the house numbers on its right side. Find the 

number of this house. 
Solve the problem. 


3.6 Powerful numbers. An integer is called powerful if the square of 


each of its prime divisors also divides it, i.e... no prime divides a 
powerful number to the first power only. Two pairs of consecutive 
powerful numbers are (8,9) = (23,37) and (288,289) = (2° - 
a) 17”); Prove that infinitely many pairs of consecutive powerful 
numbers exist. 


ait 


3.8 


a9 


3.10 


3.11 
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Prove that there exist infinitely many integers n such that n,n + 1, 
and n + 2 are each the sum of two squares (of integers). Example: 
8 = 274+ 27,9 = 37+ 07,10 = 32 + 12. [This is problem A2 from 
the 2000 William Lowell Putnam Competition. | 

Observe that 3? + 42 + 5% = 77 +1. Are there other instances 
where the sum of the squares of three consecutive integers is one 
more than a square? 

In Table we see that T, = 37, and Tg = 375, two instances 
where one triangular number equals three times another. Are 
there other instances of this phenomenon? 

In Example f1.2| we observed that 225 is both square and octagonal. 
Show that there are infinitely many square octagonal numbers. 
At the end of Section we introduce the star numbers S, = 
12T,_-; +1 = 6n(n—-1) +1, and illustrated S$; = 121, ie. Ss; = 
11*. Are there other square star numbers? 


3.12 Almost square triangular numbers. In the game of pool one is given 


3.13 


16 balls arranged in a square tray, as shown on the left in Fig- 
ure When the white cue ball is removed from the tray, the re- 
maining 15 balls can be placed ina triangle, as shown on the right, 
since 42 -1=15=Ts. 


FIGURE 3.13 


More generally, an almost square triangular number is a triangular 
number that differs from a square by 1, such as T, = 3 = 27 — 1, 
T, = 10 = 37+ 1,7; = 15 = 42 — 1 (as in the game of pool), and 
Tog = 325 = 18° + 1. Show that there are infinitely many almost 
square triangular numbers. 

In 13—03 = 12 and83—73 = 13° wesee that the difference of con- 
secutive cubes is sometimes a square. Are there other instances of 
this phenomenon? 


3.14 Show that there are infinitely many positive integers k such that 


both k + 1 and 3k + 1 are perfect squares. 
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3.15 In Figure we saw that T,_; + T, = k?, that is, every pair 
of consecutive triangular numbers has a square sum. Observe that 
Ts + Tg + T7 = 15+ 21+ 28 = 64 = 87. Are there other sets of 
three consecutive triangular numbers that have a square sum? 

3.16 Theinradiusr ofatriangle is the radius ofits inscribed circle. Prove 
that the inradius of an almost equilateral Heronian triangle (see 
Example is an integer. [Hint: First show that for an arbitrary 
triangle, K = rs where K is the area and s is the semiperimeter.| 

3.17 Anurn contains k balls, of which r are red and k —r are white. Two 
balls are drawn at random without replacement, and the proba- 
bility that both are red is 1/2. What are the possible values of k 
andr? 


CHAPTER 4 


Pythagorean Triples 


The Hypotenuse has a square on, 

which is equal Pythagoras instructed, 

to the sum of the squares on the other two sides 
If a triangle is cleverly constructed. 


Richard Digance 


All is number. 


Attributed to Pythagoras 


The Pythagorean theorem, relating the squares of the lengths of two 
sides of a right triangle to the square of the length of the third side, may 
well be the best-known theorem in geometry—maybe in all of math- 
ematics. Of interest in number theory are Pythagorean triples (PTs), 
triples (a,b,c) of positive integers with a2 + b? = c%, so that a, b, 
and c can be side lengths in a right triangle. Examples include (3, 4, 5), 
(6, 8, 10), and (5, 12,13). When a, b, and c are relatively prime—as in 
(3, 4,5) and (5, 12, 13)—we call the triple a primitive Pythagorean triple 
(PPT). Note that if (x, y,z) is a PT where the greatest common divisor 
of x, y, and z is d, then ee a “) is a PPT. 

It is easy to show (see Exercise 4.1) that if (a, b,c) is a PPT, then a 
and b have opposite parity, and c is odd. For convenience in what fol- 
lows, we will usually write (a, b,c) with a odd and b even. A Pythago- 
rean triangle (PA) is a triangle whose side lengths are the numbers ina 
PT. In Figure we see the PA (3, 4,5) as represented in the Zhou 
Bi Suan Jing (el##5-), a Chinese text dating from the Zhou dynasty 
(1046-256 BCE); in Figure /4.1(b)|on a Greek postage stamp from 1955; 
and in ee on a Philippine postage stamp from 2001. 

Figure is an ancient but simple visual proof of the Pythagorean 
theorem based on the Zhou Bi Suan Jing. 
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FIGURE 4.1 


FIGURE 4.2 


4.1. Euclid’s formula 


In Book X of Euclid’s Elements, Lemma 1 preceding Proposition 29 
states: To find two square numbers such that their sum is also square. The 
procedure Euclid employs to accomplish this is equivalent to, in modern 
notation, the following theorem. 


Theorem 4.1. Euclid’s formula for PPTs. Letm and n be relatively prime 
positive integers with opposite parity with m > n. Then (m? — n?,2mn, 
m? +n?) is a PPT, and every PPT has this form. 


Figure 4. 3jis an illustration of how the pair (m,n) generates the PPT 
(m2 — n?,2mn,m? + n?). 

We call (m,n) the generator of (m2 — n?,2mn,m? + n?). A proof 
that the formula generates PPTs and that every PPT is generated by the 
formula can be found in most elementary number theory textbooks, e.g., 
[Marshall et al., 2007]. 

Euclid’s formula can also be illustrated using the double angle for- 
mula for the tangent function from elementary trigonometry. If we let 
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FIGURE 4.3 


6 denote the angle marked in the left-most triangle in Figure 4.3, then 


ania ips Une 
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4.2. Pythagorean triples and means of odd squares 


There are many ways to compute the average of two positive num- 
bers x and y. Such averages are called means, from the French moyen, 
“medium” or “middle.” Two of the most common are the arithmetic 
mean = and the geometric mean ./xy. These means are related by the 
arithmetic mean-geometric mean inequality, which states: Ifx and y are 
positive, then see > ./xy. A simple proof of the inequality follows from 
the observation that the two means are the hypotenuse and one leg of 
a right triangle, and that the hypotenuse is always the longest side. See 
Figure 

Now let x = s* and y = t*, where s and t are relatively prime 
odd positive integers with s > t, as shown in Figure Then 


FIGURE 4.4 
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s?-t? s%+t?\., : : : : er 
is ) is a PPT, since this form is equivalent to Euclid’s for- 
mula in Theorem upon setting s = m+nandt = m—n. This 
is the form for PPTs that appears in Euclid’s Elements. Consequently 
PPTs arise from computing the geometric and arithmetic means of two 
relatively prime odd squares (the even leg is half the difference of the 
squares). For example, the squares 25 and 49 yield the PPT (35,12,37). 


(se, 


4.3. The carpets theorem 


The carpets theorem is a little known but quite useful theorem in 
number theory and combinatorics. We shall use it in several sections in 
this chapter. 


The Carpets Theorem 4.2. Place two carpets ina room. The area of the 
floor equals the combined area of the carpets if and only if the area of the 
overlap equals the area of the uncovered floor. 


See Figure 4.5] for an illustration of two elliptical carpets in a rectan- 
gular room. 


FIGURE 4.5 


Proof. Expressing the area of the floor using inclusion-exclusion yields 
(where A denotes area) 


A(floor) = A(carpet X) + A(carpet Y) 
— A(overlap) + A(uncovered floor), 
and thus 
A(floor) = A(carpet X) + A(carpet Y) 
© A(overlap) = A(uncovered floor). 
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The shapes of the room and the carpets are arbitrary, and the result 
holds for more than two carpets as long as at most two carpets overlap 
simultaneously. In the next chapter we use the carpets theorem to give 
proofs that V2 andv3 are irrational. 

4.4. Pythagorean triples and the factors of even squares 


Several of the results we are about to derive for PTs and PPTs follow 
from a simple two-part lemma. 


Lemma 4.3. For positive x, y, and z we have 
(4.1) (i) «@+yt+z)* =(«4+2)? 4+ (vy +2)" © 2 = 2xy, 
(4.2) (ii) «+ytz)? =2x4+z2vtaeoexrty? =z". 


Proof. In Figure h.6lwe see two possible placements of overlapping car- 
pets in a square room. 
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FIGURE 4.6 


Applying the carpets theorem to Figure 4.6(a)| yields (4..1)), and ap- 
H.6(b) 


plying the carpets theorem to Figure yields (4.2). 


We now have a second method to generate PTs and PPTs in the fol- 
lowing theorem. 


Theorem 4.4. There exists a one-to-one correspondence between PTs 
and factorizations of even squares of the form z* = 2xy. 


Proof (Gomez, 2005]). Let (a,b,c) bea PT, and set x = c—b,y = c—a, 
andz=a+b-—c. Thenx+y+z=c,x+z=a,andy+z=b,and 
yields c? = a? + b? if and only if z* = 2xy. Note that the PT 
(a,b,c) = (x+z,y+z,x +y +2) is primitive if and only if x and y are 
relatively prime. 
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Here are some examples of the use of Theorem hal 


Example 4.1. There are three factorizations of 36 in the form z* = 2xy, 
yielding three PTs: 62 = 2- 1-18 yields the PPT (7,24,25); 62 = 2-2-9 
yields the PPT (8,15,17); and 67 = 2-3-6 yields the PT (9,12,15). Note 
that when x and y are relatively prime, the PT is primitive. 


Example 4.2. Since (2k)* = 2-1- 2k? for every positive integer k, we 
have the PPT (2k +1,2k? + 2k,2k?+2k+1), PPTs of the form 
(a, b,b + 1), i.e., the hypotenuse c is one more than the even leg b, and 
a* = 2b+1. The generators are consecutive integers, i.e., (m,n) = 
(k + 1,k). Furthermore, every odd number greater than 1 appears 
as a leg of some PPT, and if c and d are consecutive positive integers, 
then (c + d,2cd,c? + d*) is a PPT (the one above when c = k and 
d=k-+1). 


In Figure we illustrate a?+b? = (b+ iy for the PPT (a, b, b+1) 
with a” = 2b + 1 (the area of the darker gray L-shaped region). 


b* 


FIGURE 4.7 


Example 4.3. Since (2k)* = 2-2-k? for every odd integer k (so that 
2 and k? are relatively prime) we have the PPT (k? + 2k,2k + 2, 
k?+2k+ 2), PPTs of the form (a,b,a + 2), i.e, the hypotenuse c is 
two more than the odd leg a, and b* = 4a + 4. The generators are 
(m,n) = (k + 1,1) with k odd. 


In Figure /4.7(b)| we illustrate a2 + b? = (a + 2)” for the PPT (a,b, 
a+ 2) with b? = 4a + 4 (the area of the dark gray region bordering the 
light gray square). 


Example 4.4. The factorization n? = 2k(k + 1) is equivalent to 
(n+2k+ 1° =(n+ k)? +(n+k+ i Note that the legs n + k and 
n +k +1 are consecutive integers (which are relatively prime). In this 


4.5. ALMOST ISOSCELES PRIMITIVE PYTHAGOREAN TRIPLES 61 


case we refer to the PPT as almost isosceles. For example, 27=2-1-2 
yields the PPT (3,4,5), and 12” = 2-8-9 yields the PPT (21, 20, 29). 
Are there other almost isosceles PPTs? If so, how many? 


4.5. Almost isosceles primitive Pythagorean triples 
To answer the questions at the end of Example h.4} we have 
Theorem 4.5. Infinitely many almost isosceles PPTs exist. 


We present two proofs, the first using the carpets theorem, the sec- 
ond via one of the Pell equations we solved in Chapter B} 


Proof 1. Let (a,a+1,c) or (a+1,a,c) be an almost isosceles PPT. 
Place carpets with areas (3a + 2c + 1) and (3a + 2c + 2y in aroom 


with area (4a + 3c + 2)”, as shown in Figure 
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FIGURE 4.8 


The carpets theorem yields 
(4a + 3c + 2)? = (3a + 2c +1)? + Bat 2c +2) 
© (2a+c4+1)*? =2(at+c)(atct1). 
Now place two carpets each with area (a+ c)(a+c+1) ina room 


with area (2a +c +1)’, as shown in Figure The carpets the- 
orem yields 


(2a+c+1) =2(ato(atc+1) ea? +(at1) = 
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Combining the two equivalences yields 
(4.3) a2+(at1)? =c? 
@= (3at2c+1) + (Bat 2c+ 2)" = (4at3c4 2)’, 


i.e., each almost isosceles PPT yields a larger one. This generates the se- 
quence (3, 4, 5), (21, 20, 29), (119, 120, 169), (697, 696, 985), etc. 


Proof 2. Let (a,a+1,c) or (a+ 1,a,c) be an almost isosceles PPT. In 
Figure we show that two copies of this triangle (in light gray), an 
isosceles right triangle with legs c (in white), and a right triangle with 
legs 1 and 2a + 1 and hypotenuse cv2 (in darker gray) partition a 
(2a + 1) X (a + 1) rectangle. 
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FIGURE 4.9 


Hence (2a + 1)* + 1 = 2c?, equivalent to (2a + 1)* — 2c? = —1, 
the second Pell equation x? — 2y? = —1 we solved in Section B.4 with 
x =2a+1andy =c. Each solution (%y, y,) forn = 1 in TableB.4/ has x 
oddand y = 5, and hence leads to an almost isosceles PPT. For example, 
(x1, ¥1) = (7,5) leads to the PPT (3,4,5), (x2, y.) = (41, 29) leads to the 
PPT (21,20,29), (x3, ¥3) = (239, 169) leads to the PPT (119,120,169), 
etc. 


The next two corollaries exploit relationships between almost 
isosceles PPTs and triangular numbers. 


Corollary 4.6. There exist infinitely many PTs in which the legs are con- 
secutive triangular numbers. 


Proof. In Figure 4.10] we illustrate with Pythagorean triangles, expand- 
ing the (a,a + 1,c) triangle by a factor of 2a + 1 to form a Pythagorean 
triangle with legs of length T>, and T>g41. 
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? De (2a+1) a(2a+1)=T>, c(2a+1) 
. ‘ ——? 
a+] (2a+1)(a+1)=Tq4} 
FIGURE 4.10 


We note that it is possible to find a PT whose three sides are trian- 
gular numbers, e.g., (1132, 7143, 1164). However, it is unknown whether 
the number of such PTs is finite or infinite. 

In Chapters and Blwe proved that infinitely many square trian- 
gular numbers exist. Here is yet another proof using almost isosceles 
PPTs. 


Corollary 4.7. Infinitely many square triangular numbers exist. 


Proof. In the first proof of Theorem 4.5|we showed that a2 + (a + 1)? = 


c* is equivalent to (2a + c + 1)* =2(at+c)(atct1) = 47 q4¢. Since 


2 
. ct+1 P 
c + 1 is even, we have T,,, = (a + —) . Thus there is a one-to-one 


correspondence between almost isosceles PPTs and square triangular 
numbers. Since there are infinitely many almost isosceles PPTs, there 
are also infinitely many square triangular numbers. 


Thus 32 + 42 = 5? implies Tg = 62, 207 + 217 = 297 implies Tyg = 
357,119? + 120” = 169” implies Tygg = 204”, etc. 

We conclude this section by showing how almost isosceles PPTs 
solve the Pell equations x? — 2y* = +1 encountered in Section B.4, 


Theorem 4.8. If (m,n) is the generator of an almost isosceles PPT 
(a,b,c), where a and b differ by 1, then (x,y) = (m — n,n) is a solu- 
tion to either x* — 2y* = +1 or x* — 2y? = -1. 

Proof. Ifa — b = +1, then (m2 —n?) — 2mn = +1 or, equivalently, 
(m—n)* — 2n? = 41. 


It remains to show how to find the generators of the almost isosceles 
PPTs. Setting 2a + 1=a+(a+1) = a+b in the hypotenuse term 
4a + 3c + 2 in (4.3) yields 


4a+3c+2=2a+2b+3c =5m2+4mnt+n? = (2mtn)* +m, 


so if (m,n) generates one almost isosceles PPT, (2m + n,m) generates 
the next one. Hence the sequence of generators is (m,n) = (2,1), (5,2), 
(12,5), (29,12), and so on, yielding the sequence of solutions (x,y) = 
(1,1), (3,2), (7,5), (17,12), etc., to x? — 2y? = +1 (the signs alternate, 
beginning with —). 
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4.6. A Pythagorean triple tree 


There are several ways to arrange all the PPTs into a trinary tree, a 
structure where the PPT (3, 4,5) “branches” into three different PPTs, 
and each of those branches into three different PPTs, and so on. Here is 
one example of such a tree (Price, 2011], showing the first four “gener- 
ations”: 


(31,480,481) <= 
(15,112,113) (255,32,257) << 
(29, 420,421) <= 


(175,288, 337)< 
63,16,65 207,224,305)< 
cane = Pe 
= 
13,84 195,28.197 
eeaee) } 25,312 313) a 
(37° 416 425) <= 


(39,80,89) 24 96,265) <= 

(69,260,269) = 

(135,352,377) — 

(3,4,5) wan 48,73) (231,160,281)<= 


(7,24,25) 


(85,132,157) <= 
(51,140,149) = 
(21,20,29) pace 60,109 )—= 


( 
(11,60,61) (145,24, 24,145) <= 


( 
(5,12,13) (35,12,37) (110 120,169) — 


8,53) —== 

(19,180,181) <= 

(9,40,41) (99,20,101) —= 
(17,144,145) = 


FIGURE 4.11 


The tree in Figure is generated by three linear transformations 
of a triple (a, b,c), one, 


U (a,b,c) = (2a—b+c,2a+ 2b + 2c,2a+b + 3c) 
for each “upward” sloping 7 branch, another, 
A(a,b,c) = (2a+b+c,2a—2b+ 2c,2a—b+3c) 
for each “across” — branch, and a third, 
D (a,b,c) = (2a+ b—c,—2a+ 2b + 2c,-—2a+b+3c) 


for each “downward” sloping \ branch. To show that U, A, and D gen- 
erate PTs, we need to show that U (a,b,c), A (a,b,c), and D (a,b,c) are 
PTs whenever (a, b,c) is a PT. The proof employs the carpets theorem. 
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Theorem 4.9. Let a, b, c be positive numbers. Then 
(4.4) 
c2 =a2 +b? & (2a+b+3c) =(2a—b+c) + (2a+2b+ 2c)’, 


(4.5) 
c2 =a? +b? & (2a—b+3c)’ = (2at+b+c) + (2a—2b + 2c)’, 


and 
(4.6) 
c2 =a2+b? © (-2a+b+ 3c) = (2a+b—c)°+(—2a + 2b + 2c)’. 


Consequently U (a,b,c), A (a,b,c), and D (a,b,c) are PTs if and only if 
(a, b,c) is a PT. 


Proof. (i) Place carpets with areas (2a — b+ c)? and (2a + 2b + 2c)? 
in a room with area (2a + b + 3c)’, as shown in Figure 


<— 2a+b+3c ———> <— 2a-—b+3c ———> 


2a 2a - 


c-b 2a 2c+2b c+b 2a 2c—2b 
(a) (b) 


FIGURE 4.12 


The carpets theorem yields 


(2a+b+3c)’ =(2a—b+c) + (2a+ 2b + 2c) 
© (2a)* = 2(2c + 2b)(c — b). 


But (2a)* = 2 (2c + 2b) (c — b) is equivalent to a? = c? — b?, from 
which follows. Thus U (a, b,c) is a PT if and only if (a, b,c) is a PT. 

(ii) Place carpets with areas (2a + b + c)? and (2a —2b + 2c) ina 
room with area (2a — b + 3c)’, as shown in Figure 4.12(b)} The carpets 
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theorem yields 
(2a—b+3c)° = (2a+b+c) + (2a—2b + 2c) 
© (2a)* = 2(2c — 2b)(c +b). 
But (2a)* = 2(2c — 2b) (c + b) is equivalent to a? = c? — b2, from 
which (4..5)) follows. Thus A (a, b,c) is a PT if and only if (a, b,c) is a PT. 


(iii) Place carpets with areas (2a + b — c) and (—2a+ 2b + 2c)" 
in a room with area (—2a + b + 3c)’, as shown in Figure 


— — —) 


4c ta) 
2a+2b—2c a+b-—c 
c-a 


a da 
5 sGhBs ihe 


(a) (b) 


FIGURE 4.13 


The carpets theorem yields 
(-—2a+b+3c)* = (2a+b—c)* + (—2a + 2b + 2c)? 
© (2a + 2b — 2c)? = 2(c — b) (4c — 4a) 
or, equivalently, 
(-—2a+b+3c)* = (2a+b—c)* + (—2a + 2b + 2c)? 
© (a+ b—c)? =2(c—b)(c — a). 


Now place carpets with areas a? and b? ina room with area c?, as shown 
in Figure 4..13(b)| The carpets theorem yields 


c2 =a? +b? @(atb—c) =2(c—b)(c—a), 


and combining the two equivalences yields (4.6). Thus D(a, b,c) isa PT 
if and only if (a, b,c) is a PT. 


See [[Price, 2011|] for proofs that only PPTs appear in the tree and 


that every PPT appears exactly once. For another trinary tree of 
Pythagorean triples, see [[Brousseau, 1972] . 
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4.7. Primitive Pythagorean triples with square sides 


Examination of the tree in Figure shows that some PPTs have 
an odd square leg, e.g., (9,40,41) and (25,312,313); some have an even 
square leg, e.g., (63,16,65) and (77,36,85); and some have a square hy- 
potenuse, e.g., (7,24,25) and (119,120,169). We now show that there 
are infinitely many PPTs of each type, using Lemma 


Theorem 4.10. There are (a) infinitely many PPTs with an odd square 
leg; (b) infinitely many PPTs with an even square leg; and (c) infinitely 
many PPTs with a square hypotenuse. 

Proof. (a) Let (a, b,c) bea PPT (with a odd), and set (u, v) = (c?, b”) in 


2 
Lemma B.4) so that (c? — b2)" + 4b2c2 = (c? + b?) . Hence (a?, 2bc, 
b* + c?) is a PPT with an odd square leg. The PT is primitive since its 
generator is (m,n) = (c,b). 

(b) Let k be a positive integer, and set (u, v) = (4k*,1) in LemmaB.4] 
so that (4k* — 1)° + 16k = (4k4 + 1)”. Hence (4k4 — 1, 4k2, 4k4 + 1) 
is a PPT (since its generator is (m,n) = (2k?,1)). For another family of 
PPTs with an even square leg, set (u,v) = ((2k + 1)*, 2) in Lemma.4 
yielding the PPT ((2k + 1)* — 4,4(2k + 1)”, (2k + 1)* +4). 

(c) Let (a,b,c) be a PPT, and set (u,v) = (b?, a?) in LemmaB.4} so 
that (b? — a?) +4a2b? = (b? + a”). Hence (|b? — a?|,2ab,c”) isa 
PPT with a square hypotenuse since its generator is (m,n) = (b,a) or 
(a, b). 
Corollary 4.11. There are infinitely many solutions to the Diophantine 
equations x* + y? = z? (eg. (x,y,z) = (2,3,5) or (3,40,41)) and 
xn? y= 2° (eg, 067,2) = (7,245): 

See Exercise 4.14 for solutions to the Diophantine equation x* + 

Bi 
y* =z”. 
Are there PTs with two or three square sides? The answer is no; see 


[Sierpinski, 2003)] for a proof. 


4.8. Pythagorean primes and triangular numbers 


If the hypotenuse c in a PPT (a,b,c) is prime, then it is called a 
Pythagorean prime. The Pythagorean primes less than 100 are 5, 13, 
17, 29, 37, 41, 53, 61, 73, 89, and 97. It follows from Euclid’s formula 
that a Pythagorean prime p equals m? + n?, where m and n have op- 
posite parity, hence p = 4N + 1 for some N. We now show that N is 
the sum of two triangular numbers (one of which may be Ty = 0). Al- 
though we won’t prove it here, if the hypotenuse c in a PPT (a,b,c) is 
composite, it must be a product of Pythagorean primes. 
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Theorem 4.12. Letp = 4N+1= (2u)? 120 a 1) u2=i,v=0,be 

a Pythagorean prime. Then N is the sum of two triangular numbers, i.e., 
(1) N = Tysy + Ty-y-1 ifu > v, and (2) N = Tysy + Ty-y, ifu S v. 

Proof. First note that N = u? + v(v + 1). Figure 4.14] shows thatu >v 


implies N = T,4, + Ty-y—1, and Figure shows that u < v implies 
N = Ty+y + Ty-u- 


PAP PPP HH DH 
5S SR OS SS SS 


v+l1 v+l 


u+vV 


FIGURE 4.15 


A converse of Theorem also holds. 


Theorem 4.13. If N is the sum of two triangular numbers, then 4N + 1 
is the sum of two integer squares with opposite parity. 


Proof. Let N = T, + T, withu => v = O. Simple algebra establishes 
4N+1=4(7,+7,)+1=(utvt+ 1)° +(u- vy’, which we illustrate 
in Figure 14 for (u,v) = (8,3). Observe that u+ v +1 and u—v have 
opposite parity. 
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FIGURE 4.16 


Note that 4N + 1 need not be a Pythagorean prime in Theorem /4.13} 
For example, when N = T3 + T, = 16,4N + 1 = 65. However, it can be 
shown that each composite hypotenuse in a PPT is asum of two squares 
with opposite parity, and that Theorem also holds in this case. 

In the following corollary we employ the preceding theorems to gen- 
erate PPTs using triangular numbers. The proof follows from using m = 
u+v+1landn=u-— vin Euclid’s formula. 


Corollary 4.14. [fu +v+1andu-—vwithu > vare relatively prime 
positive integers, then ((2u + 1) (2v+1),4(%, —T,),4(T, + Ty) + 1) 
isa PPT: 


Example 4.5. Setting v = 0 in the above corollary yields the PPT 
(2u + 1,2u? + 2u,2u? + 2u + 1), a member of the family of PPTs of 
the form (a, b, b + 1) from Example Setting v = u — 1 in the corol- 
lary yields the PPT (4u? — 1,4u,4u? + 1), a member of the family of 
PPTs of the form (a, b,a + 2) from Example 44.3, 


The triangular numbers T,, and T, also appear in the odd leg a = 
(2u + 1) (2v + 1) of the PPT in Corollary since Lemmal2.4|implies 
that a? = (87,, + 1)(87, + 1). 


4.9. Divisibility properties 


A cursory examination of the triples in the tree in Figure leads 
one to suspect that in every PT one side is divisible by 3, one side is 
divisible by 4, and one side is divisible by 5. We now use the result in 
Corollary to show that this is indeed true. 


Example 4.6. The divisibility properties mentioned in the preceding 
paragraph may involve all three sides of a PPT, such as in (3,4,5) and 
(63,16,65); involve only two sides, such as in (15,8,17) and (7,24,25); 
or just one side, such as in (11,60,61) and (119,120,169). 
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Theorem 4.15. In every Pythagorean triple, one leg is divisible by 3, one 
leg is divisible by 4, and one of the three sides is divisible by 5. 


Proof. It suffices to prove the theorem for PPTs, since every PT is ei- 
ther a PPT or an integer multiple of one. Consider the PPT (a, b,c) with 
generator (m,n) with m > n, m and n relatively prime and of oppo- 
site parity, and setu = (nm +n-—1)/2 andv = (m—n-—1)/2. Then 
m=u+v+tiandn=u-—vso0 that 


(a,b,c) = (2u + 1)(2v + 1),4(T, — T,),4(T, + T,) + 1) 


as in Corollary Clearly 4|b. If uw = 1 (mod 3), then 2u + 1 
0 (mod 3) so that 3|a. Similarly if v = 1 (mod 3), then 2v+1 = 
(mod 3) so that 3|a. Ifboth u # 1 (mod 3) andv £# 1 (mod 3), then as 
shown in Theorem both T,, and T, are congruent to 0 (mod 3), and 
hence 3|b. 

Ifu = 2 (mod 5), then 2u + 1 = O (mod 5) so that 5a. Similarly 
ifv = 2 (mod 5), then 2v + 1 = 0 (mod 5) so that 5|a. If both u # 2 
(mod 5) and v # 2 (mod 5), then as shown in Theorem each of T,, 
and T, is congruent to 0 or to 1 (mod 5). If both T,, and T, are congruent 
to 0 or both are congruent to 1, then 4 (T,, — T,) = 0 (mod 5) so that 
5|b. If one of T,, or T, is congruent to 0 and the other is congruent to 
1, then 4 (T,, + T,) + 1 = 0 (mod 5) so that 5|c, which completes the 
proof. 


© Ill 


4.10. Pythagorean triangles 


We begin by considering the inradius of a triangle, the radius of its 
incircle (the triangle’s inscribed circle), as illustrated in Figure for 
a right triangle, and show that the inradius r is a linear function of the 
lengths of the three sides. 


r b-r 


FIGURE 4.17 
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Lemma 4.16. The inradius r of a right triangle with legs a and b and 
hypotenuse c isr = =(a b=); 


Proof. See Figure which illustrates how the angle bisectors of the 
triangle and radii of the incircle partition the triangle into six smaller 
triangles, which are congruent in pairs. Consequently the hypotenuse c 
equals a + b — 2r, from which the result follows. 


Theorem 4.17. The inradius of a Pythagorean triangle is an integer. 


Proof. Consider the PA (ka, kb, kc), where (a, b, c) isa PPT and k a posi- 
tive integer. Then its inradius r = “(a+b —c) isan integer sinceat+b—c 


is even. If the generator of (a, b,c) is (m,n), then the inradius r of the 
PA (ka,kb, kc) iskn(m—n). 


Example 4.7. In the proof of Theorem 4.4] we encountered PTs of the 
form (x + z,y+z,x+y+z) which arise from factorizations of even 
squares of the form z? = 2xy. Here the inradius of a PA isr = 2/2 
(an integer since z is even). So itis a simple matter to find all PA’s with 
a given inradius: Factor the square of the diameter z = 2r, e.g., factor- 
ing 62 yields three PA’s with r = 3: (7,24,25), (15,8,17), and (9,12,15). 
See Figure 


12 15 24 


FIGURE 4.18 


In addition to an incircle that is tangent to each of the three sides, 
every triangle possesses three excircles, each tangent to one side and to 
extensions of the other two sides. In Figure we see the three ex- 
circles, with their excenters I,, I,, I., and exradii rz, %, T-, the subscript 
indicating the triangle’s side of tangency to the excircle. We have illus- 
trated the circles for a right triangle, although all triangles possess the 
three circles, centers, and radii. 
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FIGURE 4.19 


Theorem 4.18. The three exradii of a Pythagorean triangle are integers. 


Proof. In Figure we see a general triangle ABC with a portion 
of its excircle with excenter J, and exradius 7,. The area K of ABC is 
equal to the sum of the areas of triangles AJ,C and AI, B minus the area 
of triangle BI,C, hence 


1 1 1 bEcCH=_ 
K= Pa + hats — 7a = a 
and hence my = 2K/(b+c-—a). Similarly 7 = 2K/(a—b-+c) and 
% =2K/(at+b—c). 

Now consider a primitive PA (a, b,c) with generator (m,n), perime- 
terP =a+b+c = 2m(m+n),andarea K = ab/2 = mn(m? — n’?). 
Then 

ab 
~ b+c-a 
Similarly n, = a =n(m+n) andr, = ae = m(m-+n), hence 
all three exradii are integers. In a right triangle the exradius 7, of the 
excircle on the hypotenuse equals the semiperimeter s = P/2 of (a,b,c). 
It now follows that every PA has integer exradii. 


Example 4.8. The PA in Figure 4.19(a)lis the (3,4,5) triangle with gen- 
erator (m,n) = (2,1) and thus (7, 7%, %, 7%) = (1, 2, 3,6). 


See Exercise 4.10 for some remarkable relationships among the in- 
radius r, the three exradii, the perimeter, and the area of a PA. 


Ty =m(m-—n). 
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4.11. Pythagorean runs 


The PPT (3,4,5) is the only triple where the sum of two consecutive 
squares is the next square. But what about sums of three consecutive 
squares? Four consecutive squares? Observe: 

37447 = 5?; 
10° + 11° +12 = 137 + 14°; 
21° + 227 + 23° + 247 = 25° + 26° + 27°; etc. 

What’s the pattern? A little reflection yields the observation that each 
number squared immediately to the left of the equals sign is four times 
a triangular number, ie., 4 = 4(1), 12 = 4(1 + 2), 24 = 4(1 + 2 +3), 
etc., which leads to the identity 

(4.7) (47, —n)° +--+ (4T)? = (4T +1) +2 + 4T, +0)”. 


The identity can be easily verified by mathematical induction, but Fig- 
ures and for the case n = 3) may better explain why the re- 


lationship holds [JBoardman, 2000]. In Figure we use the fact that 
4T, = 4-1+4-2+4-3 to dissect 4T3 into three sets of vertical strips. 


S, 
i 
++ 
-- 


FIGURE 4.20 


In Figure we use the strips to create “borders” around the re- 
maining three squares on the left side of the equals sign in to cre- 
ate the three squares on the right side of the equals sign, yielding alge 
227 + 237 +24” = 25" + 267 + 27”. 


4.12. Sums of two squares 


The problem of representing an integer as the sum of two squares 
has a long history, dating back to at least the time of Pythagoras. In 
Section 4.8) we discussed Pythagorean primes—primes p = 1 (mod 4) 
that can be written as the sum of two squares. It is easy to show that no 
integer n = 3 (mod 4) can be so represented. More generally we can 
ask: Which integers n can be represented as the sum of two squares and 
how many representations are there? 
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os re = 25° 
ml | = 26 
217 He - 27° 


FIGURE 4.21 


Let 7, (n) denote the number of ways to represent an integer n as the 
sum of two squares (positive, negative, or zero), i.e., %(n) denotes the 
number of solutions in integers (x, y) to the equation x? + y* = n. For 
example, % (13) = 8since the solutions to x?+y? = 13 are (2, 3), (3, 2), 
(—2, 3), (3, —2), (2, —3), (—3, 2), (—2, —3), and (—3, —2), as illustrated 
in Figure Since rz (n) = 0 whenever n = 3 (mod 4), %)(n) is a very 
erratic function, as illustrated in Figure 4.22] for 0<n< 30. 

But we can ask: What is the average value of %(k) for1 < k < n? 
We define R,(n) to be the number of solutions in integers to x?+y? < n, 
and then the average of 7,(k) forl1 < k < nis 


mA)+nZ)t+-+m™M) _ RM) 
n a 
Computation of N, (n) and N, (n)/n yields Table 


TABLE 4.1 


n |1 2 3 4 5 10 20 50 100 
Ron) |5 9 9 13 21 37 69 161 317 
R,(n)/n|5 4.5 3 3.25 4.2 3.7 345 3.22 3.17 
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15 


FIGURE 4.22 


The average value of 72 over all the positive integers is the limit of 
R2(n)/nasn — ov if the limit exists. It does, and we have 


Theorem 4.19. lim,_,.. a 1 


Proof. The proof is based on a geometric interpretation of Rz (n) and is 
due to Gauss. R>2 (n) is the number of lattice points (points with integer 
coordinates) in or on the circle x? + y* = n. For example, R2 (10) = 
37 since the circle centered at the origin with radius V10 contains 37 
lattice points, as illustrated in Figure If we draw a square with 
area 1 centered at each of the 37 lattice points, then the total area of the 
squares (in gray) is also R2(10). Thus we would expect the area of the 


FIGURE 4.23 
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squares to be approximately the area of the circle, or in general, R2(n) 


to be approximately n(n) = 71N. 

If we expand the circle with radius /n by one-half the length of the 
diagonal (V2/2) of a square with area 1, then the expanded circle con- 
tains all the squares. If we contract the circle by the same amount, then 
the contracted circle is contained in the union of all the squares, as il- 
lustrated in Figure Thus 


n(n —V2/2) < Ry (n) <n(vn + v2/2) 


Dividing each term by n and applying the squeeze theorem for limits 
yields the desired result. 


2 


We now return to the first question at the end of the introductory 
paragraph in this section: Which integers can be represented as the sum 
of two squares? In Table 1.2] we present the prime power factorizations 
for the positive integers considered in Figure 4.22] in two groups, those 
for which 7, (k) > 0 and those for which 7, (k) = 0. 


TABLE 4.2 
%(k) > 0 %(k) = 0 
1=1 16=27 3=3 21=3:-7 
2=2 17=17 6=2-3 22=2-11 


A= 2? 18=2-37 | 7=7 23 = 23 
B=5 2022765 | 11] 411 24=2°.3 
g=2° as Wears 273° 
G29? 262213 | 14]2*7 282277 
10=10 29=29 15=3:5 30=2:-3:°-5 
13 = 13 19=19 


The data in Table 4.2) leads to the statement of the following theo- 
rem, whose proof can be found in most number theory textbooks. 


Theorem 4.20. A positive integer is not the sum of two squares if and 
only if each of its prime factors congruent to 3 (mod 4) appears to an 
odd power. 


4.13. Pythagorean quadruples and Pythagorean boxes 


A Pythagorean quadruple (PQ) is a quadruple (a, b,c, d) of positive 
integers such that a? + b? + c* = d?. PQs can easily be constructed 
from PTs. For example, the PTs (3,4,5) and (5,12,13) yield the PQ 
(3,4, 12,13). But this procedure only generates a few PQs. We present 
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an extension of Euclid’s formula that generates all the PQs [Sierpinski, 
2003]. But first, we need the following lemma, from Problem 19 in Book 


III of Arithmetica, written by Diophantus of Alexandria. It was in the 
margin of one page in his copy of this book that Pierre de Fermat stated 
without proof his famous “Last Theorem.” The Arithmetica was trans- 
lated from Greek to Latin by Claude Gaspard Bachet de Méziriac in 1621. 
Below is the cover of an edition with commentary by Fermat published 
by Fermat’s son Clément-Samuel de Fermat in 1670. 
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Lemma 4.21. Diophantus of Alexandria’s sum of squares identity. If 
two positive integers are each sums of two squares, then their product is 
a sum of two squares in two different ways, i.e., 


(4.8) (a? + b?) (c?2 + d?) = (ac + bd)’ + (ad — bc)” 
and 
(4.9) (a? + b?) (c?2 +d?) = (ad + bc)’ + (ac — bd)’. 


Proof. See Figure for a proof of for the case ad > bc (the 

case ad < bc is similar), using the Pythagorean theorem to show that 
2 

(va? + b2Vc2 + a?) = (ac + bd) + (ad — be)’. Exchanging c and d 

in the figure yields a proof of (4.9). 


We now state and prove an identity to generate PQs. 


Theorem 4.22. For integers m, n, p, and q, 


(m2 + n? — p* — q2) + (2mp + 2nq)* + (2mq — 2np)* 
= (m? +n? + p? +42). 
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ac+bd 


ad—bc 


be 


FIGURE 4.24 


Proof. Setting u = m2 + n? and v = p? + q? in Lemma.4|yields 
(4.10) (m2 +n? — p? —q2) +4 (m? + n?) (p? + q?) 
= (m2 +n? +p? + q?) - 
Applying to (m? + n?) (p? + q?) and multiplying by 4 yields 
4.(m? +n?) (p? + q2) = (2mp + 2nq)* + (2mq — 2np)’, 


which completes the proof. See [Spira, 1962) for a proof that this iden- 
tity generates all PQs. 


A Pythagorean box, a three-dimensional analogue of a Pythagorean 
triangle, is arectangular box whose edges and interior diagonal are pos- 
itive integers, and hence it can be represented by a PQ (a,b,c,d). A 
nearly cubic Pythagorean box is one given by (a,a,a + 1,d) or 
(a,a+1,a+1,d). Examples include (1,2,2,3) and (6,6,7,11), as well 
as the degenerate box (0,0,1,1). We begin with a lemma relating one 
Pythagorean box to a larger one. 


Lemma 4.23. a? + b? +c? = d? if and only if 
(at+b+d)°+(b+c+d)+(atc+d) =(atbt+ct+2d)’. 


Proof. In this proof we place equilateral triangular carpets in an 
equilateral triangular room. Place carpets with areas (a+ b+ d)’, 
(Deeb dy, and (at+ct d)* in aroom with area (a+b+c+t+ 2d)’, 
as shown in Figure (here each small A and V triangle has area 
1, and we recall that each large triangle contains a square number of 
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+ q+ b+d——_> 


~« a+b+c+2d ————___—__» 


FIGURE 4.25 


small A and V triangles). The result now follows from the carpets the- 
orem. 


Theorem 4.24. Infinitely many nearly cubic Pythagorean boxes exist. 


Proof. Repeated applications of Lemma 4.23] yield the following infinite 
sequence of PQs representing Pythagorean boxes: 


(0,0,1,1) > (4,2, 2,3) > (6,6, 7, 11) 
— (23,24, 24,41) > (88, 88, 89, 153) > +. 


4.14. Exercises 


4.1 Show that if (a, b,c) isa PPT, then a and b have opposite parity and 
c is odd. 

4.2 Show that the even leg of a PPT is always a multiple of 4, and that 
every multiple of 4 appears in some PPT. 

4.3 Show that if (a,b,c) is a PPT, then 2c is a sum of two odd 
squares. 

4.4 The PTs (6,8, 10) and (27, 36,45) share the interesting property 
that one leg is an integral cube while the other leg and the hy- 
potenuse are consecutive triangular numbers. Are there other PTs 
with this property? If so, are any primitive? 
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4.6 


4.7 


4.8 


4.9 


4.10 


4,11 


4.12 


4.13 
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In Exercise 2.8 we encountered the Fermat numbers, integers of 
the form f, = 22" +1 forn > 0. Show that for n > 1, every Fermat 
number appears as the hypotenuse of some PPT. 

Show that T,, = k(k + 1) is an oblong triangular number if and 
only if (m,m + 1,2k + 1) is an almost isosceles PPT. 

Show that there are always two different PPTs with a given odd 
prime number as inradius. [Hint: Solve p = n(m — n) for m and 
n.] 

Show that there exist infinitely many PPTs whose area is the prod- 
uct of three consecutive integers. 

In Chapter {1) we encountered the triangular numbers {7T;,} = 
{1, 3, 6,10, ...} as sums of integers and the pyramidal numbers 
{Pyr,} = {1,5,14, 30, ...} as sums of squares of integers. These 
numbers appear in the perimeter P and area A of some primitive 
PA’s, as Table illustrates. 


TABLE 4.3 
PPT | P | A 
(3,4,5) | 12 = 2-T3| 6 = 6-Pyr, 
(5,12,13) | 30 = 2-T, | 30 = 6- Pyr, 
(7,24,25) |56 = 2-T, | 84 = 6-Pyr, 


Are there other PA’s with this property? 

Prove the following identities for the inradius r, the three 
exradii 7%, %, %, the perimeter P, semiperimeter s, and the area 
k of a primitive PA: 

(a) r+y EH = Te 

(b) r+m+%+%7 =P, 

(c) 1 = i + i + > 

Tr Tq Th Tc 

(d) t%% = 77% = K and hence rm%pr_ = K?, 

(e) 44% 4+hE Pie = 5. 

Prove that there are infinitely many Pythagorean boxes whose 
edges are integer squares. [Hint: Let (a,b,c) be a PT and show 
that the Pythagorean box with edges (ab)’, (bc), and (ac)* has 
interior diagonal c* — a7b?.] 

Prove that if an integer is the sum of three positive integer 
squares, then so is its square. 

Three squares in arithmetic progression. (1, 25, 49) and (289, 625, 
961) are two examples of three squares in arithmetic progression. 
Are there others? [Hint: See Figure 4.26] 
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FIGURE 4.26 


4.14 The Diophantine equation x* +y* = z3. Sometimes the sum of two 
squares can be a cube, e.g, 52 = 10° + 52 = 2? + 11 and 13° = 
46° + 92 = 39° + 26”. Are there other solutions to x2 + y2 = z3? 


CHAPTER 5 


Irrational Numbers 


To attempt to apply rational arithmetic to a problem in geometry 
resulted in the first crisis in the history of mathematics. The two 
relatively simple problems—the determination of the diagonal of 
a square and that of the circumference of a circle—revealed the 
existence of new mathematical beings for which no place could be 
found within the rational domain. 


Tobias Dantzig 


While the primary objects of study in an elementary number the- 
ory course are the integers, real numbers formed from integers are of- 
ten considered as well. The rational numbers are, as the name sug- 
gests, ratios of integers. Real numbers that are not rational—the irra- 
tional numbers—are perhaps of more interest than the rational num- 
bers. This may follow from the definition of an irrational number, telling 
us what it is not (an irrational number is a real number that is not ra- 
tional) rather than what it is. 

We confine ourselves to the positive quadratic irrational numbers, 
those that are positive irrational roots of quadratic equations with in- 
teger coefficients. We begin with some visual demonstrations of the ir- 
rationality of simple quadratic irrationals such as V2, V3, V5, and the 
golden ratio. We illustrate rational approximations to these numbers 
using solutions to Pell equations and certain triangles with integer 
sides, and express some quadratic irrationals in terms of continued frac- 
tions. 


5.1. The irrationality of V2 


The discovery that V2 is irrational was probably due to the Pythago- 
reans, perhaps even Pythagoras himself. This may have occurred by 
observing that the hypotenuse of an isosceles right triangle is not com- 
mensurable with the legs (two non-zero real numbers are commensu- 
rable if and only if their ratio is rational). Since “irrational” means “not 
rational,” virtually all proofs of irrationality are indirect, showing that 
the assumption that the number is rational leads to a contradiction. 


83 
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Theorem 5.1. V2 is irrational. 

Proof 1 ([Apostol, 2000). If V2 were rational, then V2 = a/b in lowest 
terms for positive integers a and b. Then 2b? = a?, so thata right trian- 
gle with legs b and hypotenuse a is the smallest isosceles right triangle 
with integer sides. However (see Figure 5.1), 


S 
A 
N 
b a 
D 
aR 
~~) 
b b 
if this is an isosceles right triangle then there is a smaller one with the 
with integer sides, same property. 


FIGURE 5.1 


Indeed, if the sides of the larger triangle above are b, b, and a, then 
the sides of the smaller gray triangle on the right are positive integers 
a —b,a-—b,and 2b — a. Thus the assumption that V2 = a/b in lowest 
terms is false, and V2 is irrational. 


The above proof appears in (Bloom, 1995] as a one-sentence proof 
of the irrationality of V2: If V2 were rational, say V2 = a/b in lowest 
terms, then also V2 = (2b — a)/(a — b) in lower terms, giving a contra- 
diction. 

Here is a second proof, created by the American mathematician 
Stanley Tennenbaum (1927-2005) in the 1950s [|Conway, 2005] . Item- 
ploys the carpets theorem from Section 


Proof 2. Assume that V2 is rational, and write V2 = a/b in lowest 
terms for positive integers a and b. Then a? = 2b? = b? + b?. So place 
two b x b carpets on the floor in an a X a room, as shown in Figure 

By the carpets theorem, the area (2b — a)’ of the overlap (in dark 
gray) equals the combined area 2(a — by? of the uncovered floor (in 
white). But these squares have positive integer sides with 2b — a and 
a — b smaller than a and b, respectively, a contradiction. Hence V2 is 
irrational. 
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a+_————_ a ——_ 


FIGURE 5.2 


Hippasus and the irrationality of V2 


Proving theorems in mathematics is not usually dangerous. Indeed, it is of- 
ten quite rewarding, but in Pythagoras’ time, the situation may have been 
different. Legend has it that Hippasus of Metapontum (c. 500 BCE) dis- 
covered the unexpected irrationality of V2, explained his result, and con- 
sequently was thrown overboard at sea by his fellow Pythagoreans. The 
discovery had destroyed the ideal of commensurability in geometry. For- 
tunately, history records no other such deaths as a result of proving theo- 
rems. 


We conclude this section with four examples of rational and irra- 
tional numbers constructed from V2. 


Example 5.1. An irrational number to an irrational power may be ratio- 


V2 
nal. To prove this statement, we need only give an example. If V2 _ is 


V2 V2 vR 
rational, that is our example. IfV2. is irrational, then (v2 = 2is 


our example. 


Example 5.2. An irrational number to an irrational power may be ir- 


V2 V2 
rational. If V2 _ is irrational, that is our example. If V2. is rational, 


V2+1 v2 
then V2 nae V2  -V2is our example (since the sum or product of a 


rational number and an irrational number is irrational). 


2 
An irrational number to a rational power may be rational (V2 ) or 


1 
irrational (V2 ). The case of arational number to an irrational power is 
examined in Exercise 5.2. 
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We proved the statements in Examples §.1|and 5.2| without knowing 
V2 
whether V2 _ is rational or irrational. We used only the fact that V2 


is irrational. In fact, io" is irrational, as it is the square root of the 
Gelfond-Schneider constant 2v2, which is known to be transcendental. 
A (possibly complex) number is transcendental if it is not algebraic; a 
number is algebraic if it is the root ofa non-zero polynomial with integer 


coefficients; all (real) transcendental numbers are irrational. However, 
V2 


2 
it is unknown whether V2 is rational or irrational. 


Hilbert and the Gelfond-Schneider constant 2¥2 


Early in Chapter Bl we mentioned Hilbert and the address with 23 prob- 
lems that he gave at the Second International Congress of Mathematicians 
in Paris in 1900. In his seventh problem, he conjectured that “the expres- 
sion a, for algebraic base a (not equal to 0 or 1) and an irrational alge- 
braic exponent D, e.g., 2v2 always represents a transcendental or at least 
an irrational number.” The conjecture was proved by A. O. Gelfond and T. 
Schneider independently in 1934. However, it is still unknown whether 
a” is transcendental when both a and b are transcendental. 


Example 5.3. A tower of square roots of 2. What is the meaning of an 
expression such as 


« 


The ellipsis “.’” indicates infinitely many copies of V2 as exponents, 
so we consider the limit of a sequence. Set x, = V2, define xX,4, = 


V2" for n > 1, and then the question becomes: does the sequence 
x, } converge, and, if so, to what? The sequence converges since it is 
n g q 8 


increasing (the function f (x) = 2 is increasing) and bounded above 
by 2 (since V2 < 2 we have 


woe” ea er od) 


Hence the sequence {x,,} has a limit L < 2, and 


limpsoo Xn L 


L = lim x_4, = lim V2" = V2 = 2. 
n—-0co n-c 
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FIGURE 5.3 


By inspection (see Figure the equation L = Wa has solutions L = 2 
and L = 4. Thus the limit is 2 since we know that L < 2. 


Example 5.4. Francois Viéte’s remarkable formula for m. In 1583 the 
French mathematician Francois Viéte (1540-1603) discovered the fol- 
lowing infinite product for 2 employing only 2 and v2: 


; See bese prvi 
ee 


1A 2 


We prove this formula as Viéte did, with a geometric argument based on 
ratios of the areas of regular polygons inscribed in the unit circle. First 
notice that each term in the infinite product is a cosine obtained from 
cos 1/4 = 2/2 by repeated use of the half-angle formula cos(@/2) = 
(1+ cos @)/2 = <2 + 2 cos @. It follows that the terms in the prod- 
uct are given by cos(z/2") forn = 2. 

Let V, = cos(m/4) - cos(7/8) --- cos(z/2") forn = 2. To show that 
the limit of V, as n tends to infinity is 2/7, we inscribe regular poly- 
gons with 2” sides in the unit circle. It is easy to show that the area 
of a regular polygon with k sides inscribed in the unit circle has area 
(k/2) sin(2m/k). So if A, denotes the area of a polygon with 2” sides, 
then A, = 2-1 sin(/2"~+) for n = 2. Hence the ratio of the areas of 
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two successive such polygons is 
An 2° sin(a/2"*) 
Ansa 2 sin(ae/2") 
2” sin(/2") cos(1/2") 
OR sin(e/2") 
= cos(m/2") 


forn => 2. Thus 
=> Az Az An-1 . An > Az 
. A3 Ay An An+1 An+1 
so that limy_.~ V, = 2/m since A, = 2 (the area of a square inscribed in 
the unit circle) and limy,_,.. An41 = 7 (the limiting area of the inscribed 
polygons is the area of the unit circle). 


See Exercise 5.7 for a related result involving V2. 


5.2. Rational approximations to V2: Pell equations 


Tennenbaum’s proof that V2 is irrational (proof 2 of Theorem 5.1) 
can be modified to provide a procedure for generating rational approxi- 
mations to V2. The procedure is very similar to the approach to solving 
Pell’s equations in Chapter 3} We set x = 2b —-aandy = a—bin 
Figure in the preceding section and employ the carpets theorem to 
prove the next theorem. 


Theorem 5.2. Let x and y be positive integers. Then 
(5.1) (x +2y)? —2(x + y)? = +1 ifand only if x? — 2y? = F1. 
Proof. See Figure and compute the area of the large square in two 
ways to yield (x + 2y)* = 2(x + y)* — x2 + 2y2. Hence 

x? — 2y? = 2(x + y)* — (x + 2y)?, 


from which the desired result follows. 


Set 
(5.2) Xn41 = Xn + 2Vy and Yn+1 = Xn + Yn 


in (5.1) yielding x2,, — 2y2,, = +1 ifand only if x2 — 2y? = #1. Now 
set (1,1) = (1,1), and hence x2 — 2y?2 = (—1)” foralln > 1. Division 
by y? yields for alln > 1, 


2 n 
(5.3) (=) ee 
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x+y y 


+ x+ 2y ———> 


FIGURE 5.4 


Thus x, /¥p will approximate V2 (larger than V2 when nis even, smaller 
when nis odd) with an absolute error less than 1/2y,7. This error bound 


follows from 
1 x. 
a (= oe v2) 


yr \Yn 
since ] + V2 >1+-2 > 2. See Table for the numerators and 


Xx. 
aoe) 


Yn 


x. 
eae) 


Yn 


ears 


denominators of the first few approximations to V2 generated by The- 


orem 5.2|with (1,1) = (1,1). 


TABLE 5.1 
n | 1 23 4 5 6 7 
MeL ea 7 «AZ AL 99° 239 
Y,\1 2 5 12 29 70 169 


Example 5.5. A geometric interpretation of the error. A triangle 
(Vy Yn» Xn) with side lengths from Table is isosceles and almost a 
right triangle since x2 = yZ + yZ + (—1)". If a, denotes the size of 
the angle opposite the longest side x,, then the law of cosines yields 
x2 = 2y? — 2y2 cosa, so that cosa, = —(—1)"/(2y,). Thus the error 


term (—1)"/y; in (5.3) is —2 cos ay. 


5.3. Rational approximations to V2: Almost isosceles PPTs 


In Theorem we encountered a relationship between almost 
isosceles PPTs and the Pell equations x? — 2y? = +1. Geometrically, 
the right triangle whose sides are the elements of an almost isosceles 
PPT is “almost similar” to an isosceles right triangle. Consequently, we 
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can use ratios of functions of the sides of such PPTs to approximate v2, 
via the following theorem. 


Theorem 5.3. Let (a,b,c) with |a — b| = 1 be an almost isosceles PPT. 
Then (x,y) = (a + b,c) is a solution to x* — 2y? = —1 and (x,y) = 
(a+b+2c,a+b +c) isa solution to x* — 2y? = +1. Hence (a + b)/c 
and (a+ b + 2c)/(a + b +c) are rational approximations to V2. 


Proof. We consider the case b = a + 1; the case a = b + 1 is similar. If 
(a, b,c) isan almost isosceles PA, then computing the area of the square 


in Figure .5|in two ways yields c? = 4- sa(a +1)+1=2a(a+1)+1. 


FIGURE 5.5 


Completing the square produces (2a+ 1) —2c* = -1,or (a+b)* = 
2c? = —1. Thus (x,y) = (a+ b,c) is a solution to x — 2y* = —1, and 
consequently (a + b)/c is an approximation to V2. From we con- 
clude that (x, y) = (a+b+2c,a+b+c) isa solution to x? —2y? = +1, 
and consequently (a + b + 2c)/(a + b + c) is also an approximation to 


V2: 


In Section 4.5} we encountered the first four almost isosceles PPTs: 
(3, 4,5), (21, 20, 29), (119, 120, 169), and (697, 696,985). These yield 
the following approximations to V2: From (a+b)/c we have 7/5, 41/29, 
239/169, and 1393/985 (see Table aan and from (a+ b+ 2c)/ 
(a+b +c) we have 17/12, 99/70, 577/408, and 3363/2378 (see Ta- 
ble B.3). The numerators and denominators of these approximations 
also appear in Table 


5.4. The irrationality of V3 and V5 


We mimic the two proofs of the irrationality of V2 to give two proofs 
that V3 is irrational. 
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Theorem 5.4. V3 is irrational. 


Proof 1. Assume that V3 is rational, and write V3 = a/b in lowest terms 
for positive integers a and b. Then a? = 3b? or, equivalently, a? + b? = 
(2b)°. Now draw a right triangle with hypotenuse 2b and legs a and b, 
as shown in Figure Ed 
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FIGURE 5.6 


We now construct the smaller right triangle shaded gray with in- 
teger hypotenuse 4b — 2a and integer legs 2b — a and 2a — 3b, so 
that (2a — 3b) = (4b — 2a)’ — (2b—a)* = 3(2b—a)’. Thus V3 = 
(2a — 3b)/(2b — a), acontradiction since 2a—3b and 2b—a are smaller 
than a and b, respectively. Hence V3 is irrational. 


Proof 2. Assume that V3 is rational, and write V3 = a/b in lowest 
terms for positive integers a and b. Then a? = 3b? = b? + b? + b?. 
We place three equilateral triangular carpets with side b in an equilat- 
eral triangular room with side a, as shown in Figure 

The carpets overlap in three equilateral triangles with side 2b — a 
(shaded dark gray) and leave uncovered an equilateral triangle portion 
of the floor with side 2a — 3b (in white). By the carpets theorem, the 
combined area 3k(2b — a)’ equals the area k(2a — 3b) of the floor 
(k = V3/4 whether V3 is rational or not). Thus V3 = (2a — 3b)/ 
(2b — a), a contradiction since 2a — 3b and 2b — a are smaller than 
a and b, respectively. Hence V3 is irrational. 


Example 5.6. Rational approximations to V3. Analogous to what we did 
in Sections 5.2| and |5.3| to approximate V2, we can use solutions to the 
Pell equation x* — 3y* = 1 to approximate V3. If (x»,y,) is a solution 
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FIGURE 5.7 


to x2 — 3y? = 1 (see Section 3.5), then when n > 1 we have 


(3) 395 
Yn Yn 


Thus x,,/y,, yields an approximation to V3 (larger than V3) with an ab- 
solute error less than 1/3y,. The error estimate follows from 


F(eoayenajena 
Yn Yn Yn Yn 


since - + ¥3 > ¥3+V3 > 3. Inthis case we do not obtain approxima- 
tions less than V3, since the Pell equation x? — 3y? = —1 has no integer 
solutions (to see why, examine both sides of the equation modulo 3 and 
recall Example B.2). 

Analogous to our use of almost isosceles PPTs to approximate V2 in 
the previous section, we can use almost equilateral Heronian triangles 
from Example B.10|to approximate 73. In that example we showed that 
if (Xn, Yn) isa solution to x? —3y? = 1, then (2x, —1, 2x, 2x, +1) isan 
almost equilateral Heronian triangle with semiperimeter s = 3x, and 
altitude (to the even side 2x,,) h = 3y,. Thus the ratio of semiperimeter 
to altitude is x,,/y, and approximates V3. 


We conclude this section with a single proof that V5 is irrational. 
Theorem 5.5. V5 is irrational. 


Proof. Assume that V5 is rational, and write V5 = a/b in lowest terms 
for positive integers a and b. Then a? = 5b?, or equivalently, a2 = 
b? + (2b). Now draw a right triangle with hypotenuse a and legs b and 
2b, as shown in Figure 
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FIGURE 5.8 


We now construct the smaller right triangle shaded gray with in- 
teger hypotenuse 5b — 2a and integer legs a — 2b and 2a — 4b, so 
that (5b — 2a)? = (a—2b)* + (2a— 4b)” = 5(2b—a)’. Thus V5 = 
(5b — 2a)/(a — 2b), acontradiction since 5b—2a and a—2b are smaller 
than a and b, respectively. Hence V5 is irrational. 


For another proof of the irrationality of V5 utilizing five pentagonal 
carpets in a pentagonal room, see [Miller and Montague, 2012]. See Ex- 
ample 6.12) for rational approximations to V5 using Fibonacci and Lucas 
numbers. 


5.5. The irrationality of Vd for non-square d 


We now modify the first proof of the irrationality of V2 to show that 
Vd is irrational when d is not the square of an integer. In the proof we 
interpret Vd as the slope of a line through the origin, as shown in Fig- 
ure 


Theorem 5.6. [fd is not the square of an integer, then Vd is irrational. 


Proof. Assume Vd = a/b in lowest terms. Then the point on the line 
y = Vdx = (a/b) x closest to the origin with integer coordinates is 
(b, a). Let k be the greatest integer less than Vd, so thatk < Vd < k+1. 
Then the point with integer coordinates (a—kb, db—ka) lies on the line 
and is closer to the origin since (a/b) (a — kb) = a*/b—ka = db—ka, 
andk <a/b <k+1implies0 < a—kb < band0O < db—ka < a. Thus 
we have a contradiction, and Vd is irrational. When d is the square of 
an integer, i.e., d = k?, then a/b =k/1anda—kb =0=db-—ka. 
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y y=Vdx=Fx 


FIGURE 5.9 


5.6. The golden ratio and the golden rectangle 


The golden ratio (also known as the golden section, the golden mean, 
and the divine proportion) is a ratio that occurs in a variety of mathemat- 
ical settings. Euclid provides a definition and a construction in Book VI 
of the Elements. Paraphrasing Euclid, we say that a line segment AB 
is divided by a point C into the golden ratio if AB/AC = AC/BC, or 


A as 7 See Figure 


a 


at+b 


a 


A Cc B a b 


FIGURE 5.10 


Using the segment AB as the base of a rectangle with height equal to 
AC yields a golden rectangle (see Figure E-10(b), a rectangle that has 
the property that when a square is cut off, the remaining rectangle (in 
gray) is similar to the original. 
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Using @ to denote the golden ratio, we have 


eee eee eee 
=p g T1tGtlts 


so that p* = ~ + 1, and thus g is the positive root of x? — x —1 = 0. 

To show that g is irrational, we need only solve x* — x — 1 = 0 for 
its positive root using the quadratic formula: g = (1 + V5)/2. But it is 
not necessary to solve the quadratic equation, as the following theorem 
shows. 


Theorem 5.7. ¢ is irrational. 


Proof. Assume @¢ is rational, and write g = a/b in lowest terms for pos- 


ee b. 2 

itive integers a and b. Then @ = = — implies thata(a—b) =a — 
b as 

ab = b”. Thus Q= f= —., acontradiction, since b and a — b are pos- 


b a—b 
itive integers smaller than a and b, respectively. Hence ¢ is irrational. 


See Figure 5.11] 


at+b 


b 
a 


a-b 


a b 


FIGURE 5.11 


In addition, using the quadratic formula to obtain g = (1 + V5)/2 
from g = 1+ * is not necessary either. In the following theorem we 
dissect a square with area 5 to obtain the result. 


Theorem 5.8. gp =1+ - with p > 0 implies p = ae 


2 


Proof. See Figure 5.12|to see that (29 — iy = 5, from which the result 
follows. 
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~ lip 
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FIGURE 5.12 


5.7. The golden ratio and the regular pentagon 


The golden ratio g may have been the second number shown to be 
irrational by the Pythagoreans—not from a golden rectangle, but from 
a regular pentagon. Since the diagonal of a square with side 1 is the 
irrational number V2, it is natural to also consider the length of the di- 
agonal of a regular pentagon arr a 1. We now show that the length 
of this diagonal is @. See Figure 5.13} 5.13 


AG 


FIGURE 5.13 


Let x denote the length of the diagonal of a regular pentagon with 
side 1, as shown in Figure 5.14(a)l Two diagonals from a common ver- 
tex and the opposite side form an isosceles triangle (sometimes called 
a golden triangle), which can be partitioned into two smaller triangles, 
as illustrated in Figure 

The angles labeled a in Figure 5.14(b)|are equal since each subtends 
an arc equal to 1/5 the circumference of the circumscribing circle. Thus 
the striped triangle is also isosceles, 6 = 2a, and so a = 36° since the 
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FIGURE 5.14 


sum of the angles in the original triangle yields 5a = 180°. It now fol- 
lows that the unlabeled angle in the gray triangle is 6 = 72°, and hence 
the gray triangle is also isosceles and similar to the original triangle. 
Thus x/1 = 1/(x — 1), and so x is the positive root of x? — x —1 = 0. 
Thus x = gas claimed. 

The regular pentagon can also be used to give another proof of The- 
orem 5.7 (q is irrational). Assume that ¢ is rational and write g = a/b 


in lowest terms for positive integers a and b. Draw a regular pentagon 
with sides b and diagonals a as shown in Figure 


ZBS 
x 


ef 
(a) (b) 


FIGURE 5.15 


a—b 


As shown in Figure |5.14(b)} the gray triangle in Figure 5.15(b)) is 
B15] 


similar to the gray triangle in Figure 5.15(a), and hence g = b/(a — b), 
a contradiction, since b and a — b are positive integers smaller than a 
and b, respectively. Hence ¢ is irrational. 
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5.8. Periodic continued fractions 


An infinite simple continued fraction is an expression of the form 


1 
[g, Az, Az, Az, ...] = Ag + — —+{_- 


a, + 
az + 


a3+ * 
where dg is an integer and Qj, 2, Q3,... are positive integers. When a 
block of a;,’s repeats over and over, the continued fraction is periodic. 
For example, 


1 
[a] = [a,a, a, | =a + —_— 
at+ I 
at+ 
at+*. 
and ; 
[a, b] = [a,b,a,b,...)=at : 
b+ 
oT ee 


are periodic (the vinculum or overbar indicates the repeating block of 
integers) with periods 1 and 2, respectively. 

In most number theory texts you can find theorems that state that 
every infinite simple continued fraction represents an irrational num- 
ber, and those that are periodic represent positive quadratic irration- 
als—numbers of the form p + qvk, where p and q # 0 are rational and 
k is a non-square integer (i.e, p + qvk is a positive irrational root of a 
quadratic equation with integer coefficients). 

We now evaluate the irrational number represented by the periodic 
continued fraction [ a,b ] (note that [ a] with period 1 is a special case 


of [ a, b | with period 2). 


Theorem 5.9. For positive integers a and b, 


[a5] = 5[a + Va + 4(a/)] 


Proof. Let x = [a, | =art — 


Then bx +1 = ab + — +10 that 


1 . 
x 

x=at = and hence x + = = 2x — a. We now evaluate the area of 
a square with side 2x — a in two different ways (see Figure 5.16). Thus 


x= a)? = a* + 4(a/b) so that x = sla + ja? + 4(a/b)]. 
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alb 


alb 


alb |x 


albx alb 


x albx 


+ x+a/bx = 2x —-a—> 
FIGURE 5.16 


Example 5.7. The golden, Sivel, and bronze ratios. If we seta =b=1 
in Theorem 5.9] we have [1] = =(1 + V5) = q, that is, the golden ratio 
has what may well be the simplest periodic continued fraction: 


In general, the simple periodic continued fractions [n] with period 1 are 
given by 


1 1 
Yn = [nN] =O gO 


(where ~, = @). These are the so-called metallic ratios, and correspond 
to ratios of the dimensions of rectangles analogous to the golden rectan- 
gle. For example, a silver rectangle is one with the property that when 
two squares are cut off (see Figure 5.17(a)), the remaining rectangle (in 
gray) is similar to the original. 

Thus the silver ratio equals - = 2 eS ~ =[2] =~, =1+ 


V2. Similarly, a bronze rectangle is one eit the property that when 
three squares are cut off (see Figure |5.17(b))), the remaining rectan- 


gle (in gray) is ae to the original. Thus the bronze ratio equals 
3x+1 


> =3+- ~ =[3] = = 3 = =(3 + ¥13). The remaining metal- 
lic ratios ~,, are obtained similarly. 


100 5. IRRATIONAL NUMBERS 


2x + 1 3x +1 


(a) (b) 
FIGURE 5.17 


Theorem also yields periodic continued fractions for some 
square roots. For example, for n a positive integer, 


[2n] =n+V¥n2+1=2n+ 


2n + 
ant. 
and subtracting n from both sides yields 
1 —_ 
n?+1=nt+ tT = [1 2n]. 
Zn + 
an-+™. 


Continued fractions for some other square roots are obtained from con- 
tinued fractions with period 2 in the next example. 


Example 5.8. Some periodic continued fractions with period 2. Replac- 
ing a with 2a in Theorem 5.9 yields 


[2a,b] =a+ Ja? +2a/b =2a+ 
b+ 


and subtracting a from both sides yields 


1 
a2 +2a/b=at ——— [a, b, 2a]. 


b+ 


pes 
oS aE, 


Some special cases are 
(a,b) = (n,n): ¥n2 +2 = [n,n,2n], 
(a,b) = (n,2): Vn? +n = [n,2,2n], 
(a,b) = (n,1): Vn? +2 = [n,1,2n], 
(a,b) = (2n, 4): V4n?2 +n = [2n,4,4n], etc. 


5.1 


52 


5.3 


5.4 


5.5 


5.6 


5.7 


5.8 
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5.9. Exercises 


Use the result in Example about integer squares modulo 3 to 
show that V3n — 1 is irrational for n = 1. (Hint: First show that if 
V3n — 1 = a/b for aand b relatively prime, then a? + b* = 3nb?.) 
(a) Show that a rational number to an irrational power may be ir- 
rational. 

(b) Show that a rational number to an irrational power may be ra- 
tional. 

Prove that Vn? — 1 is irrational for n > 3 an integer. [Hint: Modify 
proof 1 of Theorem 5.4] 

Prove that Vn? + 1 is irrational for n > 3 an integer. [Hint: Modify 
the proof of Theorem 5.5] 


Show that (a) p? + (1/7) = 3 and (b) y?—(1/g?) = 4. [Hints: (a) 
See Figure and (b) construct a similar square with side length 


g+1+ Q@.] 
For the metallic ratios g, = [n] in Example show that (a) 


g2 = [1+ n2,1,n?] and (b) 92 = Gp343n = [n3 + 3m]. 
In Example }5.4| we observed that 


1 1 | 
COSonti — 2 2+ {2+ 2 Pciald 


with n nested square roots. Is there a similar expression for 
r TT 

Sin Sa? 

In Example and Exercise 5.7 we encountered finitely many 

nested square roots. What meaning can be attached to infinitely 


many nested roots, such as 


PH ran Pree or [1+ j1+Vit~2 


[Hint: Examine the more general case /a+Va+vVat-: with 
a > 0 first.] 


CHAPTER 6 


Fibonacci and Lucas Numbers 


Learn the particular strength of the Fibonacci series, 
a balanced spiraling outward of shapes, 

whose golden numbers which describe dimensions 
of sea shells, rams’ horns, collections of petals 

and generations of bees. 

A formula to build your house on, 

the proportion most pleasing to the human eye. 


Judith Baumel 
The elements of the sequence 
{Fajnei = {1,1,2,3,5,8, 13,21, 34555, 89) 0} 


defined by F, = F, = land F, = F,_; + F,-2 forn = 3 are known as 
Fibonacci numbers, after the Italian mathematician Leonardo of Pisa (c. 
1170-1240) or Leonardo Fibonacci (for Filius Bonacci, son of Bonacci 
who first studied them (when convenient we set fy) = 0). In 
Thomas Koshy describes the Fibonacci sequence and its twin, the 
Lucas sequence (see Section as “the two shining stars in the vast 
array of integer sequences. They have fascinated both amateurs and 
professional mathematicians for centuries, and they continue to charm 
us with their beauty, their abundant applications, and their ubiquitous 
habit of occurring in totally surprising and unrelated places. They con- 
tinue to be fertile ground for creative amateurs and mathematicians 
alike.” The journal The Fibonacci Quarterly, first published in 1963, is 
devoted to the study of the properties of these sequences. 

On the left below we see a 2013 postage stamp from the Principal- 
ity of Liechtenstein. Across the bottom of the stamp are the first 19 Fi- 
bonacci numbers and the first three digits of F,, = 6765 (the 100 at the 
end is the stamp’s denomination). On the right is a 1999 stamp from 
the Caribbean island nation of Dominica with a likeness of Leonardo. 
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Leonardo Fibonacci and his sequence 


6.1. The Fibonacci sequence in art and nature 


In this section we illustrate with two examples the use of the se- 
quence in art and craft, and one example of its occurrence in nature. 


A Fibonacci chimney 


In the photographs below we see the Fibonacci sequence in neon on a 
chimney at the Turku Energia plant in Turku, Finland. It is a creation of 
the Italian artist Mario Merz (1925-2003) who has a similar installation 
on a chimney at the Center for International Light Art (formerly the Lin- 
denbrauerei) in Unna, Germany, and on the dome of the Mole Antonelliana 
(now the Museo Nazionale del Cinema) in Turin, Italy. 


a ee ae (b) 
The Fibonacci chimney in Turku, Finland 


There are a number of occurrences (or alleged occurrences) of the 
Fibonacci numbers in nature. The one we present in the following ex- 
ample concerns the ancestry of male honeybees. 


Example 6.1. The family tree of a male honeybee. Male honeybees are 
born from unfertilized eggs (parthenogenesis), while female honeybees 
come from fertilized eggs. So males have mother but no father, while the 
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females have both a mother and father. In Figure [6.1] we see the family 
tree of a male honeybee for six generations. The generations are num- 
bered as follows: 1 is the male bee, 2 is its mother, 3 is its grandparents, 
etc. The numbers along the right edge of the tree are the number of bee 
ancestors in each generation. 


of ot ie 08 
ee ee a a 
es i Sepa 

‘ an a ae 

of | Q ---2222222-- 2 


? Female + precio sees encase eeesia 1 
o Male Cie eee sheaee anes 1 
FIGURE 6.1. A portion of the family tree of a male honeybee 


Let Q,,, Dn, and T,, denote, respectively, the number of female bees 
(queens), male bees (drones), and the total number of bees in genera- 
tion n in the ancestry of the male in generation 1, where we have Q, = 0 
and D, = 1. Then T, = Qy + Dy, Dyn = Qn_1 (since only queens have 
fathers), and Q, = Qn—-1 + Dn_1 (since both queens and drones have 
mothers) forn = 1. Hence Q, = Qn_-1 + Qn_2 with Q, = 0 and Q, = 1, 
thus Q, = F,-1- Then Ty, = Qn + Qn-1 = Qn41 = Fn, So that the number 
of bee ancestors in generation 7 is F). 


The Fibonacci sequence in textile design 


The Fibonacci sequence appears to be popular with designers of textiles 
such as tartans, rugs, and quilts. In the tartan design below on the left, the 
widths of the stripes are elements of the sequence. Notice that the rectan- 
gles bordering the largest light gray square approximate the golden rec- 
tangle. In the center is a rug design with two copies of the sequence, one 
ascending, one descending. The quilt design on the right uses the same 
idea in both the vertical and horizontal directions. 


A Fibonacci tartan, rug, and quilt. 
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6.2. Fibonacci parallelograms, triangles, and trapezoids 


There exist many lovely identities relating the Fibonacci numbers. 
The first three we illustrate concern sums of Fibonacci numbers, specif- 
ically the sum of the first n, the sum of the first n with odd subscripts, 
and the sum of the first n with even subscripts: 


n 
(6.1) aR tht th = st 
k=1 
n 
(6.2) Poy-y = Fy + Fg to + Fon-a = Fan 
k=1 
n 
(6.3) >, Fie = Ft Fito + Ban = Finer — 1b 
k=1 


We illustrate these formulas with parallelograms, triangles, and trape- 
zoids whose side lengths are Fibonacci numbers, as shown in Figurel.J 
[Walser, 2001]. Note that the sides of the parallelogram are consecutive 
Fibonacci numbers, and the top, the slanted sides, and the base of the 
trapezoid are three consecutive Fibonacci numbers. 


FIGURE 6.2 


In Figure 6.3|we see an equilateral triangle subdivided into anumber 
of small A and V triangles with side lengths 1. The lengths of the left side 
and the base are the same, hence 1 + F, + F3 + ++: + F, = Fy42, which is 


equivalent to (6.1)). 

In Figure |6.4| we recolor and relabel the previous figure, and gather 
together the light gray pieces into an equilateral triangle, and similarly 
for the dark gray pieces. In the light gray equilateral triangle, the lengths 
of the left side and the base are the same, hence F, +f, +°*:+Foyn_1 = Fon, 
which is (6.2). Similarly, in the dark gray equilateral triangle, the lengths 
of the left side and the base are the same, hence 1+h, +k, +--+ Foy = 


Fon+1) Which is equivalent to (6.3). 
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VV = 
WAV iF ANA, JN / 
7.7, HO 

\/\/ K AY, / 
ce KK KK 


FIGURE 6.4 


6.3. Fibonacci rectangles and squares 


In this section we illustrate several identities involving products of 
Fibonacci numbers using areas of squares and rectangles. We apply the 
identities to show that consecutive Fibonacci numbers are relatively 
prime, construct Pythagorean triples, introduce the Fibonacci spiral, 
and establish some congruence results. 


Example 6.2. An interesting Fibonacci identity results from comparing 
the square of one Fibonacci number to the product of its neighbors in 
the sequence, i.e., 2-5-3? = +1,3-8—5? = —1,5-13—87 = +1,andso 
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on. This leads to an identity named for the Italian mathematician and 
astronomer Giovanni Domenico Cassini (1625-1712) who discovered 
itin 1680. 


Theorem 6.1. Cassini’s identity. For alln > 2, Fy-1Fy44 —Fe = (-1)". 


Proof. Evaluating the area of the region in Figure (6.5(a)jin two different 
ways yields 
Fn-1F +1 + Fn-2hn = Fy ag as 

which is equivalent to F,_1Fy41 — BP = —(Fy-2, — E2_,). Thus the 
terms in the sequence (Fu—1F 41 - Es have the same magnitude 
and alternate in sign. So we need only evaluate the base case n = 2: 
FF, — F? = 150 F,_1Fy41 — FE? is +1 when nis even and —1 when n is 
odd, ie., Fra Fra, — F2 = (-1)”. 


FIGURE 6.5 


With Cassini’s identity we can easily prove the following important 
property of the Fibonacci numbers. 


Corollary 6.2. Pairs of consecutive Fibonacci numbers are relatively 
prime, i.e. gcd (Fy, Fn41) = 1 forn = 1. 


Proof. Suppose a prime p divides both F, and F,,,. Then by Cassini’s 
identity, p divides (—1)", which is impossible. Hence gcd (F,, Fn41) = 
1. 


Example 6.3. A result similar to Cassini’s identity holds for four con- 
secutive Fibonacci numbers when we compare the product of the two 
outer numbers to that of the two inner numbers, i.e., 2:-8—3-5 = +1, 
3-13-—5-8=—-1,5-21—8-13 = +1, etc. 


Theorem 6.3. For all n > 3, Fx-2Fr41 — Fa-1F, = (-1)" 
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Proof. Computing the area of the region in Figure in two ways 
yields F,_2F i414 + Fy_-3F) = Fn_-1 Fy + Frn-2Fn-1, equivalent to F,_»F,44— 
F,_4F, = —(F,-3F, — Fy-2F,-1). As in the proof of Theorem the 
terms in the sequence {F,_2F,41 — je ae have the same magni- 
tude and alternate in sign. Evaluating the base case n = 3 yields F, F, — 
F,F; = 1,80 F,_2Fy41 — F,_-1F, is +1 when n is odd and —1 when n is 
even. 


Corollary 6.4. Pairs of Fibonacci numbers whose subscripts differ by 2 
are relatively prime, i.e., gcd (F,_1, Fy41) = 1 forn = 2. 


Proof. Replacing n by n + 1 in Theorem 6.3] yields Fens — Png = 
(—1)" for n = 2. Ifa prime p divides both F,_, and F,,,,, then p divides 
(—1)", which is impossible. Hence gcd (F,_-1, Fy41) = 1. 


In the next example we illustrate a variety of identities for the square 
of a Fibonacci number. 


Example 6.4. Identities for F2,,. In Figure 6.6| we illustrate two differ- 
ent ways to partition a square with area F?,,. 


(a) (b) 


FIGURE 6.6 


Figure Ollerton, 2008] illustrates the following five identities 


for the square of a Fibonacci number: 
Fist = 20 algae Ea 
SOF ea Pd lg 
(6.4) =F2+F2,+2E,F-1 
= FaaaFy—1 + FrFy—1 + Fr 
(6.5) = Pai, bg a 


110 6. FIBONACCI AND LUCAS NUMBERS 
Similarly, Figure illustrates the following four 
identities: 
Frat = 4FyPy—1 + Fa-2 
= 2F2 + 2F2.,—F2, 
= 41 + 4Fy-1Fa—2 + Fr_2 
= 4F? —4F,_,F,-2 — 3F7_>. 


Similar figures produce many other identities. 


Theorem 6.5. For all positive integers n and k, 


(6.6) FroaFcta + Fake = Fr+kt1: 
Proof. Evaluating the area of the region in Figure (6.7|in two ways yields 
(6.7) Freaheta + FaFe = PrFet2 + Fr-1Pe41- 
Fx 
Frx2 
Fr] 
Fy, Fr} 
FIGURE 6.7 


Observe that in the right-hand side of (6.7), n has been replaced by 
n — 1 and k has been replaced by k + 1, as compared to the left-hand 
side. Iterating (6.7) until n decreases to 2 yields 


FusaPesa + Pe = FaFet2 t+ Fn-1Pe41 
= Fy_-1Fie43 + Fy-2Fe+2 
= Fy_-2Fe44 + Fy-3Fe+3 


= FoFnsk + Fy Fasr—1 = Fasksy 


where in the final step we use F, = F, = 1 and Fyax + Fyan—-1 = Pyar: 


Corollary 6.6. Each Fibonacci number is either the sum or the difference 

of squares of two smaller Fibonacci numbers, i.e., for n positive, 2... + 
= 2 Oe. ee 

Fi = Fony1 and Fry1 — Fra = Fan 
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Proof. Setting k = nin yields F2,, + EB? = Fon41. Identity 
derived from Figure 6.6alis equivalent to F2,, —F2_, = Fra1F, + FyF-1, 
and setting k = n—1in yields F,44F, + F,xFn-1 = Fon, So that 
First — ae = Fon. 


The relationship F2,, + FE? = Fyn41 for the Fibonacci sequence mir- 
rors the relationship 77 + T7_, = T,2 for the sequence of triangular 
numbers encountered in Example |1.5; The sum of the squares of two 
consecutive members of the sequence is another member of the se- 
quence. 


Corollary 6.7. [fn divides m, then F,, divides Fy. 


Proof. Settingk =n—1in yields Fi44F, + Fy Fn-1 = Fon, So that F, 
divides F,,. Setting k = 2n — 1in yields Fn41/on + FrFon-1 = Fan 
so that F, divides F;,,. Continuing in this fashion shows that F, divides 
Fyn for any positive integer k. 


As a consequence of Corollary F; = 2 implies that every third 
Fibonacci number is even (i.e., F3, = 0 (mod 2) fork = 1); Fy = 3 
implies that every fourth Fibonacci number is a multiple of 3 (i.e., Fy, = 
0 (mod 3) fork = 1); & = 5 implies that every fifth Fibonacci number 
is a multiple of 5 (i.e., i, = 0 (mod 5) fork = 1); F, = 8 implies that 
every sixth Fibonacci number is a multiple of 8 (i.e., f, = 0 (mod 8) 
for k = 1); and so on. 


Example 6.5. Fibonacci Pythagorean triples. Since F, and F,., fork = 2 
are relatively prime, we can use them to generate PTs (which will be 
primitive when F, and F,,, have opposite parity). Setting (m,n) = 
(Fie Fe41) for k = 2 in Euclid’s formula yields a = F2,, —F2 = Fy_1F+2, 
b = 2FyFyay, and c = F2 + F2,, = Fox. Thus every Fibonacci num- 
ber with an odd subscript, beginning with F, = 5, is the hypotenuse of 
some PT, and the area of the PT is k,_1 Fy Fy4iF the product of four 
consecutive Fibonacci numbers. 

Similarly, F,_, and F,,, are relatively prime, so setting (m,n) = 
(Fe—1, Fe41) in Euclid’s formula yields a = F2,, — FL, = Fox, b = 
2Fy—-1Fk41, and c = F2_, + F2,,. Thus every Fibonacci number with 
an even subscript, beginning with F, = 3, is the leg of some PT. These 
results are not surprising, since every number greater than 2 appears 
in some PT. 


K+2? 


Example 6.6. Sums of Fibonacci squares and the Fibonacci spiral. A lit- 
tle experimentation leads to the hypothesis that the sum of the first n 
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Fibonacci squares is a product of consecutive Fibonacci numbers: 
174+17=2=1-2, 
124+17+4+2%=6=2:-3, 
17 4+17+ 27+ 37 =15=3:-5, 
174+ 12+27+ 37452 =40=5:°8, 
174+ 12+27+37+45? + 8? = 104 = 8 - 13, etc. 
Theorem 6.8. F? + F? + +--+ FB? = F, Fiat. 


Proof. See Figure 6.8a Brousseau, 1972] . 


(a) 


FIGURE 6.8 


If we rearrange some of the squares in Figure and draw a quar- 
ter circle in each one, we create a portion of the beautiful Fibonacci spi- 
ral, as illustrated in Figure This spiral appears on the 1987 Swiss 
postage stamp, a 10 litas Lithuanian gold coin minted in 2007, and in 
a floor tiling in the Julian Science and Mathematics Center at DePauw 
University in Greencastle, Indiana, as illustrated in Figure 


FIGURE 6.9 
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Example 6.7. Sums of Fibonacci rectangles. If we replace the Fibonacci 
squares in the preceding example with rectangles, we obtain two more 
identities. 


Theorem 6.9. 
(a) FLFR + RFs t+ + Fon-2Fon-1 = Fin-1 — 1, 
(b) FR +R Fy +++ Fon-1Fin = Ein 
See Figure 6.10] where we illustrate the identities for n = 4. 


(a) (b) 
FIGURE 6.10 


6.4. Diagonal sums in Pascal’s triangle 


In Figure 6.11] we see a portion of Pascal’s triangle, presented left- 
justified so that the binomial coefficient (") appears in the nth row and 
kth column for both n and k in {0, 1, 2,3, ...}. The gray arrows indicate 
sums of the elements on the rising diagonals, which appear to yield 
Fibonacci numbers. We also indicate the use of Pascal’s formula (also 
known as Pascal’s rule or Pascal's identity) by the boxed L-shaped region 
for adding two consecutive elements in one row to obtain an element in 
the next row, e.g., 15 + 6 = 21. 

Let D,, denote the sum of the elements in the nth diagonal forn = 0. 


Then Dy = (2) =1,D,= (*) = 1,and 
[n/2] 


= ()*(5 Jena (9) 


k=0 


114 6. FIBONACCI AND LUCAS NUMBERS 


13 
x {2 
“Ce. 


1040 5 1 
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FIGURE 6.11 


forn => 2. If wecan show that D, + Dns = Dniz for n = 0, then we can 
conclude that D, = Fy+;1 forn = 0. 

In Figure we use Pascal’s formula to illustrate D, + Dns, = 
Dn+2 for n odd, e.g., for n = 5; and in Figure for n even, e.g., for 
n = 6. We have used the fact that @) is defined to be 0 whenever k < 0 
or k > n,as indicated by the gray 0’s. 


15 20 15 6 1 
7 2) 33 ap: 21 7 1 


1 8 28 56 70 56 28 8 1 
(a) 


FIGURE 6.12 
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6.5. Lucas numbers 


Asequence closely related to the Fibonacci sequence is the sequence 
of Lucas numbers, 


{Ln} = (13,4, 7,11, 18, 29, 47, 76, 123, 199, ...} 


defined by L, = 1,L, = 3,andL, = Ly_1+Ly_2 forn = 3. When conve- 
nient we set Ly = 2. This sequence is named for the French mathemati- 
cian Francois Edouard Anatole Lucas (1842-1891), who made many 
contributions to the study of the Fibonacci numbers and is responsi- 
ble for giving that sequence its name. Since the Lucas sequence has the 
same recursion formula as the Fibonacci sequence, many of its identi- 
ties are similar. Here are some examples for the Lucas numbers, and for 
both Lucas and Fibonacci numbers. 


Example 6.8. The Lucas analogue of is 
n 
(6.8) Ly = Ly t+ Lg +0 + Ly = Lni2 — 3. 
earl 
The visual proof using Fibonacci trapezoids, as in Figure is some- 


what cumbersome for this identity. But converting the sum into a tele- 
scoping sum succeeds nicely: 


n n 
>, be = >. Cita — Lice) = Ena — La = nse — 3. 
k=1 k=1 


Example 6.9. Cassini's identity for the Lucas numbers is given by 
Ln-tbn41 — L2, = 5(-1)"** for n > 2. The proof of this identity is 
similar to the proof of Theorem 6.1. Relabeling the image in Ce ae 
with Lucas numbers yields Ly_1Ln41 + Ln—2ln = L’, + 14,1, which is 
equivalent to Ly_yLn41 — L4, = —(Ln-2ln — L4_1). Thus the terms in 
the sequence are _ ie have the same magnitude and alter- 
nate in sign. So we need only evaluate the base case n = 2: L,L; —L4 = 
1-4 —3% =—5 so Ly_1Ln41 — L’, is-5 when n is even and +5 whenn 
is odd, ie, Ln-aLn41 — L2, = 5(-1)""’. 


In Exercise 6.4 you can prove the Lucas analogues of identities 
and (6.3); and in Exercise 6.6 you can derive the Lucas analogue of the 
identity in Theorem 

There are many lovely identities involving both the Fibonacci num- 
bers and the Lucas numbers. Here is one example. 
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Theorem 6.10. For all positive integers n and k, 
Fettbnti + Fela = Ln+kts- 


Proof. Evaluating the area of the region in Figure in two ways 
yields 


(6.9) Fytibner + Felbn = Felbns2 + Pe-sbn4s- 
Frey Fn_2 
Fret F, 
Fy Fj 
Ln l Ly Ly Ly 1 
(a) (b) 
FIGURE 6.13 


Observe that in the right-hand side of (6.9), k has been replaced by 
k — 1 and n has been replaced by n + 1, as compared to the left-hand 
side. Iterating until k decreases to 2 yields 


Feptbnss + Pebn = Fyebnse + heals 
= Fe-qlns3 + Pe-2bn+2 
= Fy_-glinyg + Fe-zhnss3 


= Folysk + ylnte-1 = Lntkty 
where in the final step we use Fy = Fy = Land Lys, tl nspe1 = Iyer: 


The image in Figure a relabeling of Figure can be used 
to prove the following theorem, a Cassini-like identity for products of 


Fibonacci and Lucas numbers. 
Theorem 6.11. For all n > 2, Fy-1Ln41 — Fyn = (-1)". 
Compare this identity to those in Theorem iw ae = 
(-1)"] and Example 6.9 [Ln—1Ln41 — 12 = 5(-1)""’]. 
Proof. Evaluating the area of the region in Figure in two different 
ways yields 
Fr—tensi + Fa-2bn = Faln + Pr-ibn-1, 
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which is equivalent to F,_4ln44—Faln = — Fy-2lbn — Fy—-1Ln_1). Thus 
the terms in the sequence {F,,_Lyj44 — 18 oa have the same magni- 
tude and alternate in sign. So we need only evaluate the base casen = 2: 
F,L3 — Lz = 180 Fy_yLy41 — FyLn is +1 when n is even and —1 when 
nis odd, i.e., Fy-aLn41 — Faly = (-1)". 


With the identity in Theorem we can easily prove the follow- 
ing property of the Lucas numbers, analogous to Corollary for the 
Fibonacci numbers. 


Corollary 6.12. Pairs of consecutive Lucas numbers are relatively prime, 
ie, gcd (Ly, Lnai) = 1forn = 1. 


Proof. Suppose a prime p divides both L, and L,,,. Then by Theo- 
rem p divides (—1)", which is impossible. Hence gcd (Ly, LZn41) = 


6.6. The Pell equations x? — 5y” = +4 and Binet’s formula 


In Corollary 6.6 we established the identity F2,, + F2 = Foy44, and 
the Lucas version L2,,, + 14. = 5F)n+, follows from setting k = n 
into the identity in Exercise 6.6 at the end of this chapter. Combining 
these yields L2,, + L4, = 5F2,, + 5F? or, equivalently, 17, —5F7,, = 
— (12, — 5F?). Thus the terms in the sequence {L%, — SF} have the 
same magnitude and alternate in sign. Again we need only evaluate the 
base case n = 1: L3 — 5K? = —4; hence L2 — 5F? = 4(—-1)" forn > 1. 
This identity was discovered in 1950 by P. Schub 

Thus the pairs (Ly, F,) are solutions to the Pell equations x?—5y? = 
+4. We now solve these Pell equations using the methods presented 
in Chapter 3) to illustrate how the pairs (L,,, F,) arise as solutions, and 
proceed to find explicit formulas for the Fibonacci and Lucas numbers. 

In Section 5.2) we established the identity x? — 2y? = 2(x + y)*— 
(x + 2y)? using Figure This is the d = 2 case for the more general 
identity (d — 1)(x2 — dy?) = d(x + y)* — (x + dy)?, which holds for 
all positive integers x, y, and d. When d = 5 we have 4(x2 — 5y*) = 
Six + y) — (x + 5y)*, which we illustrate in Figure in the form 
(x + Sy)? + 4x? = 20y? + 5(x +y)’. 

The identity 4(x? — 5y2) = 5(x + y)° — (x + 5y)’ is equivalent to 


x+y\? x +5y\ 
2 ea? ae 
ae 5( 2 ( 2 } 


This establishes the following theorem, an analogue of Theorem 5.21 
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ILI 
LI 


Theorem 6.13. Let x and y be positive integers. Then 


FIGURE 6.14 


x+5y\ x+y)? = — 
(6.10) (~<*) -5(=5*) = £4 ifand only ifx? — Sy? = 4. 


Note that solutions to x? — 5y? = +4 must have x and y both even 
or both odd, hence (x + y)/2 and (x + 5y)/2 are integers. Set 


1 1 
(6.11) Xn4+1 = 5 On + 5Vn) and Yn+1 = zn + Xn) 


in (6.10) yielding x2,, — 5y2,, = +4 ifand only if x2 — 5y? = ¥4. Set- 

ting (x1, ¥,) = (1,1) yields x? — 5y? = —4, hence x2 — 5y? = 4(-1)". 

In Table el we have the solutions (xp, ¥,) to x2 — 5y2 = 4(—1)” gen- 
(6.1 1)) 


erated by for n from 1 to 9. 
TABLE 6.1 
nj1 234 5 6 7 8 9 
X,\/1 3 4 7 11 18 29 47 76 
¥,|1 12.3 5 8B 13 21 34 


It appears that x, = L, andy, = Fy, but to verify this we must show 
that Xn41 = Xn +Xn—-1 aNd Yp41 = Vn + Yn—1 for n = 2. From (6.11) we 
have 


1 
Xn4+1 = 5 On + 5Yn) 


1 5/1 1 
Sa ain-1 a pyn-1 


2 2 
1 5 1 
= z*n + Ties + 4 (2xn tna) 


= Xn + Xn-1 


6.6. THE PELL EQUATIONS x? — 5y? = +4 AND BINET’S FORMULA 119 

and 

1 

Yn+1 = 5 On + Yn) 
1 1/5 1 
= 5Yn - 3 e= as stn-1] 

1 5 1 
= 59n + G¥n-1 + F 2m — Yn-1) 
= Yn + Yn-1: 


Since x; = 1,x2 = 3,and y, = y2 = 1, we have x, = L, and y, = F, for 
n> 1,and (6.11) yields the identities 


1 1 
(6.12) Lmta = 5(Ln + SF) and Fre1 = 5 (Fa + Ln): 


Example 6.10. A square root test for Fibonacci-ness. Is 25,075 a Fi- 
bonacci number? Is 75,025 a Fibonacci number? If y is a Fibonacci 
number, then by the identity L2 — 52 = 4(—1)", either ./5y2 — 4 or 
J 5y? + 4 is an integer (in fact, a Lucas number). When y = 25,075, 


neither ./5(25, 075)? — 4nor /5(25, 075)? + 4 is an integer, so 25,075 
is not a Fibonacci number. But when y = 75,025, /5(75,025)? —4 = 
167,761, so 75,025 is a Fibonacci number (it is Fy; and 167,761 = 


Los). 


Example 6.11. Binet’s formula. We now find explicit expressions for 
F, and L,,. Evaluating the generating functions Gf(t) = yee Ft” and 
G,(t) = ae L,t” for the sequences {F,,} and {L,,} using Fys2 = Fya4 + 
F, and Lny2 = Lng, + Ly forn = 0 yields Gp(t) = t/( —t — t?) and 
G,(t) = (2 —t)/(1 —t — t?), from which we obtain 


_adfrtvs\ (1-v5\" 
‘n= Te 2 ae 


_ (145  fiaasy 
ln =| hs 


using partial fractions and geometric series. The formula for F, is 
known as Binet’s formula, after the French mathematician Jacques 
Philippe Marie Binet (1786-1856) who discovered the formula in 1843. 
However, it was known a century earlier to another French mathema- 
tician, Abraham de Moivre (1667-1754). 


and 
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. 1tv5. ; 1-5 1 
The fraction ws is, of course, the golden ratio g, and ms =—-_, 


? 
so we have simple expressions for F,, and L,, in terms of the golden ratio: 


1 —n _ on yon 
F.= Tele —(-g) "| andL, = g" + (-9) 


; at 
or, since l-—@= er 


1 n = n = n 
i= ele —(1-9)"| andL, = 9" + 1-9)”. 


Example 6.12. Rational approximations to V5. Since Fibonacci and Lu- 
cas numbers are solutions to the Pell equations x? — 5y? = +4, their 
ratios approximate V5. Dividing L2 — 5F2 = 4(—1)” by F? yields 


2 
Ln n 4 
2) 25 S4) 
(] ( ) Rp 


Thus lim),-5co Ln/F, = V5, and L,,/F, approximates V5 (larger than V5 
when nis even, smaller when n is odd) with an absolute error less than 
4/3F2. This error bound follows from 


4 (Ly 
—=(—+Vv5 
ao (et) 


since = + V5 21+ V5 >3. 


Example 6.13. Limits of F,41/F, and Lyj+41/Ly. There are various ways 
to examine the behavior of F,4,/F, and Lyj41,/Ly as n increases with- 
out bound. The values of the ratio F,,,,/F, for various values of n can 
be seen in Figures 6.8} and where the rectangles with dimensions 
F,41 X F, appear to be approximately “golden” (see Figures and 
5.11). 

Asimple way to evaluate lim,,_,,. F,41/F, utilizes the second identity 


in (6.12) and Example 6.12} 


F, 1 i 1 
lim = = im 5 (14) = (1+ v5) = 9 
n 


noo fF, noo 2 i 2 
and similarly limy50. Ln41/Ln = from the first identity in (6.12). 


For a comprehensive overview of Fibonacci and Lucas numbers and 


their properties, see [[Koshy, 2001]. 
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6.7. Exercises 


Use illustrations similar to Figures and to show that 

FP, + RF, +++ FinFonge = + FE + Fe to + Fins 

= Fonti1Pans2 — 1. 

[Hint: Partition an Fyn41 X Fon42 rectangle (with one square miss- 
ing) in two different ways.] 
In Example(6.4}we partitioned a square to illustrate (6.4), the iden- 
tity F2,, = B? + F?_, + 2F,F,-1. A cubic version is F?,, = F? + 
F3_, + 3Fy416,Fy-1. Can you partition a cube to illustrate this? 
Show that the Lucas sequence can be obtained from the Fibonacci 


sequence and vice-versa: F,_1 + F,41 = Lyand Ly_1+Ly41 = 5F. 
Prove the Lucas analogues of identities (6.2) and (6.3), i-e., 
n 
(a) Liat boR-1 = Lan — 2 and 
(b) Mins Lox = Lenss — 1. 


Show that (a) FL, = Fj, and (b) FL, + BL, = 2Fysx. 

Prove the Lucas analogue of Theorem 6.5: Lys ly, + Lyle = 
OF +k+1: 

Prove the Lucas analogue of Theorem 6.8; 12 + 13 +++ 12 = 


Lylgag = 2. 

The numbers y” and (—¢) ” appear in formulas for F, and L,, in 
Example 6.11} When n = 1 these numbers are the roots of the 
quadratic equation x? — x — 1 = 0. Find a quadratic equation with 
roots gp” and (—~) ". 

We showed that F;, = 0 (mod 2) fork = 1 as a consequence of 
Corollary (6.71 Now show that L3, = 0 (mod 2) fork = 1. 

Show that (a) the arithmetic mean of F, and L,, is F,4,, and (b) the 
harmonic mean of F, and Ly, is Fon /Fy+1- 


CHAPTER 7 


Perfect Numbers 


Perfect numbers, like perfect men, are very rare. 


René Descartes 


Man ever seeks perfection but inevitably it eludes him. He has 
sought perfect numbers through the ages and has found a very few. 


Albert H. Beiler 


A perfect number is a positive integer n that is equal to the sum of 
its proper divisors (those divisors less than n). For example, 6 is perfect 
since 6 = 1+ 2 + 3, and 28 is perfect since 28 =1+2+4+7+414. 
Perfect numbers have been studied since the time of the early Greek 
mathematicians. Perfect numbers appear in Definition 22 in Book VII 
of Euclid’s Elements [Heath, 1956] as “A perfect number is that which is 
equal to its own parts.” The Greeks were aware of the first four perfect 
numbers, 6, 28, 496, and 8128, and in the Middle Ages mathematicians 
found the next three: 33,550,336, 8,589,869,056, and 137,438,691,328. 
Two of the oldest unanswered questions in mathematics are (1) do in- 
finitely many perfect numbers exist, and (2) are there any odd perfect 
numbers? 


7.1. Euclid’s formula 


In addition to the definition mentioned above, Euclid provides a for- 
mula for perfect numbers in Proposition 36 in Book IX of the Elements: 
“Ifas many numbers as we please beginning from a unit be set out in dou- 
ble proportion, until the sum of all becomes prime, and if the sum is mul- 
tiplied into the last make some number, the product will be perfect.” The 
expression “double proportion” refers to a geometric progression with 
common ratio 2. For example, 1+2+4 = 7and7 is prime, so 4-7 = 28is 
perfect. The numbers constructed by this formula are necessarily even. 

Since 1+2+4+-+2*-1 = 2* — 1 (a formula known to the early 
Greeks), Euclid’s formula states that if 2* —1 is prime, then 2*~1(2* —1) 
is perfect. But 2" — 1 prime implies that k is prime, so we have 


123 
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Theorem 7.1. Euclid’s formula for even perfect numbers. [fp and q = 
2? — 1are prime, then N, = 2?~1q is perfect. 


Proof. The proper divisors of N, are 1, 2,4,...,2?~1, q, 2q,4q, ...,2?-*q. 
Figure Goldberg] shows that rectangles with areas equal to the 
proper divisors sum to N, = 201g, 


<—_—___- q= i 
FIGURE 7.1 


Nearly 2000 years after Euclid, the Swiss mathematician Leonard 
Euler (1807-1783) proved that every even perfect number has the form 
expressed by Euclid. 


0.2 2.222 2.2.8 22.8.8 2 O20 m. 


Leonhard Euler, Swiss postage stamp, 2007 


Mersenne primes 


Mersenne primes, named for the French mathematician and theologian 
Marin Mersenne (1588-1648), are primes of the form q = 2? — 1. These 
are precisely the primes q used to construct even perfect numbers in 
Theorem Thus the search for even perfect numbers is the search for 
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Mersenne primes. As of 2018, 50 Mersenne primes are known, the largest 
277,232,917 _ 4 has 23,249,425 digits. The 2004 Liechtenstein postage 
stamp shown below celebrates the 39th Mersenne prime, discovered in 
2001. The 15 largest known Mersenne primes were discovered by the 
Great Internet Mersenne Prime Search (GIMPS) at www.mersenne.org. 


Marin Mersenne and the 39th Mersenne prime 


7.2. Even perfect numbers and geometric progressions 


Theorem 7.2. The even perfect number N, is the sum of p terms of a 
geometric progression with first term 2°~1 and common ratio 2, i.e., Ny = 
Cae cet ah eee 

For example, when p = 5 we have Ng = 496 = 16+32+64+128+ 
256. 


Proof. We slice the rectangle representing N,, in Figure 7.1 vertically, as 
shown in Figure to form p rectangles whose areas are the terms of 
the geometric progression. 


FIGURE 7.2 
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The algebraic proof is equally simple: 


Np = 2P-1 (2? — 1) = 2°71 (142444 ep DPT), 


7.3. Even perfect numbers and triangular numbers 
Each even perfect number is also the sum of an arithmetic progres- 
sion; e.g., 
N2=1+2+43, 
N3=14+2+34+4+4+5+6+/7, 
Ng =142434+*=+30 431, etc 
Each sum is a triangular number of the form N, = Tzp_;, as shown 
in Figure (for pp = 5, Ns; = 16-31). Here we represent N, as a 


(P< (2 — 1) rectangular collection of dots, and rearrange the dots 
as indicated. 


op-| 


Ny = 2?-1(2? — 1) 


FIGURE 7.3 


The algebraic proof is equally simple: 
ue Mee ed). 


Np = 2-4 (2? — 1) = ——]—— = Tp. 


Hence we have 


Theorem 7.3. Every even perfect number is a triangular number. 
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Even perfect numbers and Pascal’s triangle modulo 


Pascal’s triangle modulo 2 is easily constructed using addition mod 2: 
0+0=0,1+1 = 0,and1+0 = 1. Here are the first 17 rows (the 
rows are numbered 0 to 16 counting from the top): 


10001 
1.1 Gee 1.1 
1010101 
11iiaiii1i1i1 
100000001 
1.1 GROROROROeG 1 1 
10100000101 
1111000031111 
1000100010001 
110031%10031%10012~1 
10101010103101041 
Tit idr~ii~iiriiriitirtdrzd1diit 
1 GROROROROSOSCROROTOROROSORORO 1 


Central inverted triangular blocks of zeros (shaded gray) appear with the 
base in row n = 2*, and the number of zeros in the gray triangle is T>x_1, 
an even perfect number when k is prime. In the figure we see the perfect 
numbers T>2_, = 6 and T>3_, = 28, and the top row of Tp4_, (which is not 
perfect). A row with exactly two 1’s interpreted in binary (i-e., 11 equals 3, 
101 equals 5, 10001 equals 17, etc.) are the Fermat numbers (see Exercise 
2.8). 


In Chapter |1} we showed that every hexagonal number Hy is a trian- 
gular number by illustrating H, = Tz,_, = n(2n— 1). Hence each even 
perfect number is also a hexagonal number: 


Ny = o0-* (2? = 1) = Top_4 = Hop-1. 


In the next four sections we explore some of the many consequences 
of ae oe rep cee of N, as a triangular number. The results in Sec- 
tions|7.4 7.4, |7.5} and|7.6|and Exercise 7.4 are motivated by the data in Table 
where we evaluate Ny modulo 7, 9, and 12 for the first seven even 
perfect numbers. 
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TABLE 7.1 
p Ny | Ny (mod 7) | Ny (mod 9) | N, (mod 12) 
Z 6 -1 6 6 
3 28 0 1 4 
5 496 -1 1 4 
7 8128 1 1 4 
13 33550336 1 1 4 
17 8589869056 -1 1 4 
19 | 137438691328 1 1 4 


7.4. Even perfect numbers modulo 9 


Here is an interesting observation about even perfect numbers: If 
you sum the digits of any even perfect number (except 6), then sum the 
digits of the resulting number, and repeat this process until you get a 
single digit, that digit will be 1. For example, for the sixth even perfect 
number 8,589,869,056 we have 


8,589, 869, 056 > 64 > 10 > 1, 
and for the seventh even perfect number 137,438,691,328 we have: 
137, 438,691,328 > 55-10-11. 


Summing digits of a number corresponds to evaluating the number 
modulo 9 since each power of 10 is congruent to 1 modulo 9. See Ta- 
ble {7.1| for some values of N, (mod 9), which leads to the following the- 
orem. 


Theorem 7.4. Every even perfect number greater than 6 is congruent to 
1 (mod 9). 


Proof. First note that when p is an odd prime, 2? — 1 = (—1)? +2 = 
1 (mod 3), so that N, = Tpp_1 = T3n41 for some n. See Figure 
for a visual proof that T3,,,; = 1 (mod 9). We also saw this result in 
Theorem 


7.5. Even perfect numbers end in 6 or 28 


Take another look at the seven even perfect numbers in Table 
each one ends in 6 or 8, and those that end in 8 actually end in 28. We 
now show that this observation holds for every even perfect number. To 
do so, we prove the following lemma concerning 75,41 and T5n+2.- 


Lemma 7.5. 75,4, = 1 (mod 5) and Tsy42 = 3 (mod 25). 
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FIGURE 7.4. Ny = T3n41 = 1+ 9Tp 


Proof. In Figure (7.5(a)| we show that T5,4; = 1 (mod 5) by partition- 
ing a triangular array of dots representing 7;,,,; into 20 copies of T,, (in 
two shades of gray), five copies of T,_, (white dots), and a single black 
dot. In Figure we show that 75n42 = 3 (mod 25) by partition- 
ing a triangular array of dots representing 75,42 into 25 copies of T, 
(in two shades of gray) and three black dots (we saw these results in 
Theorem 2.3). 


(a) Teng = 20T, + 5T,-1 +1 (b) Tne2 = 257, +3 


FIGURE 7.5 


Theorem 7.6. Every even perfect number ends in 6 or 28. 


Proof. Since N, = 6 we need only consider N, for p odd. When p is odd 
there are two cases: p = 4k + land p = 4k +3. Ifp = 4k + 1, then 
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2?-1=2-16"-1=1 (mod 5), ie. 2? —1 = 5n+1 for some positive 
integer n. Hence N, = T2p_1 = Tsn41 = 1 (mod 5), so that in base 10, 
Np ends in 1 or 6. Since N,, is even, it ends in 6. 

If p = 4k + 3, then 2? —1 = 8-16" —1=2 (mod 5), ie, 2? —1= 
5n + 2 for some positive integer n. Hence N, = Tpp_1 = Tsn42 = 3 
(mod 25), so that in base 10, N,, ends in 03, 28, 53, or 78. Since N, is a 
multiple of 4 for p = 3, it ends in 28. 


7.6. Even perfect numbers modulo 7 
With the lone exception of Nz = 28, every even perfect number is 
congruent to 1 or —1 modulo 7 (see Table ap. To prove this result 
(Wall, 1984] , we use the fact that every prime p ¥ 3 is congruent to 1 
or 2 modulo 3. 


Theorem 7.7. Every even perfect number N, with p # 3 is congruent to 
1 or —1 modulo 7, i.e., 


p = 1 (mod 3) = N, = 1 (mod 7), and 
p = 2 (mod 3) > N, = —1 (mod 7). 


Proof. If p = 3k + 1, then 2? —1 = 2-8* —1 = 1 (mod 7), so that 
N3x41 = T7n41- In Figure '7.6(a)| we show that 77,4; = 1 (mod 7) by 
partitioning a triangular array of dots representing 77,4, into 35 copies 
of T,, (in two shades of gray), 14 copies of T,_1 (white dots), anda single 
black dot. Thus N3;,4, = 1 (mod 7). 
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(a) Trnga = 357, + 14T,-1 +1 


FIGURE 7.6 
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Ifp = 3k +2, then 2? —1 = 4-8" —1 = 3 (mod 7), so that N3z42 = 
T7n+3- In Figure [7.6(b)| we see that T7,43 = 6 (mod 7) by partitioning 
a triangular array of dots representing 77,43 into 49 copies of T, (in 
two shades of gray) and six black dots. Hence N3;,42 = 6 (mod 7) or, 
equivalently, N3;,, = —1 (mod 7). 


7.7. Even perfect numbers and sums of odd cubes 


Observing that Nz = 28 = 13 + 33, N, = 496 = 13 + 33 +53 +73, 
and N, = 8128 = 13+ 33 +53 +--+ 15° yields the hypothesis for the 
following theorem. 


Theorem 7.8. Every even perfect number N, with p = 3 is the sum of 
the first n odd cubes forn = 2°-/?, i.e, Ny = 13 +33 +--+ (2n— bie 
Note that the last integer cubed is 2®*1)/2 — 1, 

Proof. When n = 2-1/2 we have 2n” = 2?, so that Ny = Top_1 = 
T,,,2_1. In Figure we show that T>,2_, = 1: 4? 983" ep 
(n= 1) «(2h 1) (for p = 5 and n = 4), representing k? as k copies 
of k? dots. The formula for T2,2_, as a sum of odd cubes can also be 


proved by mathematical induction. Hence N,, (for p odd) is the sum of 
the first n = 2-1/2 odd cubes. 


FIGURE 7.7 


7.8. Odd perfect numbers 


As noted at the beginning of this chapter, whether or not odd perfect 
numbers existis an open problem. Nonetheless, in this section we prove 
the following simple result about odd perfect numbers: Jf N is an odd 
perfect number, then N has at least three distinct prime factors. 
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It is common in number theory to let a(n) denote the sum of all the 
divisors of n, including n itself. Hence n is perfect if only ifo (n) = 2n 
or, equivalently, o(n)/n = 2. 

To show that an odd integer n with just one or two distinct prime 
factors cannot be perfect, we use the following lemma. 


Lemma 7.9. If r is a real number such that0 < r < 1 and k is any 
positive integer, then 


1 
Ltrtretetrk < —. 
L=7 


Proof. See Figure and observe that the area of the darker gray rec- 
tangle is less than the area of the entire rectangle in two shades of gray. 


FIGURE 7.8 


Theorem 7.10. JfN = p" or N = p"q™, where p and q are distinct odd 
primes and m and n are positive integers, then N is not perfect. 


Proof. Let N = p”. Then the divisors of n are 1, p, Dp, ...,p”, and since 
0 <1/p < 1 with p = 3 we have 


o(p") 1 
of ae eee) 


(5) e+) 
=1+—4+(-]} +--+ [- 

Pp Pp Pp 

1 1 3 
< —_. < ——__ = - 
1—(1/p)~ 1-1/3) 2 


so that N = p” is not perfect. The o function, like many number theo- 
retic functions, is multiplicative, that is, if m and n are relatively prime 


<2; 
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positive integers, then o(mn) = o(m)o(n). Now let N = p"q”™. If p 
and q are distinct odd primes, then one of them (say p) is at least 3, and 
the other (q) is at least 5. Hence 


o(q™) 1 Z 1 _ a 
qm ~1-(/q)~1-G/5) 4 
consequently 
nym n m 3.5 15 
COED, 2 OGD 8 ae 
prq™ p" gq 2 4 8 


and thus N = p"q" is not perfect. 
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7.9. Exercises 


Show that every even perfect number greater than six is the sum 
of three distinct triangular numbers, e.g., Nj = 496 = 2314+ 210+ 
55 = To1 + T29 + Tyo. 

Here is another pattern relating even perfect numbers and trian- 
gular numbers: 


No =6=2:°3=3:-14+3=37,+T), 
N3 =28=4-7=3-64+10=373+7, 
Ns = 496 = 16-31 =3-120 4+ 136 = 375 + Ty, 
N7 = 8128 = 64-127 =3- 2016 + 2080 = 3763 + Té64, etc. 
State and prove a theorem about even perfect numbers. 
Here is a third pattern: 

N3 = 28 =1-:10+3-6=7,T, + ToT, 

N, = 496 = 6: 36+ 10: 28 = T3Tg + T,T7, 

N7 = 8128 = 28-136 + 36-120 = T7Ty¢ + TgTjs, etc. 
State and prove a theorem about even perfect numbers greater 
than 6. 


If p is an odd prime, show that N, = 4 (mod 12). [Hint: See The- 
orem [7.41] 

As a consequence of the preceding exercise, in base 12 every even 
perfect number greater than 6 ends in 4. Show that more is true— 
in base 12 every even perfect number greater than 28 ends in 54. 
Express N,, in binary (base 2) notation. 


Show that if Ny is an even perfect number, then ,/1 + 8N, is an 
integer, e.g,V1+8-6 = 7,v¥1+8-28 = 15,v1+8- 496 = 63, 


etc. 
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7.8 Prove that if N, ends in 6 for p = 3, then the digit preceding 6 must 
be odd. 
7.9 Even perfect numbers are constructed from the Mersenne 
primes, primes of the form M, = 2? — 1 where p is also prime. 
The first few are 3, 7, 31, 127, 8191.,.... 
(a) Show that for p=3, M, =1 (mod 6) and M, =7 (mod 24). 
(b) Show that for p = 3, the units digit (in base 10) of M, is 1 or 7. 
7.10 Can a Mersenne prime be the hypotenuse element of a PPT? 


Solutions to the Exercises 


Chapter 


1.1 See Figure§1.1|for illustrations of (a) 42+0, = Tg and (b) 52+0, = 
7: 


(a) (b) coo00 

leoReRen®) OB Hi & eoo0o°0 

oo0o00 0 O & & ocoo0o°o 

oRekene) ooo & eooo0°o 

oo0o0°0 oo0o°0 ecoo0o0o —s 00000 
eoo00°0 - 00000 eo0000 eo00000 
eo0000 000000 eo0000 eo000000 
eoo0°0 eooo0000 eoo0°0 e©0000000 
eo0o0°0 eoooo00o00 e0o0000 eo00o00o000000 

FIGURE S1.1 


1.2 See Figure $1.2] where we illustrate 1 + 9 + 92 = Ty 43432. 


AR 


FIGURE S1.2 


1.3 Yes, there are infinitely many more, since is equivalent to T, = 
2Ty @ T3n+4p+3 = 2T2n+3p4z2- 

1.4 In Figure 61.3] we illustrate (2n + 15° = T3n41 —T, forn = 3. 

1.5 The kth n-gonal number is 


k k 
5 lm —2)k—(n-4)] = 5 [n—2) (k- 1) +2]. 


Assume k > 3 and n= 3. When k is even the number factors as 
= [(n—2) (k—1)+2] with both factors integers greater than 1. 
135 
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*) 

00 
<n-> len-> e000 
$33 - 
os £ 
foe) 


— 2n+1— 


FIGURE S1.3 
When k is odd, k — 1 is even and the number factors as 
k- aa, with both factors integers greater than 1. 
1.6 See Figure $1.41 
© 
Pe © 
(?) 
.*) ee 
©e—°o® = @@ 4 2: 
©- 0-0-0 ©oe 
©-0-0 ©eo0e L020°8 
© e® e® oe” © © 
eo @ @ @ @ @ 
FIGURE S1.4 


1.7 Using the hint, we have 12 + 32 +52 ++ + (2n—1)* = 7, + 
(Tz + T3) + (Ty +75) + + + Ton-2 + Ton-1) = Tetzn-1- 
1.8 The theorem is: Four times the nth pyramidal number is the (2n)th 
tetrahedral number; i.e., 4 Pyr,, = Tetzn, and a proof is 
400 +1)(2n+1)  (n)(2n+1)@n+ 2) 
6 7 6 
You can illustrate this result by restacking the tetrahedral pile of 
cannonballs into four pyramidal piles, as indicated here for the 
eighth tetrahedral number 120: 
120 = (1+ 3) +(6+10) + (15 + 21) + (28 + 36) 
=4+16+ 36+64=4(1+4+9+ 16). 


2. : fe 
Hence Pyr,, is a of a binomial coefficient: 


ee 
yYIn = 4 Cton = z( 3 ) 


1.9 
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The nth octahedral number is the sum of the (n — 1)st and nth 
pyramidal numbers and hence it equals 


(n-—1)n(2n-1) n(n+1)(n+4+1) = n(2n? +1) 
—— es |e Se 


Chapter (2 


Placing a 1 at the end of anumber k in base 9 produces the number 
9k + 1. Ifk = T,, then 9T, + 1 = T3,4, from Theorem 2.2. Since 
1 is triangular, so are 119, 1119, etc. [Also see Exercise 1.2 and its 
solution. | 
The (2k + 1)-stn-gonal number is (2k + 1) +(n—2)To,, but To, = 
k(2k +1), so the (2k +1)-st n-gonal number is a multiple of 2k +1. 
Similarly, the (2k)th n-gonal number is a multiple of k (but not 
necessarily 2k). 
In Theorem for 

(a) set (n, a,b,c) = (3k — 1,2k, 2k, 2k), 

(b) set (n, a,b,c) = (3k, 2k, 2k, 2k), and 

(c) set (n,a,b,c) = (3k +1,2k +1,2k +1,2k +1). 
From the proof of Theorem we have T3n4, = 9T, + 1 and 
T3n-1 = T3n = 0 (mod 3). Hence for every positive integer t, t is 
triangular if and only if 9t +1 is triangular. This yields the sequence 
(a, b) = (9,1), (81, 10), (729, 91),..., Le. (a, b) = (9*, (9% —1)/8) 
fork =1,2,3,.... 
Observe that n? — nis even and p and 2 are relatively prime, hence 
n? — nisa multiple of 2p. 
We need to show that n? = n (mod 6). Set p = 3 in the preceding 
exercise. 
We need to show that n° = n (mod 10). Set p = 5 in Exercise 2.5. 
It suffices to show that 22” ends in 6 for n > 2 and that its tens 
digit is odd. Since the exponent 2” is a multiple of 4 when n = 2, 
we need only show that 2** ends in 6 fork > 1. Setting p = 5 
and n = 2* in Corollary 2.10] yields 2** = (2*)" = 1 (mod 5). So 
(a) 2** ends in 1 or 6, but since it is even, it ends in 6, and (b) since 
2** is divisible by 4, it ends with 16, 36, 56, 76, or 96, i.e., the tens 
digit is odd. 
It suffices to show that 22” = 4 (mod 12) forn = 1. In this case 
the exponent 2” is even, so 2 = ew = 1 (mod 3),so 22” =1, 
4, 7, or 10 (mod 12). But 22” is a multiple of 4, hence 22” = 4 
(mod 12). 
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Yes, since the sum in question is the nth tetrahedral number from 


Section f1.5} 


n 
Te: = ye _ n(n + 1)(n + 2) agg ee 
6 3 
k=1 
Hence a T, = 0 (mod T,,) if and only if me is an integer, i.e., 
n = 1 (mod 3). 
Since P, = T3n_1/3, we require T3,_, = 0 (mod 5). From Theo- 
rem 3n — 1 = O or 4(mod 5), hence 3n = O or 
1 (mod 5) so that n = 0 or 2 (mod 5). 
Each term in the expansion of the given determinant is, except for 
sign, the product of all possible row indices and all possible col- 
umn indices, that is, (100! ); and this is the absolute value of ev- 
ery term. Now 101 isa prime, so by Wilson’s theorem, 100! = —1 
(mod 101). Hence (100! )? = (—1)” = 1 (mod 101), as required. 
Suppose a prime q < _ p divides (p—1)!+1. Since q divides 
(p — 1)! it also divides (p — 1)!+1 — (» — 1)!= 1, which is im- 
possible. But p divides (p — 1)! +1 by Wilson’s theorem, hence it 
is the smallest prime that does so. 


Chapter 


With five sailors we have x = 5a+1,4a = 5b+1,4b = 5c +1, 
4c = 5d +1, 4d = 5e +1, and 4e = 5y + 1. Eliminating a, 
b, c, d, and e yields the Diophantine equation 1024x — 15625y = 
11529. Since (x,y) = (-—4776,—313) is a solution to 
1024x—15625y = 1, (%9, Yo) = (—4776-11529, —313-11529) isa 
solution to 1024x — 15625y = 11529. Hence the general solution 
is (x,y) = (—4776 - 11529 + 15625k, —313 - 11529 + 1024k). 
For k = 3525 we have positive solutions, the smallest of which is 
(x,y) = (15621, 1023). 

Yes. Let ee ab be a set of 1,000,000 distinct primes. Then 
p; and p* are relatively prime when i # j, so by the Chinese 
remainder theorem, there is an integer x such that x = 
—k (mod p7) for 1 < k < 1,000,000. So x + k is divisible by 
pz (ie., x + k has a repeated prime factor). Thus {x + 9 aaa 
is a set of 1,000,000 consecutive integers each of which contains a 
repeated prime factor. 

Use mathematical induction: x,y, = 6 and Xn41Vy41 = 6x2 + 
17xXnYn + 12y7?. 
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If x2 —2y? = +1, then x2y? = (2y2 + 1) i Tpy2. lfx7—277 = 
—1, then x2y2 = x2 [(x2 + 1)/2] = Te 

The problem asks for integers k and m such that 1 + 2 + ++ + 
(k-—1) = (k+1) + (kK+2)+-+-+m. This simplifies to T,_, = 
Tim — Tk, OF Tm = Te-1 + Ty so that Ty, = k?. Since 50 < m < 500, 
we see that k = 204 from Example 3.6 and Table 

Proof 1. Solutions (x, V,) to the Pell equation x? — 2y” = 1 in Ta- 
ble eee the pairs (2y7, x2) of consecutive powerful numbers, 
e.g., the solution (99,70) for n = 3 yields the pair (9800,9801). 
Similarly, solutions (x,,, ¥,) to the Pell equation x* — 2y? = —1in 
Table yield the pairs (x2,2y) of consecutive powerful 
numbers, e.g., the solution (41,29) for n = 2 yields the pair 
(1681,1682). 


Proof 2. Let T, = k? be a square triangular number, and 
consider the pair (87,,, 87, +1). In terms of n and K, this is the pair 
(8k2,(2n + 1)*) using Lemma B.4| If a prime p divides 8k? or 
(2n + 1)’, so does p’, and the result follows since there are infi- 
nitely many square triangular numbers. However, not all pairs 
of consecutive powerful numbers are obtained this way, e.g., 
(675, 676) = (33 - 52,22 - 13”). 


Let (Xn, Yn) be one of the infinitely many solutions in integers to the 
Pell equation x* — 2y% = 1 in Table Then x2 - 1 
(= y?2 +2), x2 +0, x2 +1 satisfy the requirements of the problem. 
Yes, infinitely many. The equation (y — 1)°+y2+(y + 1)* = x?41 
simplifies to x? — 3y* = 1. See Sere 

Yes, infinitely many. The equation T,, = 37; is equivalent to the 
Pell equation x? — 3y? = —2 withx = 2m+1andy = 2k +1. 
Hence the pairs (Xn, Yn) in Table B.élyield solutions, e.g., (x3, y3) = 
(71,41) yields (m,k) = (35,20), (%4,y4) = (265,153) yields 
(m, k) = (132, 76), etc. 

The kth octagonal number Oct, is k? + 47, = k(3k — 2), and 
the equation k (3k — 2) = m? yields the Pell equation (3k — 1)° — 
3m? = 1. So solutions (x, y,) in Table B.5| with Xn = —1 (mod 3) 
yield square octagonal numbers, e.g., (x3, y3) = (26,15) yields 
Octy = 157, (xs, ys) = (362, 209) yields Oct,, = 209”, etc. 

Yes, infinitely many. The equation S,,, = k* is equivalent to the Pell 
equation (2k)°-6(2m — 1) = —2. Theseed (xg, yo) = (2, 1) with 
x2 — 6y?% = —2 generates infinitely many solutions, e.g., (x1, y,) = 
(22,9) yields S; = 117, (x3,y2) = (218,89) yields S4;= 1097, 
(x3, v3) = (2158, 881) yields S44, = 1079”, etc. 
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3.12 Tm = k? — 1is equivalent to (2m + 1)* — 2(2k)* = —7. The seed 
(xo, ¥o) = (1,2) with x2 — 2y2 = —7 generates infinitely many 
solutions which lead to T; = 42 — 1, T37 = 23° — 1, etc.; while 
the seed (x9, 9) = (—1,2) leads to T> = 2? —1, Ts = 11° -1, 
etc., Tm = k2 + 1 is equivalent to (2m + 1)* — 2(2k)*? = 9. The 
seed (Xo, Yo) = (3, 0) with x2 — 2y? = 9 generates infinitely many 
solutions which lead to T, = 32+1,T 5 = 187 +1, Ty4g = 105° +1, 
etc. 

3.13 Yes. Consider the equation m? = (k + 1)* —k3 = 3k? + 3k +1. 
Completing the square yields (2m)* — 3(2k + 1: = 1. So pairs 
(Xn, Yn) in Table with x, even and y, odd yield solutions, e.g., 
(x3, ¥3) = (26,15) yields (m,k) = (13,7), (5, ¥5) = (362, 209) 
yields (m,k) = (181, 104) so that 1817 = 105° — 104? etc. 

3.14 Ify? = k+1andx? = 3k + 1, then x? — 3y* = —2, and solu- 
tions (Xn, Y,) to this Pell equation appear in Table Then k = 
(x2 — y,2)/2 leads to the solutions k = 8, 120, 1680,.... 

3.15 Yes. Tm-1 + Tin + Tma1 = k2 is equivalent to (4k)? — 6(2m + 1)” 
= 10. Solving the Pell equation x* — 6y? = 10 with initial solution 
(xo, Yo) = (4, 1) yields solutions T3+Te4+T6s = 79”, Tea7 +Teaa + 
Tg39 = 782°, etc. 

3.16 See Figure for a visual proof that K = rs, using the fact that 
the three angle bisectors in a triangle meet at the center of the in- 
scribed circle. 


FIGURE S3.1 


For an almost equilateral Heronian triangle in Example 
K = 3Xnyy and s = 3xXy, hence r = yp, an integer. 


3.17 The probability that both balls are red is 1/2 = (7)/(), so that 


2r(r—1) = k(k — 1). Completing the square yields (2k — 1)’ 
~2(2r — 1)* = -1, the Pell equation x?—2y? = —1 withx = 2k- 
1 and y = 2r — 1. Each solution (%y, y,) forn = 1 in Table 8.4) has 
x and y odd and y > 2, and hence leads to a solution. For example, 
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(x1, ¥1) = (7,5) leads to (k,r) = (4,3), (x2, y2) = (41, 29) leads to 
(k,r) = (21,15), (x3, y3) = (239, 169) leads to (k,r) = (120,85), 
etc. 


Chapter i] 


If a and b are both even, then so is c, so (a,b,c) cannot be primi- 
tive. If a and b are both odd, then a® = 1 (mod 4), 
b* = 1 (mod 4), and soc? = 2 (mod 4), which is impossible since 
even squares are congruent to 0 (mod 4). Hence a and b have op- 
posite parity, and c is odd. 

2mn is a multiple of 4 since either m or n is even. Now consider 
the generator (m,n) = (2k,1) for k = 1. 

When c = m?+n?,2c = 2(m? +n?) = (m+ n)’+(m—n)’,and 
both m + n and m — nare odd since m and n have opposite parity. 
Yes, infinitely many. The generator (m,n) = (Tx, Ty-1) fork = 2 
generates the PT (k?, T,2_1,7,2). None are primitive, since each 
side is a multiple of k. 


The generator (m,n) = Ce 1) for k = O generates the PPT 


Fret1 — 2,2(f; — 1), fe+1)- 

Tm = k(k + 1) implies m2? + (m+ 1)? = 2m? + 2m+1 = 4k? + 
4k +1 = (2k +1)’, and conversely m2 + (m+1)* = (2k +1)° 
implies m (m + 1) = 2k(k + 1) so that T,, = k(k + 1). 

The generators (m,n) = (p+1,1) and (m,n) = (p +1, p) generate 
different PPTs with inradius p. 

The generator (m, 1) generates a PPT with area 


(m — 1) m(m + 1). 


Yes, infinitely many. The PPT (2k + 1,2k? + 2k,2k? + 2k +1) 
from Example 4.1)has P = 4k? + 6k + 2 = 2T x4 and 


A=k(k+1)(2k +1) = 6Pyr,. 
Express the inradius and exradii in terms of the generator 
(m,n) and simplify. 
Using the hint we have 
(ab)* + (bc)* + (ac)* = atb* + c* (a* + b*) 
=atb*+c*(a2 + b?)" — 2c*ab? 
=e o2Cu bap 


= (c* _ a*b?)?, 
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Setting g = 0 in Theorem 4.19] yields 

(m2 +n? + p*)” = (m2 + n2 — p2)” + (2mp)” + (2np)’. 
Note: Adrien-Marie Legendre (1752-1833) proved that an integer 
is the sum of three positive integer squares if and only if it is not of 
the form 4"(8k + 7) for non-negative integers k and n. 
There are infinitely many such sets of three squares. The sides of 
the three squares in Figure 4.24] are b-—a,c,and a+b, and hence 
their areas (b — a)’ = c? — 2ab, c2, and (a+b)° = c? + 2ab 
form an arithmetic progression of three squares with a common 
difference equal to four times the area of the PA. The two examples 
in the statement of the problem are generated by the PA’s (3,4,5) 
and (7,24,25). 
Yes, infinitely many. From Lemma 4.21] we have 


(a2 + b?)? = (a2 — b?)? + (2ab)?, 
so that 
(a2 +b?) = [(a2 — 2)? + (2ab)?] (a2 + b”) 
= (a? + ab)? + (a*b + b?)? 
= (a? — 3ab’)? + (3a*b — b)?. 
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Using the hint we have a? + b* = 0 (mod 3), hencea = b = 
0 (mod 3) from Example This contradicts a and b relatively 
prime. 

(a) The Gelfand-Schneider number 2V2 is an example. 

(b) 1%? = Lisan example. As a consequence of the solution to 
Hilbert’s seventh problem, we cannot provide an example showing 
that a rational not equal to 0 or 1 to an irrational power may be 
rational. 

Assume vn? — 1 is irrational, and write Vn? — 1 = a/b in lowest 
terms for positive integers a and b. Then a2 + b? = (nb)’. Follow- 
ing the hint we modify Figure so that the large right triangle 
has hypotenuse nb, vertical leg b, and horizontal leg a. Conse- 
quently, the smaller gray triangle has legs nb — a and 
na — (n* — 1)b, and hypotenuse n*b — na. Thus Vn2-—1 = 
[na — (n? — 1)b]/(nb — a), a contradiction since na — (n? — 1)b 
and nb — a are smaller than a and D, respectively. Hence Vn? — 1 
is irrational. 
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Assume vn? + 1 is irrational, and write Vn? + 1 = a/b in lowest 
terms for positive integers a and b. Then b? + (nb)* = a2. Follow- 
ing the hint we modify Figure 5.8)so that the large right triangle has 
hypotenuse a, vertical leg b, and horizontal leg nb. Consequently 
the smaller gray triangle has legs a — nb and na — n?b, and hy- 
potenuse (n° +1)b—na. Thus vn? — 1 = [(n?+1)b—na]/(a—nb), 
a contradiction since (n? + 1)b — na and a — nb are smaller than 
a and b, respectively. Hence Vn? + 1 is irrational. 

(a) Create a square with area 3 in Figure 6.12|by drawing a diagonal 
in each of the four rectangles. 

(b) Partition the square with side length g+1+¢@ into a unit square 
surrounded by four g X (@ + 1) rectangles. Hence (2g + 1) = 
1+4(p? + ¢). But gp? = yg? +9 = 29 +150 that (p?)? = 1+ 49? 
from which the result follows. 


(a) Since o, = $n + Vn? +4), np, = +(n?+Vnt Fan) = 
[n?, 1]. Thus g2 =1+ng, =1+ [n2, 1] a [1 +n, Ln?|, 


n 


(b) Expanding and simplifying (on = =) =o yields os _ 4 = 


n? + 3n, hence 03 = @y343n = [n3 + 3n|. 


Yes. Itis 
2- 2+ J2t--+v2 


with n square root signs. This follows from sin@ = v1 — cos? 0 
with @ = —— 


2nt1* 


Using the hint we have x = ja+Vva+Vat- = Ja+x,sox? — 


x—a=0. Hence x = =[ + v1+ 4a]. Thus 


l2+j24+V2+~ = 2and 1+j14+vie~=o. 
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6.1 See Figure 
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FIGURE S6.1 


6.2 See Figure 


Fret 


FIGURE S6.2 


6.3 The identities in (6.12) yield 5F, = 2Ln41 —Lyn = Ln41 +Ly_1 and 


Ly = 2Fn4a — Fa = Fraga + Pot: 


6.4 (a) Liat LoR-1 = Liat (Lox = LoK-2) = Lon — Lo = Lon — 2. 
(b) ae Lox = eH1 (Lox41 = Lox-1) =Lons1 — Ly =Lony1 — 1. 


6.5 (a) Solve (b) first, then set k = n. 


(b) 


Fighin + Fry = Fe Fn-1 + Fra) + Fr Ce-a + Feta) 
= (Fi Fao + Fe—1Fn) + nFicta + Fr-1 Fe) 
= Fratk + Frek = 2Fr+k: 

6.6 Relabel Figure with Lucas numbers and follow the steps 
in the proof of Theorem 6.5] to reach Ly4yqoggy + Ln ly = Lelnset 
Lylnte-1 = 3Ln+K + Lntk-1 = 2Ln+K + Lntkes = Lntk + Lnskta2 
and note that Lyi + Lnskse = SFysrs. (from Exercise 6.3). 


6.7 
6.8 


6.9 


6.10 


7.1 


7.2 
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FIGURE S6.3 


See Figure (Brousseau, 1972] . 


In the quadratic equation x? — bx + c = 0, the coefficient b is the 
sum of the roots, and the coefficient c is the product. Thus a qua- 
dratic equation with roots p” and (—~) "isx?—L,x+(—1)" = 0. 
Since F3;,, is even, both F3,_, and F3,,, are odd. Hence L3;, = F3,_4+ 
F344 is even (see Exercise 6.3). 

(a) See the second identity in raDN 

(b) The harmonic mean of a and b is 2ab/(a + b), so from 
and Exercise 6.5(a) the harmonic mean of F, and L,, is 2Fy,/2Fy+1.- 


Chapter 


From the proof of Theorem 7.4] we have Ny = T3n+41 and it is easily 
verified that T3441 = Tons1 + Ton + Ty. See Figure for a visual 


proof. 


FIGURE S7.1 


The pattern is N, = Top_, = 3T2p-1_1 + T2p-1. So we need only 
show that T>,_, = 37, 1, + T, (which can be verified by simple 
algebra) and setn = 2?71. In Figure we illustrate T,,_; = 
3Tn-1 + Ty for n = 12, i.e., T>3 partitioned into three copies of T,, 
(in light gray) and one T;, (in dark gray). 
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7.4 
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FIGURE S7.2 


For p = 2k + 1 the pattern is Ny = Ty2n+1_, = Tok_4Tok+1 + 
T,xT>k+1_1. SO we need only show that Tyyn—1 = Tm—1Tn + TmTn-1 
(which can be verified by simple algebra) and set m = 2%,n = 
2**1. In Figure S7.3 we illustrate Tmn—1 = Tm—1Tn + TmTn—1 for 
(m,n) = (4,6), ie., T23 partitioned into T3 copies of T, (in dark 
gray) and T, copies of T; (in light gray). 


FIGURE S7.3 


From the proof of Theorem [7.4] N, = 1 (mod 3), and N, = 0 
(mod 4) since N,, is divisible by 4. Solving these two pone ences 
simultaneously (asi in Example B.3) ) yields N, =4 (mod 12). 

For p = 5we have N, = 0 (mod 16) andN,, = 1 (mod 9). Solving 
these two sonemuences simultaneously is in Example yields 
N, = 64 (mod 144) so that N, = 144k + 64 for some positive 
integer k. Hence in base 12 the rightmost two digits of Np are 54, 
the base 12 number equal to 64 in base 10. 

From the representation of N,, as the sum of a geometric progres- 
sion in Section we have Ny = 111:--- 100--- Opase 2, with p ones 
followed by p — 1 zeros. ; 


J1+8N, = 2?*1 —1 since 1+8- 2?-1 (2? — 1) = (2+? w 


For a visual proof, recall that N,, is triangular and see Figure 
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7.8 When p = 3, Nz is a multiple of 4. Hence when the units digit is 
6, the rightmost two digits must be 16, 36, 56, 76, or 96. All five 
possibilities are known to occur. 

7.9 (a) It suffices to show that 


2? = 2 (mod 6) and 2” = 8 (mod 24). 


As in the proof of Theorem we have 2? = (a1)" = -1 
2 (mod 3). Using the Chinese remainder theorem to solve 2” 
0 (mod 2) and 2? = 2 (mod 3) yields 2? = 2 (mod 6); simi- 
larly solving 2? = 0 (mod 8) and 2? = 2 (mod 3) yields 2? = 8 
(mod 24). 

(b) Ifp = 4k +1, then 2? = 2-42* =2-(-1)*" = 2 (mod 5), so 
that 2” ends in 2 or 7 (base 10). Since 2” is even it ends in 2, and 
hence M, = 2? —1endsin 1. Ifp = 4k + 3, then 2? = 8- 42k = 


3: (eres = 3 (mod 5), so that 2? ends in 3 or 8 (base 10). Since 
2? is even it ends in 8, and hence M, = 2? — 1 ends in 7. 
7.10 No, since 2? — 1 #1 (mod 4). 
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